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According to recent investigations by F. H. Seares,' there is considerable 
evidence in favor of the idea that the center of the larger galactic system 
is identical with the center of the system of globular clusters, placed by 
H. Shapley? approximately in galactic longitude 325° and at a distance of 
approximately 20,000 parsecs from the sun. J. H. Oort’s center of galactic 
rotation lies in the same direction, although at a somewhat smaller distance. 
His work has been amply confirmed by J. Schilt, J. S. Plaskett and others, 
so that it is safe to regard the direction ] = 325° as fundamental in our 
Milky Way system. B. Lindblad and others have advanced the hypothesis 
that the observed distribution of stars is mainly an effect of dark nebulous 
masses in the vicinity of the solar system, and that the actual center is 
hidden from us. 

H. Shapley and his collaborators* have undertaken a detailed study of 
the region near the galactic center. It appears that there is a notable 
concentration of cluster-type variable stars in the vast star-clouds near the 
center, at a distance of about 14,400 parsecs. Many variables lie beyond 
the center, and there are a number of extra-galactic nebulae within a small 
angular distance from the point 1 = 325°, b = 0° (or a = 17°95; 6 = 
—30°). Accordingly, in comparatively small angular distances from the 
center, the line of sight encounters no obstruction from dark matter. If 
we consider a series of observations covering the whole galaxy and com- 
plete to a certain given apparent magnitude, the average distance of the 
stars must be a function of galactic longitude, being greatest in] = 325° 
and least in] = 145°. The exact form of the relationship depends upon 
the distribution of dark matter as well as upon the concentration of the 
stars. toward the center, and is therefore not known, but the curve is 
probably symmetrical and its maximum and minimum are defined as 
stated. 

The determinations of the Ca+ intensities in early type stars, which I 
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made at the Harvard College Observatory,‘ are fairly uniform with respect 
to limiting magnitude. It might therefore be expected that the mean 
intensity of an interstellar line would vary with galactic longitude. The 
intensity 

I=] 7h” 


where @ is the absorption coefficient per parsec and D is the distance in 
parsecs. What we observe is 


Am = —2.5 log I/I) = 2.5 mod. B.D 


If B is constant, the observed value of Am is directly proportional to the 
distance. Stronger lines should therefore predominate near the galactic 
center. In reality we do not know whether or not @ is quite constant 
over large distances. It is probably fairly constant up to D = 1000 
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Relationship between galactic longitude and average intensity of 
interstellar calcium lines. 


parsecs, but seems to decrease for greater distances from the sun. Whether 
or not it is a function of galactic longitude cannot be determined from the 
observations. A. S. Eddington’ has shown that the degree of ionization 
of the interstellar Ca-gas is that of a gas of density pé in thermodynamic 
equilibrium, where p is the density of the interstellar gas and 6 is the 
“dilution-factor’’ determined in such a way that the density of radiation 
between y and y + dy in interstellar space is 1/5 ](yi:7)dy. The density 
of radiation increases as the galactic center is approached, so that 6 de- 
creases. But the density of interstellar calcium probably increases—at 
least this is the most natural hypothesis. Indeed, if the interstellar gas is 
directly due to expulsion from stellar.atmospheres by radiation pressure, 
we should have very nearly pS = const., and 8 would be independent of 
galactic longitude. In any case it is extremely unlikely that a change in 
8 would exactly offset the variation in mean distance. 
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The observations furnish direct evidence that 8 is really a function of 
galactic longitude. In Astrophysical Journal 67, 370, 1928, figure 2, I 
have given a graphical representation of the Ca’ intensities in four 
shadings: 


Very strong - Am 2 3.7 
Strong 3.2 < Am < 3.7 
Weak 2.8 < Am < 3.2 
Very weak Am < 2.8 


For convenience of computation I have introduced the following new units 


of intensity: 
INTENSITY IN NEW 


CHARACTER OF LINES UNITS 
Very strong +2 
Strong +1 
Weak = 4 
Very weak —2 


One unit in the new scale corresponds roughly to a step of three or four 
hundredths of a stellar magnitude. 

Means for every 30° in J, between b = +10° and b = —10°, were 
then taken. They are shown in table 1, expressed in the new units. 
Every entry of the table corresponds to the mean from six squares each 
having 10° on the side, as shown in the figure referred to above. Only in 


TABLE 1 
NO. OF SQUARES MEAN INTENSITY IN 
MEAN GAL. LONG. SED NEW UNITS 
15 4 +0.25 
45 6 +0.33 
75 6 +0.83 
105 6 —0.33 
135 5 —1.40 
165 4 +0.75 
195 6 —0.33 
225 6 —0.17 
255 6 +0.50 
285 6 +0.33 
315 6 +0.83 
345 4 +1.75 


a few cases there were blank squares in the figure, due to an insufficient 
number of stars in these regions. In such cases the number of squares 
used was smaller. It seemed best not to use any of the squares in high 
galactic latitudes because of their small number and unsymmetrical dis- 
tribution. There is also a tendency for the lines to be fainter in the 
higher galactic latitudes. 

This series of values can be represented by an equation of the form 











166 ASTRONOMY: 0. STRUVE Proc. N. A. S. 


Am = K + X cos (I — a) (1) 


where a is the direction of maximum intensity. A graphical solution 
gives for the unknowns 


a = 325° 
K = +0.2 units 
xX = 0.6 ‘i 


These values can be improved by a least-squares solution. Differentiat- 
ing (1) we obtain for our observational equations 


d(Am) = dK + cos (l — a)dX + X sin (1 — a)da. (2) 
Considering as unknowns the corrections to the preliminary values, 


Aa, AK, AX, I have made a rough computation which yields the follow- 
ing final results: 


a = 337.6° + 18° (probable error) 
K = +0.28 + 0.15 < zs 
X = 0.67 + 0.20 ty - 


The direction a very nearly coincides with the direction toward the 
galactic center. The results are shown in figure 1 in which the curve is 
drawn for ] = 325°. In the computations equal weights were given to 
all observations. If the three points determined from only four squares 
each are omitted, we get a slightly different curve, as shown by the dashed 
lines. The value of a would not, however, be appreciably affected. The 
least-squares solution reduced the sum of the squares of the residuals from 
429 to 414. The probable error of one point on the diagram is + 0.5 in 
the new units. 

While it is possible that the results of the intensity estimates are affected 
by systematic errors, I do not think that they are serious in this case. It 
so happens that the position of the center is only 30° distant from the 
equator, so that systematic differences between northern and southern 
hemispheres would not much affect the results, and these are probably 
the only ones that might be feared. 


1 Astrophys. J., 67, 123, 1928. 

2 Tbid., J., 49, 311, 1919. 

*H. Shapley, Nature, 122, 482, 1928; H. Shapley, these ProcEEDINGs, 14, 825, 
1928; H. Shapley and H. H. Swope, Jbid., 14, 830, 1928. 

4 Astrophys. J., 67, 353, 1928. ; 

5 Der Innere Aufbau der Sterne, p. 478, 1928. 
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THE LARGE RADIAL VELOCITY OF N. G. C. 7619 


By Miron L. HUMASON 
Mount WILSON OBSERVATORY, CARNEGIE INSTITUTION OF WASHINGTON 


Communicated January 17, 1929 


Centered around the region at R. A. 23" 16", Dec. +7° 50’ (1920) is 
a cluster of small spiral, globular, and elongated nebulae. Two of the 
brightest nebulae in this group are N. G. C. 7619 and 7626. The first 
of these is a globular nebula and at the Cassegrain focus of the 100-inch 
telescope appears visually to be a little brighter than N. G. C. 7626. 

About a year ago Mr. Hubble suggested that a selected list of fainter 
and more distant extra-galactic nebulae, especially those occurring in 
groups, be observed to determine, if possible, whether the absorption lines 
in these objects show large displacements toward longer wave-lengths, 
as might be expected on de Sitter’s theory of curved space-time. 

During the past year two spectrograms of N. G. C. 7619 were obtained 
with Cassegrain spectrograph VI attached to the 100-inch telescope. 
This spectrograph has a 24-inch collimating lens, two prisms, and a 3-inch 
camera, and gives a dispersion of 183 A per millimeter at \ 4500. The 
exposure times for the spectrograms were 33" and 45", respectively 
The radial velocity from these plates has been measured by Miss Mac- 
Cormack, of the computing division, and by myself, the weighted mean 
value being +3779 km./sec. The velocity of this nebula is, therefore, 
twice as large as any hitherto observed, the highest previously known 
‘being that of N. G. C. 584, for which Slipher obtained +1800 km./sec. 
Individual velocities from the two plates are +3828 km./sec. for the shorter 
exposure, and +3754 km./sec. for the longer, which is much the better of 
the two exposures and is given double weight. 

It may be mentioned that Hubble, in a paper in these PROCEEDINGS, 
gives approximate distances for 24 extra-galactic nebulae, and finds a 
marked increase in radial velocity with distance. The high velocity for 
N. G. C. 7619 derived from these plates falls on the extrapolated line which 
expresses the relationship between line displacement and distance. These 
results suggest an influence of distance upon the observed line shift—such 
as would be produced, for example, on de Sitter’s theory, both by the 
apparent slowing-down of light vibrations with distance and by a real 
tendency of material bodies to scatter in space. 

The spectral type of N. G. C. 7619 is estimated as F8, and the probable 
error of the measured velocity is presumably not greater than 100 km./sec. 
The large probable error is due to the very poor quality of the absorption 
lines, but is not of great importance when dealing with such large displace- 
ments. ‘The lines are not well defined, being rather wide and diffuse. In 
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appearance the spectrum is very much like spectra of the Milky Way 
clouds in Sagittarius and Cygnus, and is also similar to spectra of binary 
stars of the W Ursae Majoris type, where the widening and depth of the 
lines are affected by the rapid rotation of the stars involved. 

The wide shallow absorption lines observed in the spectrum of N. G. C. 
7619 have been noticed in the spectra of other extra-galactic nebulae, and 
may be due to a dispersion in velocity and a blending of the spectral types 
of the many stars which presumably exist in the central parts of these 
nebulae. The lack of depth in the absorption lines seems to be more 
pronounced among the smaller and fainter nebulae, and in N. G. C. 7619 
the absorption is very weak. 

It is hoped that velocities of more of these interesting objects will soon 
be available. 


A RELATION BETWEEN DISTANCE AND RADIAL VELOCITY 
AMONG EXTRA-GALACTIC NEBULAE 


By Epwin HuBBLE 
Mount WILson OBSERVATORY, CARNEGIE INSTITUTION OF WASHINGTON 


Communicated January 17, 1929 


Determinations of the motion of the sun with respect to the extra- 
galactic nebulae have involved a K term of several hundred kilometers 
which appears to be variable. Explanations of this paradox have been 
sought in a correlation between apparent radial velocities and distances, 
but so far the results have not been convincing. The present paper is a 
re-examination of the question, based on only those nebular distances 
which are believed to be fairly reliable. 

Distances of extra-galactic nebulae depend ultimately upon the appli- 
cation of absolute-luminosity criteria to involved stars whose types can 
be recognized. These include, among others, Cepheid variables, novae, 
and blue stars involved in emission nebulosity. Numerical values depend 
upon the zero point of the period-luminosity relation among Cepheids, 
the other criteria merely check the order of the distances. This method 
is restricted to the few nebulae which are well resolved by existing instru- 
ments. A study of these nebulae, together with those in which any stars 
at all can be recognized, indicates the probability of an approximately 
uniform upper limit to the absolute luminosity of stars, in the late-type 
spirals and irregular nebulae at least, of the order of M (photographic) = 
—6.3.!_ The apparent luminosities of the brightest stars in such nebulae 
are thus criteria which, although rough and to be applied with caution, 
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furnish reasonable estimates of the distances of all extra-galactic systems 
in which even a few stars can be detected. 


TABLE 1 
NEBULAE WHOSE DisTANCES HAVE BEEN ESTIMATED FROM STARS INVOLVED OR FROM 
MEAN LUMINOSITIES IN A CLUSTER 





OBJECT ™m, r v m, M, 
S. Mag. aie 0.032 + 170 1.5 —16.0 
L. Mag. “ 0.034 + 290 0.5 17.2 

N. G. C. 6822 se 0.214 — 130 9.0 12.7 
598 ws 0.263 —- 7 7.0 15.1 
221 ee 0.275 — 185 8.8 13.4 
224 res 0.275 — 220 5.0 17.2 
5457 17.0 0.45 + 200 9.9 13.3 
4736 17.3 0.5 + 290 8.4 15.1 
5194 17.3 0.5 + 270 7.4 16.1 
4449 17.8 0.63 + 200 9.5 14.5 
4214 18.3 0.8 + 300 11.3 13.2 
3031 18.5 0.9 =— 3 8.3 16.4 
3627 18.5 0.9 + 650 9.1 15.7 
4826 18.5 0.9 + 150 9.0 15,7 
5236 18.5 0.9 + 500 10.4 14.4 
1068 18.7 1.0 + 920 9.1 15.9 
5055 19.0 Et + 450 9.6 15.6 
7331 19.0 £3 + 500 10.4 14.8 
4258 19.5 1.4 + 500 8.7 17.0 
4151 20.0 1 af + 960 12.0 14.2 
4382 2.0 + 500 10.0 16.5 
4472 2.0 + 850 8.8 by £4 
4486 2.0 + 800 9.7 16.8 
4649 we 2.0 +1090 9.5 17.0 
Mean —15.5 


m, = photographic magnitude of brightest stars involved. 

distance in units of 10° parsecs. The first two are Shapley’s values. 

measured velocities in km./sec. N. G. C. 6822, 221, 224 and 5457 are recent 

determinations by Humason. 

m: = Holetschek’s visual magnitude as corrected by Hopmann. The first three 
objects were not measured by Holetschek, and the values of m: represent 
estimates by the author based upon such data as are available. 

M, = total visual absolute magnitude computed fromm: and r. 


en 
i il 


Finally, the nebulae themselves appear to be of a definite order of 
absolute luminosity, exhibiting a range of four or five magnitudes about 
an average value M (visual) = —15.2.! The application of this statistical 
average to individual cases can rarely be used to advantage, but where 
considerable numbers are involved, and especially in the various clusters 
of nebulae, mean apparent luminosities of the nebulae themselves offer 
reliable estimates of the mean distances. 

Radial velocities of 46 extra-galactic nebulae are now available, but 
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individual distances are estimated for only 24. For one other, N. G. C. 
3521, an estimate could probably be made, but no photographs are avail- 
able at Mount Wilson. The data are given in table 1. The first seven 
distances are the most reliable, depending, except for M 32 the companion of 
M 31, upon extensive investigations of many stars involved. The next 
thirteen distances, depending upon the criterion of a uniform upper limit 
of stellar luminosity, are subject to considerable probable errors but are 
believed to be the most reasonable values at present available. The last 
four objects appear to be in the Virgo Cluster. The distance assigned 
to the cluster, 2 X 10° parsecs, is derived from the distribution of nebular 
luminosities, together with luminosities of stars in some of the later-type 
spirals, and differs somewhat from the Harvard estimate of ten million 
light years.” 

The data in the table indicate a linear correlation between distances and 
velocities, whether the latter are used directly or corrected for solar motion, 
according to the older solutions. ‘This suggests a new solution for the solar 
motion in which the distances are introduced as coefficients of the K term, 
i. e., the velocities are assumed to vary directly with the distances, and 
hence K represents the velocity at unit distance due to this effect. The 
equations of condition then take the form 


rK +X cosa cosi+ Y sinacos6 + Z sind = v. 


Two solutions have been made, one using the 24 nebulae individually, 
the other combining them into 9 groups according to proximity in direc- 
tion and in distance. ‘The results are 


24 OBJECTS 9 GROUPS 

x — 65 = 50 + 3+ 70 

Y +226 + 95 +230 + 120 

Z —195 + 40 —133 + 70 

K +465 + 50 +513 + 60km./sec. per 10 parsecs. 
A 286° 269° 

D + 40° + 33° 

Vo 306 km./sec. 247 km./sec. 


For such scanty material, so poorly distributed, the results are fairly 
definite. Differences between the two solutions are due largely to the 
four Virgo nebulae, which, being the most distant objects and all sharing 
the peculiar motion of the cluster, unduly influence the value of K and 
hence of Vo. New data on imore distant objects will be required to reduce 
the effect of such peculiar motion. . Meanwhile round numbers, inter- 
mediate between the two solutions, will represent the probable order of 
the values. For instance, let A = 277°, D = +36° (Gal. long. = 32°, 
lat. = +18°), Vo = 280 km./sec., K = +500 km./sec. per million par- 
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secs. Mr. Strémberg has very kindly checked the general order of these 
values by independent solutions for different groupings of the data. 

A constant term, introduced into the equations, was found to be small 
and negative. ‘This seems to dispose of the necessity for the old constant 
K term. Solutions of this sort have been published by Lundmark,? who 
replaced the old K by k + /r + mr*. His favored solution gave k = 513, 
as against the former value of the order of 700, and hence offered little 
advantage. 


TABLE 2 
NEBULAE. WHOSE DisTANCES ARE ESTIMATED FROM RADIAL VELOCITIES 

OBJECT 9 %, r m, M, 
N.G.C. 278 + 650 —110 1.52 12.0 —13.9 

404 — 25 — 65 os 11.1 ea 
584 +1800 + 75 3.45 10.9 16.8 
936 +1300 +115 2.37 bE ae 15.7 
1023 + 300 -— 10 0.62 10.2 13.8 
1700 + 800 +220 1.16 12.5 12.8 
2681 + 700 -— 10 1.42 10.7 15.0 
2683 + 400 + 65 0.67 9.9 14.3 
2841 + 600 — 20 1.24 9.4 16.1 
3034 + 290 —105 0.79 9.0 15.5 
3115 + 600 +105 1.00 9.5 15.5 
3368 + 940 + 70 1.74 10.0 16.2 
3379 + 810 + 65 1.49 9.4 16.4 
3489 + 600 + 50 1.10 11.2 14.0 
3521 + 730 + 95 1.27 10.1 15.4 
3623 + 800 + 35 1.53 9.9 16.0 
4111 + 800 — 95 1.79 10.1 16.1 
4526 + 580 — 20 1.20 11.1 14.3 
4565 +1100 — 75 2.35 11.0 15.9 
4594 +1140 + 25 2.23 9.1 17.6 
5005 + 900 —130 2.06 12% 15.5 
5866 + 650 —215 1.73 1 es —14.5 
Mean 10.5 —15.3 


The residuals for the two solutions given above average 150 and 110 
km./sec. and should represent the average peculiar motions of the in- 
dividual nebulae and of the groups, respectively. In order to exhibit 
the results in a graphical form, the solar motion has been eliminated from 
the observed velocities and the remainders, the distance terms plus the 
residuals, have been plotted against the distances. The run of the re- 
siduals is about as smooth as can be expected, and in general the form of 
the solutions appears to be adequate. 

The 22 nebulae for which distances are not available can be treated in 
two ways. First, the mean distance of the group derived from the mean 
apparent magnitudes can be compared with the mean of the velocities 
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corrected for solar motion. The result, 745 km./sec. for a distance of 
1.4 X 10° parsecs, falls between the two previous solutions and indicates 
a value for K of 530 as against the proposed value, 500 km./sec. 
Secondly, the scatter of the individual nebulae can be examined by 
assuming the relation between distances and velocities as previously 
determined. Distances can then be calculated from the velocities cor- 
rected for solar motion, and absolute magnitudes can be derived from the 
apparent magnitudes. The results are given in table 2 and may be 
compared with the distribution of absolute magnitudes among the nebulae 
in table 1, whose distances are derived from other criteria. N.G. C. 404 
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Velocity-Distance Relation among Extra-Galactic Nebulae. 

Radial velocities, corrected for solar motion, are plotted against 
distances estimated from involved stars and mean luminosities of 
nebulae in a cluster. The black discs and full line represent the 
solution for solar motion using the nebulae individually; the circles 
and broken line represent the solution combining the nebulae into 
groups; the cross represents the mean velocity corresponding to 
the mean distance of 22 nebulae whose distances could not be esti- 
mated individually. 


can be excluded, since the observed velocity is so small that the peculiar 
motion must be large in comparison with the distance effect. The object 
is not necessarily an exception, however, since a distance can be assigned 
for which the peculiar motion and the absolute magnitude are both within 
the range previously determined. The two mean magnitudes, —15.3 
and —15.5, the ranges, 4.9 and 5.0 mag., and the frequency distributions 
are closely similar for these two entirely independent sets of data; and 
even the slight difference in mean magnitudes can be attributed to the 
selected, very bright, nebulae in the Virgo Cluster. This entirely unforced 
agreement supports the validity of the velocity-distance relation in a very 
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evident matter. Finally, it is worth recording that the frequency distribu- 
tion of absolute magnitudes in the two tables combined is comparable 
with those found in the various clusters of nebulae. 

The results establish a roughly linear relation between velocities and 
distances among nebulae for which velocities have been previously pub- 
lished, and the relation appears to dominate the distribution of velocities. 
In order to investigate the matter on a much larger scale, Mr. Humason 
at Mount Wilson has initiated a program of determining velocities 
of the most distant nebulae that can be observed with confidence. 
These, naturally, are the brightest nebulae in clusters of nebulae. 
The first definite result,4 1 = + 3779 km./sec. for N. G. C. 7619, is 
thoroughly consistent with the present conclusions. Corrected for the 
solar motion, this velocity is +3910, which, with K = 500, corresponds to 
a distance of 7.8 X 10° parsecs. Since the apparent magnitude is 11.8, 
the absolute magnitude at such a distance is —17.65, which is of the 
right order for the brightest nebulae in a cluster. A preliminary dis- 
tance, derived independently from the cluster of which this nebula appears 
to be a member, is of the order of 7 X 10° parsecs. 

New data to be expected in the near future may modify the significance 
of the present investigation or, if confirmatory, will lead to a solution 
having many times the weight. For this reason it is thought premature 
to discuss in detail the obvious consequences of the present results. For 
example, if the solar motion with respect to the clusters represents the 
rotation of the galactic system, this motion could be subtracted from the 
results for the nebulae and the remainder would represent the motion of 
the galactic system with respect to the extra-galactic nebulae. 

The outstanding feature, however, is the possibility that the velocity- 
distance relation may represent the de Sitter effect, and hence that numer- 
ical data may be introduced into discussions of the general curvature of 
space. In the de Sitter cosmology, displacements of the spectra arise 
from two sources, an apparent slowing down of atomic vibrations and a 
general tendency of material particles to scatter. The latter involves an 
acceleration and hence introduces the element of time. The relative im- 
portance of these two effects should determine the form of the relation 
between distances and observed velocities; and in this connection it may 
be emphasized that the linear relation found in the present discussion is a 
first approximation representing a restricted range in distance. 

1 Mt. Wilson Contr., No. 324; Astroph. J., Chicago, IIl., 64, 1926 (321). 

2 Harvard Coll. Obs. Circ., 294, 1926. 

3 Mon. Not. R. Astr. Soc., 85, 1925 (865-894). 

‘These PRocEEDINGS, 15, 1929 (167). 
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STUDIES OF THE GALACTIC CENTER. IV. ON THE 
TRANSPARENCY OF THE GALACTIC STAR CLOUDS 


By HarLow SHAPLEY 
HARVARD OBSERVATORY, CamprincE, MASSACHUSETTS 
Communicated February 14, 1929 
1. Preceding papers of this series contain reports on the progress of an 
investigation of the dimensions and structure of the Galaxy. In the 
first, a program was outlined for a fifteen-year study of the distant variable 
stars in the Milky Way. The paper incidentally reported the discovery 


during the past five years of more than a thousand faint variables, most 
of which are more than thirty thousand light years distant. 


‘ 





FIGURE 1 


Composite picture of the central Milky Way region, made from plates taken at 
Arequipa with a three-inch Ross-Tessar lens. The rings around the brighter stars 
are merely optical. The center of the Galaxy is indicated at 327°, 0°. This 
picture covers a somewhat larger area than figures 2 and 3. 


In the second paper was presented the evidence bearing on the existence 
and dimensions of a massive galactic nucleus. The star clouds in the 
Sagittarius region of the southern Milky Way, according to the study of 
the newly discovered variable stars, are at approximately the same dis- 
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tance as the center of the Galaxy, previously deduced from globular 
clusters. ‘These rich star fields appear, therefore, to constitute a nucleus 
of the Galaxy, and around it the stars in the solar neighborhood move in 
a period of some hundreds of millions of years. 

The third note gave a new determination of the average absolute mag- 
nitude of long period variables, and indicated the probable usefulness of 
this common type of star in future measurements of galactic dimensions. 

2. The occurrence of dark nebulosities in the direction of the galactic 
center, shown in the composite photograph reproduced in figure 1, raises 
some doubt as to the security of photometric estimates of distances in 
that direction. It is probable that these obscuring nebulosities, roughly 
sketched in figure 2, are mainly within a thousand parsecs of the earth. 
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FIGURE 2 


A rough sketch to indicate the general distribution of obscur- 
ing nebulosity shown in figure 1. The patch around 324°,+15°, 
is the Rho Ophiuchi nebulosity. 


If there are similar nebulous regions intermingled with the star clouds or 
lying beyond the nucleus, they may make completely hopeless the attempt 
to measure the bounds of the Milky Way beyond the center. 

Fortunately we can test the transparency of the Galaxy in any specific 
direction with the aid of the extra-galactic nebulae; at the same time we 
can make a definite contribution to the question of the distribution of 
these outside galaxies. 

3. Long-exposure photographs with the Bruce 24-inch reflector are 
available for nearly half of the central region. Each plate covers about 
thirty-five square degrees, and all are capable of showing nebulous objects 
for which the integrated photographic magnitude is as bright as 175; 
for most of the plates the limit is approximately 18"0. The centers of 
these long-exposure plates are indicated by crosses in figure 3, which 
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shows also the distribution of the extra-galactic nebulae that are now 
known for this region. Plates of shorter exposure (60” and 120”) cover- 
ing the whole of the central region have been examined, but most of them 
are limited to photographic magnitudes brighter than 16.5 and therefore 
are relatively impotent for the majority of extra-galactic nebulae. 

4. The systematic examination of plates for new nebulae in the region 
of the galactic center has been carried out mainly by Professor S. I. Bailey,' 
Dr. D. H. Menzel,’ and Miss Sylvia Mussells (results unpublished). The 
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Extra-galactic nebulae in the central Milky Way region. Large dots represent 
N. G. C. objects; small dots, those in the I. C. or in published Harvard lists; 
small circles, new nebulae discovered recently at Harvard. The same symbols 
crossed by diagonal bars represent nebulae in each of the groups which have 
been definitely identified as spirals. The object at 353°, —20°, is N. G. C. 
6822. Crosses show the centers of long-exposure plates. 


total number of extra-galactic nebulae, as shown in figure 3, is 331, in- 
cluding those which fall somewhat outside the boundaries of the central 
region. All have been re-examined by the writer. Although most of the 
objects in lowest latitude are doubtful, as indicated by colons, there can 
be no doubt of the existence of a number of extra-galactic nebulae within 
twenty degrees of the galactic center and within ten degrees of the central 
line of the Milky Way. The object at 353°, —20°, is N. G. C. 6822, a 
well-known extra-galactic system of the Magellanic Cloud type. 

Of the 331 objects plotted in figure 3, about twenty per cent are definitely 
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spirals; an additional twelve per cent are probably spirals; and eight per 
cent are described doubtfully as extra-galactic nebulae. The remaining 
sixty per cent are unspecified extra-galactic nebulae, largely of the elliptical 
type or spirals in plan; their detailed classification will appear in a forth- 
coming catalogue of new nebulae. ; 

The doubtful objects may be remote star clusters, or planetary and 
irregular nebulae; or perhaps some are but fortuitous aggregations of 
faint stars. Larger scale photographs will be necessary to distinguish 
their character. 

5. From the distribution of long-exposure plates and recorded nebulae, 
shown in figure 3, the incompleteness of the data is obvious. We may 
reasonably assume that when plates are available many more extra- 
galactic objects will be found in low latitude, at least in the interval between 
ten and twenty degrees of latitude on both sides of the galactic circle. The 
effective obscuration is apparently limited to easily recognized dark 
nebulosities. 

Dependable magnitudes are not yet available for the extra-galactic 
nebulae of this region, but a preliminary examination shows that the 
relation of angular diameter to apparent magnitude is not greatly different 
from that in Coma-Virgo. For the 2775 objects in the Coma-Virgo 
region (high galactic latitude) for which measures are available, we have 
recently shown that the diameter-magnitude relation indicates that space 
is effectively transparent up to a distance of something like a hundred 
million light years. ‘The same now appears to be true in low latitudes as 
well, for in the places where the outside galaxies are seen the apparent 
magnitudes are not seriously dimmed by either general or differential 
space absorption. 

In summary, we conclude that though readily-observed obscuring 
matter effectively conceals some of the galactic star clouds in the central 
region, and even hides the center itself, the Galaxy is completely trans- 
parent at the borders of these obscuring nebulae, permitting the observa- 
tion of not only the most distant stars of our system but also the outside 
galaxies, probably millions of light years beyond the farthest edge of our 
Galaxy. 

1 Harv. Ann., 72, No. 2, 1913. 


2 Harv. Ann., 85, No. 6, 1924. 
3 Shapley and Miss Ames, Harv. Obs. Bull. 864, 1929. 
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THE VISUAL ACUITY OF THE BEE AND ITS RELATION 
TO ILLUMINATION 


By SELIG HEcHT AND Ernst WoLF* 
LABORATORY OF Bropuysics, COLUMBIA UNIVERSITY, NEw YorRK 


Communicated February 6, 1929 


1. Two reasons prompted us to investigate the relation between visual 
acuity and illumination in the bee. We wished to know, first, whether 
the variation of visual acuity with illumination as it exists in the human 
eye! is to be found in so differently constructed an organ as the insect 
eye; and second, whether the explanation suggested by one of us? for 
this relation in the human eye is of sufficiently general significance so 
as to be adequate for the bee’s eye as well. 

In the human eye visual acuity is poor at low illuminations; as the 
intensity increases, visual acuity increases with it at first rapidly and 
then slowly; and finally at high illuminations further increases in intensity 
produce no change in acuity. Visual acuity, since it represents the 
resolving power of the retina, depends on the distance which separates 
the receiving elements in the retina. In order to make this distance 
vary with illumination it has been assumed that the minimum illumina- 
tion necessary to stimulate the individual cones and rods—their threshold, 
in short—is distributed in the usual probability or statistical manner of 
populations. When expressed quantitatively these ideas describe in 
detail the data of human visual acuity. 

It is proposed to see whether these facts and ideas have any general 
validity in the physiology of the visual process. In order to do so we had 
to develop a method for the investigation of the vision of animals other 
than man. Starting with the common observation that animals with 
eyes respond to a sudden movement in their visual field, we converted 
it, in terms of the following considerations, into a method of measuring 
visual acuity. If the visual field of a sensitive animal is made up of a 
pattern of dark and illuminated bars of equal size, the animal will respond 
to a displacement of this field only when it can distinguish the components 
of the pattern. In case the animal cannot resolve the black and white 
bars, the field will appear uniformly illuminated and displacement of the 
pattern will elicit no response. If visual acuity varies with illumination, 
then the capacity to respond to these movements in the visual field will 
depend on the illumination and on the size of the pattern. 

The bee is sensitive to changes in its visual field, and responds by a 
reflex, sidewise movement of the head and thorax. If the bee is crawling 
on an inclined, transparent surface, below which is the luminous visual 
pattern, the response to a movement of this pattern becomes evident by 
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a sudden change in the direction of its progression, which is opposite 
in sign to the movement of the pattern. We prepared a series of plates 
composed of equally wide opaque and translucent bars, each plate having 
a different size of bar. The experiments then consisted in determining 
for each size of pattern the minimum illumination at which a bee will 
just respond to a movement of that pattern. The reciprocal of the 
visual angle subtended by each size of bar is then the visual acuity of the 
eye at the corresponding illumination. 
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FIGURE 1 
Relation between visual acuity and illumination. An effort has been made to repre- 
sent with a dot each measurement made with each bee. Since this is obviously 
impossible when the measurements come close together, the plot shows the general 
distribution of the data. 


2. We made measurements with 91 normal, worker bees. The data 
are given in figure 1, where visual acuity is plotted against the logarithm 
of the illumination. Each dot represents, as well as is graphically possible, 
a single determination with a single bee. It is at once apparent that 
visual acuity in the bee’s eye varies with the illumination in much the 
same way as in the human eye. 

Certain comparisons may be made between the two. The disposition 
of the data in figure 1 indicates that the experiments cover the whole 
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range of visual acuity of which the bee is capable. The maximum lies 
between 0.016 and 0.017. This is below the visual acuity of the human 
eye at the lowest perceptible illuminations. Our maximum visual acuity 
is about 1.5; our minimum is about 0.03. Incredible as it may seem, 
the bee’s greatest capacity for the optical resolution of its environment is 
never better than ours is at our worst. 

The disparity in the visual acuity of our eyes and those of the bee is 
even greater than this and is brought out by considering the visual acuities 
at the same illuminations. The maximum visual acuity in both cases 
occurs at very nearly the same intensities of illuminations and corresponds 
to a brightness of between 50 and 100 millilamberts. Our maximum here 
is about 1.5; the bee’s maximum is 0.017. We can therefore resolve the 
environment about 100 times better than a bee can. That this low value 
is no laboratory product is borne out by the experiments of Baumgartner® 
in which bees on the wing in the field were shown to have a similar, almost 
negligible, form discrimination. 

3. Differences in resolving power mean differences in the distances 
which separate the centers of the receiving elements. The data of figure 1 
would then require that the ommatidia be separated by a variable dis- 
tance which depends on the intensity. Since this cannot be true struc- 
turally the results must be interpreted in such a way as to secure a func- 
tionally variable separation of ommatidia which are structurally fixed. 

Let it therefore be supposed that the receptor elements in the ocular 
mosaic do not all possess the same threshold, but that the threshold varies 
among the ommatidia as does any other characteristic in a population. 
At low illuminations then, only a few ommatidia are functional. Since 
these are distributed at random, they will be far apart and will give the 
same result functionally as if there were no receiving structures between 
them. As the illumination increases more and more ommatidia become 
functional; the distance between functional elements becomes smaller, 
and the resolving power becomes greater. This continues until an il- 
lumination is reached when all the elements are functional, and no further 
increase in visual acuity can take place. Such an explanation obviously 
describes the data. But before it can be formulated quantitatively it is 
necessary to examine the structure of the eye in some detail. 

4. In these experiments the relation of the creeping bee to the visual 
field is such that the pattern is registered across the long axis of the eye. 
At an illumination when all the elements are functional, the maximum 
visual acuity will then occur when a horizontal row of elements receives 
light, and an adjacent row receives no light, and so on. The size of the 
smallest perceptible pattern will correspond to the visual angle which 
separates the centers of two adjacent elements. 

If the elements were all the same angular size the maximum visual 
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acuity could be registered all over the eye. But since the elements in 
the bee’s eye are not uniform in angular dimension, the maximum visual 
acuity can be obtained only at that position on the eye where the angular 
separation is a minimum. In our measurements the visual acuity cor- 
responding to any illumination is always the maximum visual acuity at 
that intensity. Therefore, at any illumination, no matter how many 
functional elements it represents, visual acuity is determined at that 
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FIGURE 2 ‘ 


Relation between visual acuity and illumination for bees’ with the 
central part of the eye painted out as shown in the figure. The points 
are individual measurements. The broken curve is the normal rela- 
tion taken from figure 1. The full curve is’ made from the normal 
curve by multiplying its ordinates by 0.75. 


. area on the eye where the angular separation between functional elements , 
is at a minimum. 

This concept is so important for understanding the vision of the bee 
that we tested it in several ways. The measurements as plotted in figure 1 
show that the maximum visual acuity of which the bee is capable at 
the highest illuminations is about 0.017. This corresponds to a visual 
angle of between 0.9° and 1.0°. Since at these illuminations all the 
ommatidia are functional, this experimentally determined, minimal, 





182 BIOPHYSICS: HECHT AND WOLF Proc. N. A. S. 


angular separation should correspond to the smallest vertical separation 
between adjacent ommatidia as determined anatomically. Baumgdartner* 
has recently measured the angular separation of adjacent ommatidia in 
the bee’s eye. In vertical section the ‘smallest separation is near the 
middle, in the lower half of the eye; it includes about 20 elements and 
its value lies between 0.9° and 1.0°. 

The evidence is even better if one examines the structure of the eye in 
detail. In vertical section, according to Baumgartner, the angular 
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FIGURE 3 
Relation between visual acuity and illumination for bees with the 
anterior half of the eyes painted out. The points are individual mea- 
surements. The broken curve is the normal curve of Fig. 1. The full 


curve is constructed from the normal by multiplying its ordinates by 
0.62. 


separation of adjacent ommatidia increases from about 1° at the center 
to about 4° at the periphery. The increase is gradual, the middle half 
of the eye constituting a region of small angular separation in comparison 
with the rest of the eye—a sort of fovea. If visual acuity is always medi- 
ated by the region of minimum angular separation, then the elimination 
of this central foveal area should depress the visual acuity function to 
the level of the remaining peripheral ommatidia. We measured the 
relation between visual acuity and illumination in 21 bees with a spot 
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of opaque black paint placed in the center of each eye. Such a spot of 
paint covers about one quarter of the area of the eye, and renders non- 
functional all ommatidia whose angular separation is less than about 
1.3°. Such bees should give a visual acuity of about three-fourths of 
normal. The data secured with these animals are given in figure 2. 
The broken line is the normal visual acuity; the continuous line is drawn 
so that its ordinates are 0.75 of the normal. It obviously describes the 
data, and supports the idea that visual acuity determination is a regional 
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FIGURE 4 


Comparison between visual acuity and number of ommatidia functional in a given 
angular distance in their relation to the logarithm of the illumination. ‘The relation 
between the number functional and log I is in the nature of an integral probability or 
distribution curve. 


function, and depends on the utilization of the part of the bee’s eye which 
structurally permits the maximal resolution. 

This resolution occurs on the long axis of the eye, because the pattern 
is received as a series of horizontal dark and luminous bars across the 
eye. For the eye to perceive the pattern as a series of bars, each bar 
must stimulate at least two functional ommatidia; and since the eye is 
very nearly symmetrical these ommatidia are most likely distributed 
on either side of the eye. If one side were rendered non-functional this 
would reduce the number of functional ommatidia to half, and in order 
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to have the same number of elements determine a bar as before, the 
width of the bar would have to be increased about twice, and the visual 
acuity function would be depressed to about half. We made measure- 
ments with 19 bees in which the anterior half of each eye was painted 
out. The results are given in figure 3, where each measurement is given 
as before. Through the points there is drawn a continuous curve whose 
ordinates are 0.62 of the values for the normal, unpainted eye. 

The longitudinal use of the bee’s eye is not fortuitous. As shown in 
figures 2 and 3, the eye is about four times as long as it is wide. Further- 
more the angular separation between adjacent ommatidia is more than 
three times as great in the horizontal meridian as in the vertical meridian. 
Both these facts would tend to make the bee’s eye an organ which functions 
essentially as a linear receptor. This we found to be true experimentally 
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FIGURE 5 


Distribution of thresholds of the ommatidia in the bee’s eye. The 
curve is the first differential of the number curve in figure 4, and is a 
differential probability or distribution curve. 


because we were unable to get any responses to a pattern arranged to 
register on the eye as bars parallel to its long axis. This confirms Baum- 
gartner’s findings that bees in the field are very astigmatic, and resolve 
their environment vertically with much greater accuracy than horizon- 
tally. 

5. We have interpreted the shape of the curve in figure 1 as meaning 
that the number of elements functional varies with the illumination. If 
the ocular mosaic were uniform, it would follow that since visual acuity 
is determined by the vertical distance between elements in the region of 
maximum ommatidial density, the curve in figure 1 represents the number 
of functional elements in the vertical axis of the fovea corresponding to 
any illumination. But the angular separation between adjacent omma- 
tidia is not constant. The precise way in which it varies must therefore 
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be considered in the conversion of visual acuity data into the number 
of ommatidia functional in the eye. 

Fortunately this is possible because of Baumgiartner’s anatomical 
study. By remembering that visual acuity is always determined in the 
region of the eye where the angular separation of ommatidia is a minimum, 
one can determine the actual number of ommatidia which are included 
in vertical section in any given visual angle. There are only about 70 
ommatidia in vertical section in the lower half of the eye, and one can 
lay off the angle occupied by each ommatidium on a linear scale and com- 
pute the actual number of ommatidia in vertical projection which must 
be functional in order to give a definite visual acuity. The results of this 
computation are to be found in figure 4, where the number of elements 
in vertical section required to produce a given visual acuity are given. 
The ordinates of the number curve have been arbitrarily multiplied by 
3.5 to make the two curves comparable. 

The number curve in figure 4 resembles the usual integral distribution 
curves of the statisticians, even as its first differential, the threshold curve, 
in figure 5 resembles the more commonly encountered differential dis- 
tribution curves. Therefore, we may make our hypothesis of the relation 
of visual acuity and illumination quantitatively specific by stating it as 
follows. ‘Taking the structural relations of the ocular mosaic as given 
by Baumgartner, our data relating visual acuity and illumination may be 
described with complete fidelity by assuming a distribution of the thresh- 
olds of the various ommatidia corresponding to the population curve of 
figure 5. 

This distribution curve may be interpreted in two ways. Assuming 
that a given threshold is a permanent characteristic of a given element, 
the curve in figure 5 then represents the distribution of this characteristic 
in the population of ocular elements. However, one may conceive this 
situation purely in terms of probability. Assume that the threshold of a 
given element is not fixed, but can vary over the whole range included in 
figure 5. Then the probable number of elements which have the same 
threshold at the same time is given by the curve in figure 5. This may 
therefore be used as the basis for visual acuity in precisely the same way 
as before. 

The full details of these experiments are to appear in the Journal of 
General Physiology. 

* FELLOW, INTERNATIONAL EpucaTION Boarp. 

1 Kéenig, A., Sitzungsber. k. Akad. Wissensch., 559, 1897. 

2 Hecht, S., Proc. Nat. Acad. Sci., xiii, 569, 1927; Jour. Gen. Physiol., xi, 255, 1928. 
3’ Baumgartner, H., Zeit. vergl. Physiol., vii, 56, 1928. 
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SOME PROBLEMS OF BUSINESS FORECASTING 
By C. F. Roos 


DEPARTMENT OF MATHEMATICS, CORNELL UNIVERSITY 


Communicated January 11, 1929 


1. Introduction.—For four years the people of the United States have 
lived in a steadily increasing prosperity uninterrupted by severe industrial 
reactions. It is true that there have been business failures but these have 
been of the minor sort. Production has had its usual seasonal fluctua- 
tions, but the swing from the high point to the low point has been less 
emphatic than formerly. After recovery from these short depression 
periods orders from the trade have invariably lifted volume of production 
and distribution to new high levels. Bankers and business men are 
beginning to wonder if these conditions will continue and if it is possible 
to scientifically forecast business conditions. 

In answer to an ever-increasing demand various new and extraordinary 
methods of forecasting have been devised. Many bureaus of economic 
forecasting are endeavoring to render an honest service by plotting curves 
of existing economic data and extrapolating to forecast future conditions. 
If we fit a curve to existing economic data, we can certainly expect it to 
explain this data, but can we hope to use such an empirical curve to fore- 
cast new phenomena? ‘The obvious answer is that we may provided that 
we have taken all of the influencing factors into account. ‘Thus, if we knew 
that the demand for clothing depended only upon the price of clothing, 
we could gather sufficient data in the form of prices of clothing and total 
sales of clothing at these prices to enable us to draw a curve which would 
accurately describe the relation existing between demand and price of 
clothing. It happens, however, that the demand for an article depends 
upon the time, the rate of change of price, the rate of production of the 
article and other factors (psychological and otherwise) as well as upon 
the price.! Since it is obviously impractical to plot a curve in more than 
two dimensions, the graphical method is at a disadvantage at the beginning. 

It is true that we can eliminate all influencing factors except the time 
by means of our supply functions, cost of production functions,” etc., but 
this elimination gives the price, the demand, etc., as transcendental 
functions of the time and not the simple polynomial expressions used in 
statistical work. ‘Thus, if we start with the simple assumptions that the 
rate of demand y(t) is a linear function of the rate of supply (production) 
u(t); that the rate of demand depends linearly upon the price p(#) and the 
rate of change of price dp/dt, that is, on eliminating y(#) we have u(t) = 
ap(t) + b + hdp/dt where a, b and h are constants; that the cost of 
production is a quadratic function Au? + Bu + C + D(du/dt)? + H(du/dt) 
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+ F p(t) + Gdp/dt, where A, ..., H are constants, and, finally, that each 
producer attempts to maximize his net profit over an interval of time ¢, 
to # discounted to the time 4, we are led to the following different kinds 
of transcendal price and production curves, 


P = —Po/Lo + Kie**™ + Kee*™ + Kye*™ + Kee*™ 
? = —Po/Lo + Kie*™ + Koe*™ + e*™(K3 cos xt + Kasin xg) (1.1) 


P = —fo/Lo + e** (Ky cos xot + Ke sin xot) + e*™ (Ks cos xit + 
K4 sin xt). ‘ 


where here and following %;, x2, x3 and x4 are positive constants or zero.® 
Furthermore, a slight change in interest rate or over-production may 
eause the price curve to shift from any one of these curves to any other 
curve. There is obviously more to the problem of business forecasting 
than the plotting and extrapolation of statistical curves. 

Just as an astronomer is practically helpless without a telescope so the 
quantitative business analyst, if he is to have a workable understanding 
of his own work, is helpless without an economic theory. Although such 
a theory may not be perfect; although it may assume perfect control of 
conditions, i.e., that the psychological influences in our business system 
can be measured, and although it may need much overhauling, never- 
theless a theory is necessary. Any theory which is to be of service must 
be able to explain existing phenomena and predict new ones. In the case 
of economics it seems that the kind of theory which can best do this is a 
mathematical one. It is true that other economic theories are valuable, 
but the mathematical theories have marked advantages in the definition 
and analysis of concepts and in the discovery and presentation of a general 
view of the economic field. 

Unfortunately, dynamical or time changes which we know occur cannot 
be taken into account by the statical or motionless theories of Cournot, 
Auspitz, Lieben, Walras, Pareto and others. These men realized only too 
well this limitation of their theories, but none were able to overcome it. 
It should be apparent that there is not ever a remote hope of being able to 
forecast by means of a statical theory which does not take time changes 
into account. To bring the general theory of economic equilibrium into 
agreement with dynamical facts as we observe them it is necessary to pass 
from the statical hypothetical equilibrium of Walras and Pareto to a real- 
istic treatment of an actual time equilibrium. Such a theory I have 
recently given.! 

2. Operation of a Mine or a Machine.—Many problems of forecasting 
of a very fundamental importance which could not be touched by the 
older theories are now theoretically simple. For example, it is now possible 
to tell (theoretically) how to operate a mine or machine and when to scrap 
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it in order that a maximum profit result to the owner. Thus, if the cost 
of production is Q(u,u',p,p’,t) and the demand-supply equation 


u(t) =ap+b+hp’, (2.1) 


where primes denote derivatives with respect to time, our problem is to 
determine the rate of production u(t) and the price p(t) satisfying certain 
end conditions 


6,(u(t,) ,U(te), p(t), P(te) tte) - 0 (u = i. o 0 og M < 6) (2.2) 


imposed on u(t) and p(?) at the initial and final times ¢, and f, so that the 
total discounted net profit 


te 
r= (pu — Q(u,u',p,p’ ,t)E(t,t)dt, 


t 
where E(t,,t) is the discount factor exp(— a 5(r)dr and 4(r) is the force 
hh 


of interest, is a maximum. If we fix 4 and f, and the price and rate of 
production at these times, our problem is a Lagrange problem with fixed 
end-points.‘ In many instances we will want to leave some, or all, of 
these end values variable. Thus, if we wished to determine the best 
time to discontinue our business or scrap our machine, as the case might 
be, f2 would be variable. The most general situation is typified by the 
conditions (2.2). 

The price and production curves which give a maximum profit are in 
any case the so-called extremal curves determined from the Euler-Lagrange 
equations of the calculus of variations. If Q(u,u’,p,p’,t) takes the form 
suggested in §1, these extremals belong to one of the four-parameter 
families (1.1), the particular family depending upon the relative sizes of 
the quantities A,..., H,a,b,handé.3 The four parameters Kj, .. ., Kg 
and the times ¢, and / will in the general case be determined by the m 
end-conditions (2.2) and 6 — m transversality conditions of the calculus 
of variations. 

If we can control the quantities A, ..., H, a, b, h and 6, we can accurately 
forecast the best rates of production and prices, i.e., the rates of produc- 
tion and prices which yield a maximum profit. In actual economic prob- 
lems these quantities will remain constant for certain short periods of time 
and then change abruptly. For such cases the profit integral would be 
replaced by an expression which would be the sum of several integrals, 
each of which represented a profit over an interval of time for which our 
coefficients A, . . ., 6 were constant. The real difficulty in forecasting 
lies in the fact that we have no information to tell us when these changes 
will occur and how much they will be. If we can obtain the step-functions 
which represent these coefficients, the whole problem of forecasting will 
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reduce to the problem of maximizing an expression which is the sum of 
several definite integrals, but at present it seems that only a gypsy fortune 
teller or other clairvoyant can tell us when A, . . ., 5 will change and 
how much. 

The present state of affairs seems to be something of the following sort. 
We can expect to forecast prices and rates of production for intervals of 
time of duration short enough that all of our coefficients remain constant. 
We will miss our forecast at the critical points when these quantities change, 
unless we miraculously divine when these changes will occur and how 
much they will be. As soon as we are able to determine our new state of 
affairs after a change, we can again make forecasts which will be accurate 
until another change in the coefficients occurs. 

3. The Problem of Replacement.—In a more general situation an 
operator will desire to replace a machine which is in operation by a se¢ond 
machine whose operating expense will in general be different. If 
Q:(u,u',p,p’,t) is the operating expense of the first machine and 
O2(u,u’,p,p’,t) the operating expense of the replacement machine, then the 
total discounted net profit over a period of time including the life of 
both machines, will be 


¥ ts 
T = fpu-Qulun'p.p" DEC +f (ou-Qaleun'p.2° DEC 


where ¢, is the time at which the first machine begins operations, 4% the 
time at which it is replaced and /; the time at which the second machine 
stops operating. If we fix ¢ this problem is, of course, simply the general 
problem of the preceding section. We may, however, wish to know the 
““‘best’”’ time at which to make the replacement. In this case ¢, would be 
variable and would have to be obtained from a corner condition of the 
calculus of variations. I have obtained these corner conditions under the 
assumptions that the second machine begins at the same time and same 
rate of production at which the first machine stops.’ Henry Pixley has 
shown that the problem is still possible when these continuity conditions 
are removed and has discussed such interesting problems as the problem 
of replacement when a lapse of time is allowed at the time of replace- 
ment to allow for installation, etc.® 

4. Maximum Profit and Minimum Loss.—Conceivably the profit 
integral s might not have a maximum. ‘To assure ourselves on this point 
we need only examine Clebsch’s generalization of the Legendre condition.’ 
This condition requires that the quadratic form T°0°F/Ou’du'’ + 2TW- 
OF /Ou'dp’ + W0°F/dp’dp’ be < 0 for all real values of T and W and 
all u, p and ¢ in a region defined by the end conditions (2.2) and the 
demand-supply equation g(u,u’,p,p’,t) = 0, and for F = (pu — Q(u, ...,t) 
E(t,t) + A(é) o(u,u’,p,p’,t). For a minimum the inequality must be re- 
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versed. If we assume the cost function suggested in the first section, 
the Clebsch quadratic form becomes for a fixed 4, 


—(T?D + W°E) 


and since D and E are positive this quadratic form is definitely negative 
for all values of u, p and ¢. D is the “‘idling’’ depreciation and hence is 
positive in general, and E is the cost due to changing prices and hence is 
positive.® 

We have just shown that the rate of production u(t) and the price p(é) 
determined from the Euler-Lagrange equations can never minimize + 
and satisfy the demand equation assumed in this paper. This, however, 
does not imply that z will always be positive, that is, that there will always 
be a profit. If the maximum of 7 is negative no possible choice of p(t) 
and’ u(t) will produce a profit. In other words the business is doomed 
from the beginning. If such a business must operate because of a fran- 
chise or otherwise, the best it can hope to do is to choose u(#) and p(t) so 
that they minimize its losses. 

5. The Problems of Codperation and Competition.—The discussion so 
far has been confined to the problems of monopoly. The analysis for the 
problems of competition is only slightly more general than that given 
here and hardly more difficult, and leads to the same general kinds of price 
curves.* The problems of monopoly, however, are about the only ones 
which can be expected to be of practical value. In any actual situation 
it would be practically impossible to get all manufacturers of a like com- 
modity to give up what they consider to be business secrets, and, if they 
did pool their records, many of them would be dissatisfied with the com- 
putations. No company would be anxious to believe that they could 
make the most profit by discontinuing production immediately, yet the 
computations might rightly show this to be the case. 

On the other hand any company which is powerful enough to safely 
neglect its competitors can obtain useful information. ‘Thus it is possible 
to predict how twenty or more (or less of course) plants operating for a 
large holding company should operate over a period of time for which 
the quantities A;,..., H;, (¢ = 1, ..., 20) remain constant so that a 
maximum profit result for the company. If each of these plants operated 
at a different cost of production, the computation would show the futility 
of continuing certain plants and would advise the enlargement of others. 
This would be indicated by the fact that for some plants u(#) would be 
very much below maximum capacity (i.e., these plants should be scrapped) 
and for other plants u(t) would be above the maximum capacity of the 
plant. . 

The small producer, who cannot produce the major volume of trade in 
his commodity can also forecast his rate of production and price, provided 
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he is willing to make a still further approximation. Thus, he must be 
willing to assume that the demand for his article depends only upon his 
own price and rate of change of price. This producer will then be able 
to set a partial equation of demand which will be valid so long as his 
trading is confined to certain limits. It will then be possible for him to 
determine at what rate he should operate and what his price should be 
in order that he make as large a profit as possible over a given interval 
of time. The various other problems already suggested will also be 
possible. 

1C. F. Roos, “A Dynamical Theory of Economics,” J. Pol. Eco., Oct., 1927. 

2 C. F. Roos, loc. cit., part IV, or H. L. Moore, ““A Theory of Economic Oscillations,” 
Quar. J. Econ., Nov., 1926. 

3C. F. Roos, “A Mathematical Theory of Business Fluctuations,”” Am. J. Math., 
not yet published. . 

4 J. Hadamard, Lecons sur les calcul des variations, 1910. ia 

5C. F. Roos, “A General Problem of Minimizing an Integral with Discontinuous ' 
Integrand,” Trans. Am. Math. Soc., 31, 58-70 (1929). a 

6 Henry Pixley, “Discontinuous Solutions of the Problem of Replacement,’’ Am. 
J. Math., not yet published. 
7 J. Hadamard, loc. cit., p. 391. 
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By measuring the rate of change of electromotive force of a cell of the 
type 






He Q, QO 
Hg, HgCl Sat. KCl Buffer solution | Au 
O, (1 atm.) 
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in which the hydroquinone (HQ) and quinone (Q) is kept saturated with 
molecular oxygen and rapidly stirred, one may follow the kinetics of the 
reaction 









O:2 + He Q —> H,0: + Q (1) 











more accurately and in more acid solutions than is possible by the method 
of oxygen absorption.! The complicating difficulties of side reactions, so 
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prominent in alkaline solutions are also eliminated. Denoting by X, the 
percentage of hydroquinone oxidized we have the relation 


E = K — 0.0296 log ae (2) 


connecting the momentary electrometric force E and the concentrations 
of H.Q and Q. K is a constant which can be evaluated for any given 
buffer by measuring the e.m.f. using quinhydrone (X = 50%). The 
initial values of X were always 2%, since in this range of concentration 
the potential was reproducible and the rate of change with time was of a 
suitable velocity. 





















































A- 0-0002 MMntt pH « 6-0 

B- 0-00 Mat? pH = 5-3 

C- 00005 Mmnt* pH = 5-7 
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When the values of the velocity of reaction dX /dt are plotted against 


the value of log me , which is proportional to the free energy content 


of the system, we find a linear relationship, which may be represented by 
the equation 


dX/dt = A + Blog me (3) 


where A and B are constants. Some of the more representative data are 
plotted in figure 1. In order to avoid side reactions, it is necessary to 
work in a pH range 5.3 to 6.3 and to introduce minute quantities of MnCl, 
as catalyst to bring the reaction to a suitable velocity. 
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Table 1 gives a summary of the results in terms of equation (3). 


TABLE 1 
[MNcL2] x E 
Pu 103 A B (EQUIL.) (EQUIL) 
6.3 0.05 —0.0250 0.0644 29.0 0.0696 
6.3 0.2 0.0777 0.2409 32.2 0.0716 
6.3 0.5 0.0840 0.5171 40.8 0.0763 
6.0 0.05 —0.0366 0.0505 15.9 0.0774 
6.0 0.1 0.0272 0.0686 28 .6 0.0871 
6.0 0.2 0.0770 0.1546 24.1 0.0841 
6.0 0.35 0.0614 0.2143 34.1 0.0903 
6.0 0.5 0.0731 0.2857 35.7 0.0912 
6.0 0.7 0.0444 0.3605 43 .0 0.0952 
6.0 1.0 0.0182 0.4545 47.7 0.0976 
5.7 0.2 —0.0381 0.0842 26.1 0.1033 
5.7 0.5 0.0329 0.1586 38.3 0.1106 
5.7 1.0 0.0264 0.2886 44.8 0.1140 
5.3 0.2 —0.0079 0.0308 35.7 0.1320 
5.3 0.5 0.0264 0.0748 30.7 0.1292 
5.3 1.0 0.0071 0.1214 46.7 0.1379 


The extrapolated values for zero velocity are computed from the rela- 
tions: 
ated) «2. eet) = Ke (- 4) 
+10 38 B 
The rate of reaction in this pH range is directly proportional to the [Mn*+*] 
and practically inversely proportional to the [H*]. 

Although van’t Hoff expected that reaction velocity should be related 
to the driving force, abundant experience has shown that the velocity of 
homogeneous chemical reactions is proportional to some integral power of 
the concentrations of the reacting substances. As far as we are aware, 
this is the first case in which the velocity of a homogeneous chemical reac- 
tion has proved to be proportional to the logarithms of the concentrations 
of the initial and final states or, what is the same thing, the energy levels 
of the system. 

To account for the failure of van’t Hoff’s prediction, it has been cus- 
tomary to assume the formation of intermediary complexes (collision 
complexes), the stability of which control the rate of the reaction. It is 
suggested that in the present case no such rate-controlling complexes are 
formed and that reaction proceeds by the following steps: 


H.Q —> Q + 2H+ + 25, 


which is strictly reversible followed by the oxygen molecule accepting the 


electrons, viz., 
2E + O. —> OO" 


2H+ + O.m — H.20x. 
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Thus the total reaction is the result of two forces: the tendency of one 
system to be oxidized, and the tendency of the other to be reduced. Since 
the environment is held constant in any one given experiment, the force 
controlling the velocity of the reaction is the electron pressure E. 

1 La Mer and Rideal, J. Am. Chem. Soc., 46, 223(1924). 


[ConTRIBUTION No. 597 FROM THE DEPARTMENT OF CHEMISTRY, COLUMBIA UNI- 
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THE TEMPERATURE COEFFICIENTS OF SOME 
PHOTOCHEMICAL REACTIONS 


By GEorGE B. KISTIAKOWSKY 
LABORATORY OF PHYSICAL CHEMISTRY, PRINCETON UNIVERSITY, PRINCETON, N. J. 


Communicated February 14, 1929 


A theory of gaseous explosions, first suggested by Christiansen and 
Kramers! and recently developed by Semenoff? assumes that in many 
cases the formation and propagation of an explosion-like reaction in gases 
may be due to the thermal reaction chains becoming infinitely long. By 
infinitely long is meant here that, as the ultimate result of the reactions 
of one primarily activated molecule (or atom), on the average one or more 
than one new activated molecule (or atom) is formed, the resulting chain, 
instead of being limited in length, spreading out to infinity, theoretically 
at least. Experiments of Sagulin,* as interpreted by Semenoff,? indicated 
in particular that explosions of oxygen-hydrogen and of oxygen-carbon 
monoxide mixtures belong to this type of gaseous reactions. Marshall,‘ 
on the other hand, has shown that the quantum yield in the photochemical 
water or hydrogen peroxide formation, sensitized to the wave-length 
2536 A by excited mercury atoms, does not exceed a few molecules. 
The chains in this reaction at room temperature are therefore very short. 
Furthermore, Andrejeff> and Coehn and Grote’ have found that the 
temperature coefficient of the non-sensitized photochemical formation of 
water is only 1.04-1.06 in the range 130-240°. The temperature thus 
seems to have only a very slight accelerating effect on the rate of this 
reaction. By analogy with the data of Marshall one can expect that the 
quantum yield of the non-sensitized reaction is also of the order of unity. 
Thus, this photochemical reaction would appear to have only very short 
chains even at the explosion temperature, whereas according to Semenoff 
the thermal chains become at this temperature infinite. Such divergent 
behavior, although quite possible, is not very satisfactory from the theo- 
retical point of view and a study of the non-sensitized photochemical 
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water formation at varying temperatures has, therefore, been made. 
The rate of carbon dioxide formation from dry oxygen and carbon mon- 
oxide has been similarly studied. 

Since the experimental details will be described elsewhere, it will be 
sufficient to mention now that the rate—in the initial stages of the re- 
action only—has been followed manometrically, a thin-walled quartz 
vessel being placed in a suitable furnace with an opening for illumination. 
As the source of light, a high-powered condensed spark burning between 
aluminum electrodes and cooled by a strong blast of air was used. ‘The 
radiation from the spark on the way to the reaction vessel passed a layer 
of air of about 1 cm. thickness. Of the strong lines in the aluminum- 
spark spectrum the only lines absorbed by oxygen which come into con- 
sideration are a triplet between the wave-lengths 1719 to 1725 A and a ig 
group of lines in the range 1854 to 1862 A. The intensity of the first i 
group of lines, which are strongly absorbed by oxygen,’ is very consider- a 
ably decreased, before it reaches the reacting mixture by absorption in 
the layer of air, and in the quartz wall. The second group of lines is 
absorbed only to a slight extent by the oxygen in the reaction vessel,’ li 
the result being that it is impossible to decide at present whether oxygen ¥ 
atoms® or activated oxygen molecules are mainly formed on light absorp- 4 
tion in the experiments. — a 

Two circumstances affected unfavorably the precision of the measur- . 
ments: inconstancy of the light source and a gradually increasing opaque- 4 
ness of the quartz reaction vessel to the active radiation, the latter effect 4 
being particularly pronounced when the quartz vessel was being heated. a 
Accordingly, averages of three to four runs were taken as the reaction a 
velocity at a given temperature and the velocity was followed, going up 
and down on the temperature scale. The following tables 1 and 2 repre- 




























TABLE 2 
CaRBON MONOXIDE-OxYGEN REACTION 












TABLE 1 
HypROGEN-OxYGEN REACTION 















TEMPERATURE TEMPERATURE TEMPERATURE TEMPERATURE 
INTERVAL COEFFICIENT INTERVAL COEFFICIENT a 
i 
25-100° 1.055 25-100° 1.005 ? i 
i 









100-200° 1.045 100-200° 1.026 
200-300 ° 1.07 100-300° 1.047 

300-400° 1.07 200-300 ° 1.045 | 
400-490 ° 1.13 300-410° 1.07 
490-504° 1.25 410-480° 1.06 
290-440° 1.11 480-504 ° 1.10 
440-480° 1.12 504-527° 1.23 







480-505 ° 1.20 








sent the results, the averaged temperature coefficients kr +10/kr being 
- given together with the temperatures between which they were measured. 
All experiments were carried out with stoichiometric mixtures of the gases 
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and in all experiments, except the last three of the first table, the pres- 
sure was near to 800 mm. in those three experiments the pressure was 
360 mm. The dark velocity, an average of the rate before and after 
each illumination, was subtracted from the total rate and it has been 
ascertained that during illumination the temperature of the reacting gas 
did not rise by more than 1° over that of the containing vessel. The 
dark rate was quite small in the case of the carbon monoxide—-oxygen 
mixture and did not exceed 5 to 8 percent of the total rate even at the 
highest temperature studied. In the case of hydrogen-oxygen mixtures 
on the other hand, the dark rate amounted sometimes to more than 25 
per cent of the total and was particularly high in the three last determina- 
tions of the first table for which a different quartz vessel was used. 

These results show that the temperature coefficients of the reactions 
studied are far from following the Arrhenius equation 


dlnk 
dT 





= Q/RT?, 


but, instead, increase with temperature and the more rapidly the higher 
the temperature. Accordingly, the “‘activation energy,” constant Q of 
the above equation will be found to increase from a few hundred calories 
at room temperature to 25,000 cal. at 500°, if a definite meaning can be 
attached to this constant in the present case. Although Wartenberg® 
has shown that the absorption of light by oxygen increases with tempera- 
ture, the effect here observed can hardly be due to this cause. The shorter 
wave-lengths of the aluminum spectrum are absorbed completely already 
at room temperature and the longer ones are absorbed at this tempera- 
ture in the 5 cm. deep reaction vessel, to at least 3 to4 percent.’ Thus, 
even if these lines were completely absorbed at some higher temperature, 
the maximal reaction velocity would be only a 20 to 30 fold of that at 
room temperature. An increase up to 40 fold was obtained without 
any indications of a limiting velocity being observed. Still, a minor 
fraction of the temperature coefficient may be due to a change in absorp- 
tion and it is intended to investigate this point in detail. The main in- 
crease in the reaction velocity is undoubtedly due to an increasing quantum 
yield. The latter has been determined now at room temperature by com- 
paring the rates of water and of carbon dioxide formation with that of 
ozone formation. Assuming that the volume change in the two former 
reactions takes place according to the equations: 2H: + O2: = 2H:0 
and 2CO + O2 = 2CO, and that the quantum yield in case of ozone is 
exactly equal to two, an identical yield of about 2 molecules of H,O or 
CO; per absorbed quantum has been obtained. It becomes thus evident 
that, at the temperatures in the neighborhood of 500°, quite long chains 
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are formed on illumination in both gaseous mixtures. The discrepancy 
between the thermal reaction theories and the earlier photochemical data 
is therefore eliminated to a considerable extent. The possibility of chains 
becoming infinitely long at the point of explosion, although made plausible, 
cannot yet be answered definitely. 

The present knowledge of the mechanisms of the photochemical water 
and carbon dioxide formation is insufficient for a detailed interpretation 
of the observed abnormal behavior of the temperature coefficients. A 
general interpretation may be indicated, however. The total rate of the 
photochemical reaction may be conveniently considered as being com- 
posed of two parts: first, the immediate processes following light ab- 
sorption and requiring for their completion either none or a very small 
additional thermal activation energy; second, the secondary chain mech- 
anism which may be started by some of the primary processes if a sufficient 
additional thermal activation energy is supplied.'° This secondary re- 
action may be quite identical in mechanism with the normal thermal 
reaction. The composite reaction of this type can be described by a 
general equation of the type: 


k = IpAe~/F7(1 + Be“ 87) 


from which it will be seen that if Q; < Qs, the temperature coefficient, 
instead of following the Arrhenius equation, will rise with increasing tem- 
perature. Some preliminary observations support the idea as to. the 
reaction mechanism here expressed by indicating that the kinetics of the 
photochemical formation of carbon dioxide change at higher temperatures. 
Thus, for instance, water vapor acts apparently as an accelerator at high 
temperatures, although at room temperature it has no effect on the photo- 
chemical rate.'' Experiments on the kinetics of the two reactions here 
discussed are now in progress. 


1 Christiansen and Kramers, Z. physik. Chem., 104, 451 (1923). 

2 Semenoff, Z. Physik, 48, 571 (1928). 

3 Sagulin, Z® physik. Chem., 1B, 275 (1928). 

4 Marshall, J. Phys. Chem., 30, 34 (1926); Ibid., 30, 1078 (1926). 

5 Andrejeff, J. Russ. Chem. Soc., 43, 13842 (1911). 

6 Coehn and Grote, Nernst-Testschrift, (1912). 

7 Eucken, Z. phys. Chem., 107, 436 (1923). 

8 Birge and Sponer, Phys. Review, 28, 259 (1926). 

® Wartenberg, Phys. Z., 11, 1168 (1910). 

10 A similar idea has been communicated some time ago to the writer by Professor 
Semenoff in another connection. 
11 Coehn and Tramm, Ber., 54, 1178 (1920). 
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LINEAR TRANSFORMATIONS IN HILBERT SPACE. 
I. GEOMETRICAL ASPECTS! 


By M. H. Stone 
DEPARTMENT OF MATHEMATICS, HARVARD UNIVERSITY 


Communicated January 17, 1929 


In this and succeeding notes, we summarize and extend some recent 
developments of Hilbert’s classical theory of transformations in space of 
infinitely many dimensions. Any thorough discussion of the current 
formulation of the quantum theory depends upon the theory to be de- 
scribed.? Notable results in the direction of a complete treatment have 
been established by J. von Neumann, who has dealt with transformations 
in real space.* Since the quantum theory involves transformations of 
complex space, a more extended discussion appears to be desirable. The 
methods of von Neumann seem to be incapable of application to the more 
general problem, and cannot in any case be considered as the most power- 
ful or satisfactory for the end in view. By a suitable combination of 
fundamental methods introduced by Carleman‘ and ideas elaborated by 
von Neumann, we are able to deal successfully with transformations in 
complex space. The advance is made possible by the introduction and 
utilization, ab initio, of the inverses of certain transformations. The 
simplification of analytical details brought about in this way can hardly 
be over-emphasized. 

A complex abstract Hilbert space § is a set of elements, denoted here- 
after by f, g, . . .. which is characterized by five properties: linearity, the 
existence of a unitary metric, the non-existence of a finite basis, complete- 
ness, and separability. The metric is determined by a numerically 
valued function, Q(f,g), defined for all pairs of elements of §; this function 
is an Hermitian symmetric bilinear function of its arguments and has the 
property that Q(f) = Q(f,f) is real and never negative, vanishing if and 
only if f is the null element; the square of the distance watenen the ele- 
ments f and g is O(f — g). 

A transformation or operator T is a correspondence between two not- 
empty subsets, § and D, of § such that to each element of § corre- 
sponds one element of D and that each element of D is the correspondent 
of at least one element of §. The sets § and D will be called the field and 
domain of T, respectively. If to each element of the field corresponds an 
element of the domain which is the correspondent of no other element of 
the field, the relation between the two sets is reciprocal and defines a 
transformation R which is called the inverse of T. 

We consider only linear transformations with fields everywhere dense 
in §; a transformation T is said to be linear if its field is a linear set and 
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if T(afi + def.) = aT f, + aeTfe for all complex numbers a; and a2 and 
all elements f; and f of its field. ‘The linear transformations which satisfy 
the condition of boundedness, Q(7f) = C?Q(f) with C a positive constant 
independent of f, are of considerable importance since they are handled 
with particular ease; a bounded linear transformation whose field is 
everywhere dense in § may be extended by continuity to the entire space. 
Associated with a given linear transformation T is a unique linear trans- 
formation 7*, the adjoint of 7; if there exist elements g and g* such 
that O(Tf,g) = Q(f,g*) for all f in the field of 7, then g is in the field of 
T* and T*g = g*, and conversely. If the field of T* is everywhere dense 
in §, the adjoint of T* is defined and is either identical with T or is an 
extension of JT. We shall consider mainly two particular types of linear 
transformation: the self-adjoint transformations T, for which T and T* 
are identical; and the unitary transformations U, for which U and U* 
are defined throughout § and are inverses of one another. As special cases 
we mention the following: the null-transformation or operator O, which 
takes every element of § into the null element; the identity J, which 
takes every element of § into itself; and the projection or special operator 
E which projects every element of § onto a linear subspace, so that E°f= 
Ff. 

The fundamental problem of the theory of linear transformations is to 
prove the existence of linear subspaces of § invariant under a given 
transformation T and to determine these subspaces. We consider here 
the first part of the solution of this problem for self-adjoint transformations. 

If T is a self-adjoint transformation, the transformation 7; = T — LJ, 
where / is a complex parameter, is easily shown to have a unique bounded 
linear inverse R; whose field is §, if / is not real. More generally, the 
following theorem can be asserted: 

THEOREM.—A set of necessary and sufficient conditions that a family 
of linear transformations X; defined over § for all not-real / coincide with 
the family of inverses, R;, of some family of transformations T; = T — lI, 
where T is a self-adjoint transformation, is the following: 

(1) X;, is the adjoint of X,,, if 1 and m are conjugate complex numbers; 

(2) (lL — m)X,:X,, = X; — X-, for all ] and m; 

(3) there is at least one value of / such that X;f = 0 implies that f = 0. 
The proof is geometrical in character. By similar methods, a more 
refined result can be established, to wit: 

THEorEM.—If T is a self-adjoint transformation, the points of the 
l-plane fall into three mutually exclusive sets A, B, C, with the following 
properties: 

(1) A is an open set containing all not-real points of the /-plane and, 
possibly, some real points; if / is in A, then T; has a unique bounded in- 
verse R, with field §, and, if J is also real, R, is self-adjoint ; 
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(2) JB, unless it is empty, contains only real points; if / is in B, T, has 
a unique inverse R; which is an unbounded self-adjoint transformation 
with a proper subset of § as its field; 

(3) C, unless it is empty, is a finite or denumerably infinite set of real 
points; if / is in C, then there is at least one element of the field of T for 
which T7;f = 0 while Q(f) > 0. 

There seems to be no simple geometrical argument which will yield a 
demonstration of the fact that the set B + C is not empty; even in n- 
dimensional space an appeal to the fundamental theorem of algebra is 
necessary to avoid difficult proofs. In the present discussion we will 
take up the question at a later stage. 


1 Presented to the American Mathematical Society, October 31, 1928. 

2H. Weyl, Zeit. Physik, 46 (1927-28), pp. 146, and Gruppentheorie und Quanten- 
mechanik, Leipzig, 1928; J. von Neumann, Géttinger Nachrichten, 1927, pp. 1-57 and 
245-272. 

3 J. von Neumann, /oci citati; and an unpublished paper which is to appear in the 
Mathematische Annalen, cf. Géttinger Nachrichten, 1927, pp. 1-55, footnotes 12 and 27. 

* T. Carleman, Equations intégrales singuliéres, Uppsala, 1923. 


ON COMMUTATION RULES IN THE ALGEBRA OF 
QUANTUM MECHANICS 


By NEAL H. McCoy 
DEPARTMENT OF MATHEMATICS, STATE UNIVERSITY OF IOWA 


Communicated February 7, 1929 


The algebra of the new quantum mechanics is a special form of the 
algebra of matrices. For one pair of canonically conjugate variables 
its properties are determined by the assumption 


Pq —- I=, (1) 


where g and # are matrices which represent the coérdinate and momentum, 
respectively. The value of the constant c plays no part in the develop- 
ment of the algebra although in the quantum mechanics it is assigned the 
value h/2z77. 

The theory of differentiation of functions of a matrix has been discussed 
by Heisenberg, Born and Jordan.' Let f be a function of the argument 
matrix x. The definition of derivative used in this paper is 


Fe) _ tim } Uf(x + a1) — f(x)]. 
dx oo >0a 


Partial derivatives are defined in a similar manner. Born proved by 
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induction that ——— <s ¢) =/2 +4 yg. Iff is a function of the two matrices, 
dx 
: ; te O° °F 
x and y, it may be shown in a similar manner that = ——. These 
OxOy dy0x 


results are rigorously established for polynomials only, although they are 
formally correct for other analytic functions. 

Let f be an arbitrary polynomial in p and g. Born? obtained the follow- 
ing formulas for commuting f with p or q 


i -fp=c2 
q 

of (2) 
aoa ae g 


If g denotes another polynomial it may be shown that these formulas 
are only special cases of the more general one 
poof dof vel 
-@gf@=2.-(-—-—- = * 3 

aes pa s! ae dp? dqt apr 8) 
We have established the validity of this formula for polynomials by 
induction. It is formally correct for other analytic functions and thus 
gives correct results if the series which arise are convergent. For the 
special case, f = p", 9 = q”™, this formula gives the following result which 
was obtained directly by Born and Jordan*® 


re -> e's! (3) () pr, 


the sum being extended to the smaller of » and m. Some interesting 
formal results may be obtained by means of formula (3) and Taylor’s 
series. The following are easily verified: 


f(p,ge™ = e™ f(p + nc,q) 
f(piq + che? ** = e? fp + ,q). 


The first of these was given by Dirac‘ and used in his theory of the hydrogen 
atom. 

Functions of more than one pair of canonically conjugate variables® 
arise in the quantum mechanics. Formula (3) becomes in the case of 
two pairs of variables 


fe — of = Xs - Hof) H® (f,¢)] 


where (5) 


! 3° a 
H(g, Ss ¢g 5 


(4) 
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; (s) . : : (22 of , Ww af 

It is seen that H”’ (g, f) is symbolically the expansion of | — — + —— }, 

Og Op, Oz Ope 
which suggests a method of extension to functions of three or more pairs 
of conjugate variables. 

Proofs and further applications of the formulas appearing in this paper 
will be published elsewhere. 

1 Born and Jordan, Zeits. Physik, 34, 1925 (858-888); Born, Heisenberg and Jordan, 
Ibid., 35, 1926 (557-615); Born, M., ‘Problems of Atomic Dynamics,’’ M. I. T. Lec- 
tures, 1926. 

2 Born, loc. cit., p. 79. 

3 Born and Jordan, loc. cit., p. 873. 

4 Dirac, P. A. M., Proc. Roy. Soc., 110, 1926, p. 566. 

5 Born, loc. cit., p. 94. 
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Let 
Qin 
f= ¢' 
and = 
eq ) 
e«7e -~£) 
where 7 is a primitive root of the prime / and set 
1-3 
—2n —4n “2 —x¢€ — 3), 
E 1+ s7 + s*r +...4+5 r 


=e€ > 


the symbol s representing the substitution ({/{’) in the notation of the 
Kronecker-Hilbert symbolic powers; n = 1,2,...(1 — 3)/2. 

If the field k(¢) is irregular and none of the above units is the /-th power 
of a unit in k(¢), then this field is said to be properly irregular. I shall 
prove in another paper’ that if k(f) is properly irregular then the second 
factor of the class number of this field is prime to /.* If w is an integer 
in k(¢) and q is an ideal prime then 


w V-* & 1 (mod q) 


where N(q) is the norm of q in k({) and (w) is prime to q. Also there is 
one and but one integer b in the set 0, 1, ..., 2 — 1, such that 
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N(q)-1 
w ! =? (mod q) 


provided q is prime to (/). We write 


fe 


In the present note I shall indicate how these power characters may be 
determined, if k({) is a properly irregular cyclotomic field and w is a singular 
number in k(f), and q is any ideal in k(f). An integer w in the field 
k(£) is said to be singular if and only if (w) is the /-th power of some ideal 
in the field not principal. 

First let E,, = w, which is included among the singular numbers. Now, 
if a is a rational integer,0 <a </ — 1, 


re —] d,~ ™* log ve) 
t-2n i — 2ny ° IT — 2n 
21i+a (a + 1) ) dv 


En — pind Ey). 
{Fy 


(x) => ge @ + Vn + ind (e” + D 
h 


ind E, = od 1; (1) 





where 


q is an ideal prime in k(¢) whose norm is q', g is a primitive root of q, 
h ranges over the integers 0, 1, 2,..., ¢ — 2, excepting (q — 1)/2 if q 
is odd and excepting zero if g is even. The primitive root g is selected 
so that 

gv-D/ = ¢ (mod q). 
If 

Sh + 1 = g; (mod q)* 


0<k<d — 1, then we write 


ind (g* + 1) =k. 
The symbol 
di.~ *" log W,(e’) 
dy! — 2n 


means that the (J — 2mn)th derivative of log W*(e’) is taken with respect 
to v and v = 0 substituted in the result, e being the Napierian base.* 
Let ¢(e’) and ¢,(e’) be two rational functions of e” with rational integral 
coefficients connected by the relation 





e 


@) = a0) + (G4) x 
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such that 
$(1) = ¢i(1) (mod I * *) 
and 
¢(1) # 0 (mod J), 
then 


¢(a) = oi(a), 


and we may prove that 


_ log o(e”) Pe ” cg log ¢:(e”) 


dymii at «mod aie : (2) 





We also have 


val) = X Aes (2a) 
where c ranges over the integers in the set 1,2,...J—1, such that 


cc’ = 1(mod J), 
where 
| — ac’|] +|—c’| >1. 


| x | is the least positive residue of x, modulo / and c’ is in the set, 1, 2, 
...,l— 1; also, the q, represents the ideal obtained from q by the sub- 
stitution (¢/{°). Assume that the q belongs to the exponent jl’ ~ 1 where 
j is prime to / so that qh *e- = (w), where w is an integer in k(f), 
such that w = 1 (mod (1 — a)). Raising both sides of (2a) to the power 
gj —*(l — 1) we obtain 
=i 
(Wa(S))* 
where 7 is a unit in k(f). It is easily proved that 7 is a power of ¢£ so 
that we have the identity 


(-) = 9X a 


, an 
(Hale D = + wl") XV S 
é ee 


’ 





where V is an expression of the type 
= bye”, 


where s ranges over a finite number of positive or negative integers or 
zero and the bs are rational integers. By (2) 


mi? v 
(I ae 1)jl' -1 dy log Vole ) 
dv™ 


mit v : 
=> dy log we(e ) (mod ‘hg + 1) 
dv™ 
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and 


i. a> di! tog ale) 
dyv™ 


R ml" v 
ay cm “eee? (mod /' +4), 
¢ du™ 


After a number of transformations the expression uo c™” reduces to a 
form which yields from (1) 


ee (t — 2n)it 
r 1 B® log we") (4, od 2) 


ind. E, = 
2j — nl yt - 2 


(3) 








where ¢ = (2n — 1)I'~?. 

Since the E,, » = 1,2,3,..., (1 — 3)/2, form a system of independent 
units in k({), then any unit in this field may be expressed as products of 
powers of these E’s, the exponents being fractions whose denominators 
are prime to / since k({) is a properly irregular cyclotomic field. Hence 
the power character of any unit in the field k(¢) with respect to any prime 
ideal q prime to / may be expressed by means of (3). 

Suppose now that @ is a singular integer in k(¢) which is not necessarily 
aunit. Since the field is properly irregular, the second factor of the class 
number is prime to / so that if q is a prime ideal in k({) prime to 6 and / 
then 

(qq-1)* = (&), 


where ¢ is an integer in k({), g is a rational integer prime to /, and q_, 
represents the ideal obtained from q by the substitution (¢/¢7~ 1). Hence 
by the law of reciprocity we have 


tad "Woe er ® 


I= (1 — @). 


lath 


we note that this may be written as 


oo. 


00_ 4 = on, 





To reduce the symbol 


But since 
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where 7 is a unit in k(¢) and \ an integer in that field, we have 


6, =m e-? 


= Meats id 

q [) tq 

and the last symbol on the right-hand side is determined by means of 
(3); hence {0/q} is completely determined. 


so that (4) gives 


The relation (3) enables us to determine the composition of the irregular 
class group in a properly irregular cyclotomic field, that is, we may deter- 
mine directly by the use of rational integers only, the invariants of this 
group as well as fundamental properties of the basis. We may prove the 

THEOREM.—Consider the Bernoulli numbers B,, Bo, . . .. By~3and let 

2 


the distinct Bernoulli numbers in this group whose numerators are divisible 
by 1 be 
B,,, Ba, . . +» Ba,. 


Further, let the Bernoulli number 
thi(2a; — 1) +1 
Bz My d = 2 


be divisible by "i but not by ++. Then it is possible to find a basis for 
the irregular class group in the properly irregular cyclotomic field k(¢) 





C, C2, eo C, 
cH a1 ete, 
and such that 
1 — 2a; 
iy gt = 1 
i= 1, 2, »¢ 


Corresponding to this basis there exists a basis for the singular primary 
numbers in k(t). It consists of the units 


1-3 
l1—1—2n i—1—4n — 27 §$-1-—-(—3)n 
td > r 


E =e@t®? +s oe tS 
n 


% = Qh, Aa, . . ») Be 


If j; is an ideal belonging to the class C; and 


i” = (w) 
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with 
wo=b thao the+h_2f 7? 
and 


we’) =bo the +...4+ b, 9°? - 


Then, ifk <1 —1, 


E ne w;(e”) 


dy" | = ." . ond ) 


fork #1 — 2a,, 


and 





d' ~ *% log w,(e”) 
dv! — 2a; 


|, 0 (moa 


oe eee 


We may deduce from this theorem the following 
CorROLLARY. If (w) ts the lth power of an ideal in the properly irregular 
field k() then w is a singular number if and only if 


lL — 2a; (0? 
[? 7 | 0 ae ie Y 


re a) Np eee 





for at least one value of 1. 
Many known results concerning cyclotomic fields are included in the 

above results. For example, certain theorems of Takagi.‘ Each of the 

members of the bases of singular primary numbers in the field, the specific 

units given above, defines an absolute class field of k(¢). ‘Through some 

computations carried out recently, Miss E. T. Stafford and the writer 

have determined a number of properly irregular cyclotomic fields, namely, 

for 1 = 37, 59, 67, 101, 103, 131, 149 and 157. The work was carried 

through for all primes / less than 200. No irregular fields other than 

properly irregular fields were discovered. Full details are given in several 

papers to be presented for publication in the Transactions of the American | 

Mathematical Society. / 
1 Bull. Amer. Math. Soc., 1929. 
2 Fiieter, Synihetische Zahlentheorie, 2nd edition, p. 223; Bull: Nat. Res. Coun., 

No. 62, February, 1928, Report of the Committee on Algebraic Numbers, II, pp. 34, 

44-45. 

3? Kummer, J. Math., 44, 1852, 121-130; Vandiver, Ann. Math., 1929. 

4 Takagi, J. reine angew. Math., 157, 236 (1927). 
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THE THERMODYNAMICS OF GASES WHICH SHOW 
DEGENERACY (ENTARTUNG) 


By Grpert N. Lewis AND JOSEPH E. MAYER 
CHEMICAL LABORATORY, UNIVERSITY OF CALIFORNIA 


Communicated January 11, 1929 


In our three preceding papers, ' we have studied the fundamental relations 
of thermodynamics and statistics. In order to ascertain the thermo- 
dynamic properties of particular systems, we employed an assumption 
of extreme simplicity and generality which, so far as is now known, is 
consistent with present experimental data. Therefore while that assump- 
tion may still be regarded as not actually proved, it is desirable to obtain 
all of its important consequences. In so far as we deal with the problem 
of degeneration, this paper may be regarded as an extension, to all types 
of molecules, of the work of Bose? on photons and of Einstein* on mon- 
atomic molecules. 

The Invalidity of Boltzmann's Equation for Distribution with Respect to 
Energy.—Let us first consider a single gaseous substance whose molecules 
may exist in two different states, a and b, which have the energies ¢, and 
¢,. Thus, in addition to its translational state, each molecule may be in 
state a or b. As a simple introduction to our general problem, we shall 
separate these two kinds of freedom and, ignoring the entropy due to 
the translational states for the moment, merely calculate the entropy 
due to the existence of the two states a and b. 

In general, we regard like molecules as inherently indistinguishable from 
one another, but when two molecules are in different states they are by 
that fact distinguished, so that they may be named or numbered by the 
states in which they are. Let us assume for the moment that the gas is 
at such a small concentration that no two molecules are in the same trans- 
lational state; each is therefore distinguishable from all the others. 

If ¢, > ¢«, and M molecules are in the state b, that part of the energy 
of the system which results from the existence of these b-molecules is, 


EE a a ee eee ee ee. ee 





E = M(e, — «). (1) 


We have then, as it were, the problem of distributing M unnumbered 
units of energy among N numbered molecules, no more than one to each. 
For the number of detailed states of the system as a whole we then have 
the following expression, which also, appears in the Fermi statistics, 


is ca N! 
~ (N — MIM!" 


(2) 
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Therefore by our definition of the entrophy, and by the usual approxima- 
tion to Stirling’s theorem, 


S = klogQ = k[NlogN — (N — M)log(N — M) — MlogM]. (3) 
Now we have from the thermodynamic definition of temperature, 
ae i: & 





2 tee ae ee (4) 
T OF «—<€«,0M 
and therefore by differentiating equation 3 with respect to M, 
1 k N-M 
-= lo . 5 
T @-«) M 9s 
Writing M = N,, and N — M = N,, 
N, = Soma 
—=e kT. 6 
N, (6) 


Thus we obtain from our simple rules the familiar equation of Boltzmann, 
but at the same time we see that it is in general incorrect. If the gas had 
been more concentrated, so that sometimes two or more particles occur 
in the same translational state, they then cease to be distinguishable, and 
our whole method of counting Q falls to the ground. The correct method 
of treatment in the general case will be discussed presently. 

Under the same limiting assumptions of a complete and unambiguous 
numbering of the molecules by their translational states, we may study a 
system of molecules having a series of vibrational states, each differing in 
energy from the next by the constant amount ¢9. Here again we are led 
to the classical results. As before, take the energy as E = Me. We 
may think of M unnumbered units of energy distributed among N num- 
bered particles but now without restriction as to the number of units 
that may go to one particle. The number of possible distinct detailed 
states is now given by the equation that also appears in the Bose-Einstein 
statistics, 

_(M+N-1)! 
MI(N — 1)! — 


Ignoring unity in comparison with the large numbers M and N, and pro- 
ceeding as before, 


S = k logQ = k[(M + N)log(M + N) — MlogM — NlogN] (8) 


(7) 





and differentiating with respect to E, we find 


2 N+ mM) 
1k (NM) 9 
5 = Fog (Xt ©) 
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Replacing M by E/eo, 
FM aes , (10) 





an equation which we might also have obtained by finding the relative 
number of molecules in any two states from equation 6. This equation 
has been derived in a similar way by Planck‘ but his method gives no clue 
to the fact that the result is not true, except in the limiting case when the 
molecules are all distinguishable. The general relation which replaces 
that of Boltzmann will be given in equation 37. 

Monatomic Molecules——We may now turn from systems which are 
definitely quantized to a system which is arbitrarily quantized, and 
consider an ideal monatomic gas or a gas in which the mean free path is 
large compared with the dimensions of the enclosure, and in which there 
is no energy except translational. Although most of the results we ob- 
tain are identical with those already obtained by Einstein,’ this discussion 
will serve as an introduction to more complex systems. 

We shall consider N molecules, all of the same. substance, having the 
mass m, with the total energy E, and a volume V. The minimum num- 
ber of states (or cells) C, which we need to consider as available for each 
atom is determined by the sole condition that the three components 
of momentum /,, ~, and p, may independently assume certain discrete 
values. If Ap,, Ap, and Ap, are the intervals between these successive 


values, 
3 


h 
Ap, Ap, Ap; = - (11) 


If therefore we construct a momentum space with p,, p, and p, as 
axes, the number of states C in any region is equal to the volume of the 
region in this space divided by h*/V. 

Consider the region ¢ lying between the two concentric spheres deter- 
mined by the scalar momenta p, and p, + A,p. The number of states 
is then, 

V4xrp; A 


C, = i 


iD. (12) 





Since we are going to regard C; as a large number, A,p must be large 
compared with the quantities such as Ap, which measure the intervals 
between successive states. Since we wish to regard A,p as an infinitesimal, 
this can only be justified by noting: from equation 11 that Ap, Ap, Ap, 
may be made as small as we please by indefinitely increasing the volume. 
Strictly speaking, therefore, the results that we obtain by treating A,p as 
a differential are rigorously proved only for a system of infinite size. 
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If ¢, is the energy of a particle in any one of the states of region /, 


2 
¢ = Pr, Ae = Pt Ap (13) 
2m m 
and, 
a V 20 (2m) '? 


C; 73 


Ave. (14) 





If we have N, particles with momenta lying within this region, our 
problem is to find Q,, the number of ways that the N, particles may be 
distributed among the C, states, remembering that two distributions are 
different only when they differ in the number of particles in one or more 
of the states. It is the problem of distributing NV, unnumbered objects 
among C,; numbered boxes. Hence, 


Q = (N, + C, Rey 1)! 
: Nic, ~ Th 


The entropy of the region is, therefore, 
S; = k log Q, = k[(N, + C,)log(N, + Ct) — N, log N, — Clog C,] (16) 
and 
oS, (™ + 2) 
ame = Kloet f . 17 
ON, e\—y, (17) 


Since we take each Q; as large we may as we have seen in our first 
paper, write for the entropy of the whole system, as the conditions for 


(15) 





equilibrium, 
S=F a x. = ben ’ (18) 
t ON, ON, 
8S = 0. (19) 


Now for any system with a constant energy and number of particles, 
and with entropy at a maximum, the method of variations gives the 
fundamental equation 





oS 

am a : 20 

oN, a + Be, (20) 
where the constants a and B are determined by the total energy and 
number of particles, and § is readily shown from the thermodynamic 


definition of temperature to be, 
1 
B= 7. (21) 


Writing, merely for simplification, 


. log A = -5 (22) 
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and combining all these equations with equation 17, we have, 
ee 1 


ga ee (23) 


1 
—e 
A 

The quantity N,/C,, or the average number of particles in a single state, 
is of fundamental importance. For when this quantity is small compared 
with unity, all methods of counting © are essentially the same, but when 

(in the absence of any special inhibition such as the one suggested by 
Fermi) this average number of particles per state becomes larger, there 
is, in our opinion, which coincides entirely with that of Schrédinger,® 
only one self-consistent method of counting that has ever been proposed. 
This is the method that we have been following, which was first proposed 
by Bose and Einstein. 

Now expressing equation 14 in differential form, by writing N,/ Aj = 
ON /de, and combining with equation 23, 


ON _ V 2m (2m)* ‘Vs 
we Se 
re des 1 
This is Einstein’s generalization of the Maxwell distribution law, which 
must be correct if our fundamental statistical assumption is correct, in 
all cases where there are no special inhibitions. 
For the total number of particles and total energy we have the equations, 


a a eo Tr 7 
N < 7  & Es A OF ads. (25) 
0 Oc 0 Oc 


When we substitute from equation (24), we obtain expressions of a type 
which may be readily integrated by expansion into series, followed by 
integration term by term, using the [ functions. But the resulting series 
is not, as a rule, expressible in terms of known functions. Let us, as an 
abbreviation, use the symbolism, 


on 





(24) 


"- 





ee Fo ee eee (26) 
eg a r r” 
We find, then, 
3/3 

N= v (=) (kT)? .*4A, (27) 

3 Qam\""* 
E=-v(——) (7.4. 28 
: (=) 7) (28) 


So far our results are those of Einstein. 
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By combining equations 14, 16 and 23, we obtain, for the total entropy, 


re) 3/2 7 
S -f kV —_— e/* —+(é- In A) — In(1 — Ae) lie 
0 h® in ‘ T | 
é —_— 


(29) 


>] 


By a very complicated process of integration which, however, does not 
differ essentially from those we have mentioned, 


3/9 
SakV (2) (kT) (2 AA — 14 log A), (30) 
or from equations 27 and 28, 
5E 
S=--— —kNlogA. 31 
3T og (31) 


This is a remarkably simple equation, considering the complicated ones 
from which it is derived. 

The quantity A, which always lies between zero and unity, is smaller 
the smaller the number of particles in a single state. It may be regarded 
as a measure of what has been called the degeneracy or Entartung. These 
words now have their etymological meaning, for when two or more particles 
are in exactly the same state they lose their separate classification. In 
other words, they lose their identity, being no longer distinguishable from 
one another. It is this fact which leads to equations differing from the 
classical equations, but approaching the latter as a limit. 

Thus when A approaches zero, all of the series such as “A, “*A ap- 
proach A. In that case, we obtain the classical result for the energy per 
particle; by dividing equation 28 by equation 27 we find, 

E 


3 
Sac Sag 31 
| ee G1) 


but when A is larger, ‘*A being always smaller than “A, the energy 
per particle is less than the value given by equation 31. 

Some of the consequences of these equations are so curious as to lead 
Einstein to cast doubt upon the validity of his own deductions. If in- 
stead of a single gas, we had considered two different species of molecules, 
each species would have its own particular state. At a given total con- 
centration and a given energy, the number of cases in which two or more 
particles are in the same state would be less for a mixture than for a pure 
gas. This would give a different degree of degeneracy in the two cases 
and therefore, for example, a different temperature, and this would still 
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be true even if the two types of molecules were almost identical (they 
might, for example, have the same mass). Einstein says, “This appears, 
however, to be as good as impossible.” 

Here we must take the opposite view; we must believe that even a 
mixture of two isotopes would have a quite different equation of condition 
from either pure isotope. Indeed, in the cases we are about to consider, 
we must say the same of two kinds of molecules which are exactly alike, 
except that they are in different vibrational or rotational levels. Other- 
wise, we should be obliged to abandon the only self-consistent method 
now available for obtaining thermodynamic information from statistics. 

Molecules Which Possess Internal Energy—The state of a molecule is 
not always completely defined by the values of p,, p, and p,; besides its 
translational energy, the molecule may have various kinds of internal 
energy such as those termed vibrational and rotational. Suppose that 
all of the internal states of a molecule are arranged, as far as possible 
in the order of their energies, and numbered. 

We may, as before, consider a region / in which all the translational 
states have energies lying between ¢, and ¢, + de, and we may use the 
differential method in spite of the fact that such a region comprises a 
large number of individual translational states. 

This cannot, in general, be done for the internal states, since they may 
lie far apart. We shall, therefore, direct our attention to the group of 
molecules, whose number we may call N;,, all of which have just the same 
internal state 7, and translation states lying in the region ¢. 

The total number of individual states thus comprised is given by equa- 
tion 14 which we may now write, 


3/3 \/y 
V 2r (2m) & de, 
ht 
The energy of a particle in one of these states is the sum of its internal 


energy, which is due to its being in the internal state 7, and its translational 
energy, due to its being in translational region ?, 


Gr = & + & (33) 


We may now employ once more the fundamental equations 17 and 
20 and find for the number of particles per individual state, 


Ree 


Cys (34) 
—~e —1 


Ci = (32) 





\2 


co 
x 


A 


Combining the last three equations, we find for the total number of par- 
ticles having a given 7, regardless of their translational energy, 
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7 0 h 3 





JF a (35) 
- eT — 1] 
Integrating just as before, 
3/s _@ 
ys = v (72) en. Mae“, (36) 


It is known that in some cases, as in rotation, we have the phenomenon 
of multiple levels, namely, a group of distinct states which, in the absence 
of certain fields, have the same energy ¢;, By our method of treatment, 
these are handled as quite distinct. If we wish to consider, not individual 
states, but groups of states, then an expression such as is given in equa- 
tions 35 and 36 is merely to be multiplied by 7;, the multiplicity at the 
given energy. 

The ratio of the number of particles in two states, 7 and 7’ (since m is 
the same except for a minute relativity correction), is 


Ny _ “(Ae #) (37) 


Mi % a. i 
"(Ae *f) 


This is a generalized form of the Boltzmann relation. His equation is 
@ i’ 

obtained in the limiting case where Ae *T and Ae *7 are very small. 

Then, 


, =e 
2 =e ‘AT , (38) 





If instead of only one term of each series, we use two, we obtain as a 
second approximation, 


«’ 
Ni _ 3°" iF [ as “oe (a — “| (39) 








We may obtain a better idea of the meaning of A if we follow the com- 
mon practice, and of all the states 7 call the one of smallest energy the 
zero state, and assign it zero energy, thus ¢, = 0. Then we have, 


No =-V (z)" eT)” AA. (40) 
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This is identical with equation 27, and therefore A for the whole system 
is the same as the A for all the particles in the zero state considered as 
simple monatomic gas. If, as often happens, most of the particles are 
in the zero state, they are the ones in which degeneracy is most noticeable. 

All the molecules that are in the internal state 7 may be treated as though 
they alone were present, with the factor of degeneracy, 


A; = Ae #7, (41) 


Substituting this value in equation 28, and adding to the translational 
energy E; (transl.), the internal energy, E; (int.) = Nje;, we find for 
the total energy of all the molecules in the internal state 7, 


3/2 a 
E, = Nye; + : V es, (kT) “(Ae *), (42) 


and from equation 30, 
2am\"/* mis oie Mea a 
S; = kV ag (kT) % (Ae *&?) — “(Ae *T) log (Ae #7) |, (48) 


or from equation 31, 


Et E;(transl.) 


E,(transl.) 
3 T T 


= 5 
S; — kN; log (Ae #7) = 37 


E,(int.) 


— kN; lo A+ 
oe 


(44) 
Since both A and T are everywhere constant, we may easily sum over all 
the states 7 and obtain the total entropy of the gas in terms of the total 
number of molecules, the total translational energy, and the total internal 
energy, 


5 E( 
Sink soe — kN log A. 45 
37 og (45) 


transl.) 4 E(int.) 

T 
This simple equation is of great generality since it is valid for a perfect gas 
with any number of kinds of internal molecular energy. 

When the degeneracy is at a maximum, A = 1, that is, when the con- 
centration is very great or the temperature very low, the entropy assumes 
as simple a form as it does in the corresponding case of black body radia- 
tion (which means radiation with maximum degeneracy). 

Molecules in Electronic States.—We do not wish to enter at length into 
the discussion of the thermodynamics of gases under such conditions that 
the so-called electronic states become important, but such systems afford 
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a very interesting illustration of the principles that we adduced in our 
last paper. 

The work of Bohr showed that the spectrum of the hydrogen atom may 
be accounted for by postulating states of energy, 


es (46) 


n2 


where R is the Rydberg constant and m has integral values. Later, it was 
shown that these states are multiple, the multiplicity increasing with m. 
In order to simplify our illustration, however, we shall ignore this multi- 
plicity, and consider that there is but one state for each value of n. 

If we take the states defined by equation 46 as the ones to be used in 
calculating 2, then, since the number of these states up to a finite energy 
is infinite, every system having enough energy to raise even one atom 
to the limiting energy would have an infinite 2 and therefore an infinite 
entropy. 

Such a result is clearly impossible. Some efforts have been made to 
eliminate the difficulty by assuming that the higher quantum states are 
in some way inhibited owing to a limited volume or a limited free path. 
However, while such an explanation is possible, we shall see that it is 
by no means necessary. 

The first lines of one of the hydrogen series have a natural width which 
is very small compared with the separation of these lines, but as we pro- 
ceed to the higher members, there will be a definite point at which the 
width of the lines becomes just equal to their separation. According to 
our preceding paper this is the point at which we pass suddenly from 
definite quantization to arbitrary quantization, where the states which 
it is necessary to assume would be more or less equally spaced and a finite 
number of them would carry us to the head of the series. 

Below the point at which arbitrary quantization begins, the thermal 
distribution among the several states would be that which is ordinarily 
calculated. Above that point, the quantization should be calculable by 
extrapolation from the quantization of the ionized gas. 

The width of the spectral lines corresponding to low quantum numbers 
can be calculated from the known duration of the excited states. If this 
were constant, the point at which arbitrary quantization begins would be 
about ~ = 400. But by all classical analogies, we should expect the 
duration of the states to diminish rapidly with increasing m, so that it 
seems probable that the number of lines in the Balmer series which have 
actually been observed by astronomers are nearly all that could ever be 
observed, even with indefinite approach to experimental perfection. 

It is not merely that the thermodynamic properties would correspond 











218 PHYSICS: HULL AND LANGMUIR Proc. N. A. S. 


to the assumption of arbitrary quantization, but that by any method of 
observation the spectrum should become continuous at a point below the 
“head of the spectrum’ and that this continuous spectrum should merge 
with the continuous spectrum which is known to exist above the ionizing 
point. 

1 Lewis and Mayer, these ProckEpINGS, 14, 569 (1928); 14, 575 (1928); 15, 127 
(1929). 

2 Bose, Zeits. Physik, 26, 178 (1924). 

* Einstein, Ber., Berlin Akad., 24-25, 261 (1924); 26-27, 3 (1925). 

4 Planck, Ber., Berlin Akad., 26-27, 442 (1925). 

5 Schrédinger, Ber., Berlin Akad., 26-27, 434 (1925). 


CONTROL OF AN ARC DISCHARGE BY MEANS OF A GRID 


By ALBERT W. Hut AND IRvING LANGMUIR 
RESEARCH LABORATORY, GENERAL Evectric Co., SCHENECTADY, N. Y. 


Read before the Academy, November 19, 1928 


The problem of controlling arcs has been under investigation in the re- 
search laboratory for many years. Before describing the results of these 
investigations it may be well to call attention to the difference between 
this problem and that of controlling a pure electron discharge. 

The mechanism by which a grid controls the space charge in a pure 
electron tube is well known. It is also well known that the grid ceases 
to influence the space charge as soon as there is appreciable gas ionization. 
The reason for this failure, however, though often pointed out, is still 
sufficiently new and of sufficient importance to require restatement. 

Figure 1 shows graphically the condition that exists in such a tube con- 
taining mercury vapor. The space between filament and plate is filled 
with a mixture of electrons and positive ions, in nearly equal numbers, 
to which has been given the name plasma. A wire immersed in the 
plasma, at zero potential with respect to it, will absorb every ion and 
electron that strikes it. Since the electrons move about 600 times as fast 
as the ions, 600 times as many electrons will strike the wire as ions. If 
the wire is insulated it must assume such a negative potential that it re- 
ceives equal numbers of electrons and ions, that is, such a potential that 
it repels all but 1 in 600 of the electrons headed for it. 

Suppose that this wire, which we may take to be part of a grid, is made 
still more negative with a view to controlling the current through the tube. 
It will now repel all the electrons headed for it, but will receive all the 
positive ions that fly toward it. There will thus be a region around the 
wire which contains positive ions and no electrons, as shown diagrammati- 
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cally in figure 1. The ions are accelerated as they approach the negative 
wire, and there will exist a potential gradient in this sheath, as we may call 
it, of positive ions, such that the potential is less and less negative as we 
recede from the wire, and at a certain distance is equal to the potential 
of the plasma. This distance we define as the boundary of the sheath. 
Beyond this distance there is no effect due to the potential of the wire. 
The same space charge equations apply to this positive-ion sheath as to 
electron space-charge, except that the mass of the electron is replaced 
by that of the ion. These equations are well known. Hence it is possible, 
knowing the positive-ion current, to calculate the thickness of the sheath 
for any given potential difference between the wire and the surrounding 
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FIGURE 1 


Positive ion sheaths around grid wires in a thermionic tube 
containing gas. 


space. ‘Take, for example, a wire 1 mm. in diameter and a positive-ion 
current density of 10 m. a. per sq. cm., which corresponds to a very weak 
arc discharge. ‘The thickness is found to be 0.0035 cm. with the wire at 
—10 volts, and only 0.02 cm. at —100 volts. This sheath contains the 
whole potential difference between the space and the wire, so that beyond 
the sheath boundary there is no potential gradient due to the wire and 
nothing to indicate the existence of the wire. And since the thickness 
of the sheath is in general very small compared to the distance between 
wires, its only effect is to slightly decrease the cross-section available for 
the arc. This agrees with the experimental evidence that changing the 
potential of the grid has not the slightest effect on the current. 
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It should be pointed out that it is possible, by using closely-spaced grid 
wires and not too great currents, to make the positive-ion sheaths of ad- 
jacent wires touch or overlap. Under these conditions the potential of 
the grid does control the arc current. This kind of control, though re- 
stricted, has many interesting features, and will be reported upon at a 
future time. 

The control to be described in this communication is of a different kind. 
Although it is not possible in general to stop an arc current or to modify 
its value by varying the grid potential, it is possible to prevent it from 
starting. For before current starts 
there are no ions. The electro- 
static fields are the same as in a 
high vacuum tube, and the criterion 
for no current is the same, viz., 
that the grid voltage shall be more 
negative than a definite value. In 
figure 2 let the solid curve repre- 

, ; Gro vareGe _ sent the current as a function of 

noah ore ears grid voltage of a high vacuum tube, 

mares bs ‘ for a constant plate voltage of 100. 

Plate current as a function of grid voltage This tube is shown as having an 

in high vacuum and gas tube, respectively. 

amplification factor 10, which 

means that each volt of grid potential is equivalent to 10 volts plate po- 

tential, so that 10 volts negative grid potential just neutralizes the effect 

of 100 volts positive plate potential and reduces the plate current to zero. 

For all values of grid voltage more negative than 10 the electric field around 
the cathode is negative and no electrons escape. 

The same condition must exist when there is gas in the tube, so long 
as the grid is maintained continuously below 10 volts. As soon as the 
grid potential is raised to —9.9 volts, however, the small electron current, 
which would flow in the absence of gas, now produces ions which form a 
sheath around the grid wires as illustrated in figure 1, and the current in- 
stantly rises to a value limited only by the emission of the filament or the 
circuit resistance. The characteristic of the gas-filled tube, for increas- 
ing grid voltages, is therefore that given by the broken line in figure 2. 
It is the characteristic of an exceedingly sensitive current or voltage relay. 
For example, it is possible by this means to measure the voltage of single 
transients of only a fraction of a volt and less than a microsecond duration. 

The characteristic shown in figure 2 is not reversible. Upon decreasing 
the grid voltage below —10 the current does not fall to zero, but remains 
at its full value. Hence this characteristic represents only a one-way 
control, which is capable of turning on a current, but not turning it off. 
It should be noted, however, that the currents which can thus be turned 
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on may be very large. Hundreds and even thousands of amperes may be 
turned on by a fraction of a volt applied to the grid. 

The next step in grid control is represented by figure 3. Although the 
current cannot be stopped by varying the grid voltage, it can be stopped 
by removing the anode voltage, or reducing it below the ionizing potential. 
There are many ways of accomplishing this, of which two will be discussed. 

The first is to use alternating anode voltage, as shown in figure 3. The 
current then stops at the end of each positive half-cycle of anode voltage. 








GRID VOLTAGE CRITICAL VOLTAGE 


FIGURE 3 
Plate current as function of grid voltage 
in gas tube with alternating anode voltage. FIGURE 4 
If the grid voltage rises above the critical Alternating grid voltage 180 degrees 
value in any cycle, plate current starts and out of phase with anode voltage. No 
flows for the rest of the cycle. plate current. 














FIGURE 5 FIGURE 6 
Phase of grid voltage advanced. Cur- Phase of grid voltage further advanced. 
rent starts just before end of cycle. Current starts earlier in each cycle. 


The positive ions have ample time to diffuse to the walls, so that at the 
beginning of the next positive half-cycle the starting process can be re- 
peated. If the grid is constantly more negative than the critical or cutoff 
value during any cycle, the current will not start in that cycle, if less nega- 
tive current will start. The grid potential thus determines during how 
many of the voltage cycles in a second current will flow, and hence con- 
trols the average current. By increasing the frequency of the alternating 
voltage this average may obviously be made to approximate a continuous 
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current to any desired degree. In figure 3 the critical grid voltage, which 
is a function of anode voltage, is represented by the broken line, and the 
voltage applied to the grid by an arbitrary curve marked grid voltage. 
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FIGURE 7 FIGURE 8 
Example of control of current by phase Method of stopping current with direct 
of grid voltage. Phase is determined by current power supply. Closing switch 
rate of charging of condenser ‘C by photo- makes plate negative for an instant, thus 
electric current. allowing grid to gain control. 
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FIGURE 9 
Time required for deionization of space around grid, as 
function of pressure. 


Current begins in each cycle when the applied grid voltage becomes more 
positive than the critical value, and the total amount of current flowing 
in each cycle is represented by the shaded area. The total current per 
second is the sum of these areas. 
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There is another and better method of controlling the average current 
when alternating plate voltage is used, which consists in allowing some 
current to flow in each cycle, and predetermining, by means of grid voltage, 
exactly how much current shall flow in each cycle. ‘This is accomplished 
by using alternating grid voltage and controlling its phase with respect 
to the anode voltage. It is illustrated in figures 4 to6. In figure 4 an alter- 
nating voltage is applied to the grid which is 180° out of phase with the 
anode voltage and of sufficient magnitude to be always greater than the 
critical voltage. No current flows during the cycle. If the phase of the 
grid voltage is advanced a few degrees, as in figure 5, current will start just 
before the end of each positive half-cycle, and flow for the remainder, as 
shown by the shaded areas. If the grid 
voltage is further advanced, as in figure 
6, the current will start earlier in the cycle ANODE TO GRID 4.2 CM. 

Ol. TEMP 30°C. = 3.7 BARS 
and the average current, represented by 
the sum of the shaded areas, will be pro- 
portionately greater. 

There are many well-known ways of 
controlling voltage phase. A single ex- 
ample is shown in figure 7, in which the 
phase is determined by the rate of charg- 
ing of the grid condenser C by current 
through a photo cell. In this way the 
average current through the tube is made 
to increase continuously and uniformly 
as the illumination is increased. 


























When a direct current source of power FIGURE 10 
is used, it is necessary to employ special Variation of de-ionization time with 
means to reduce the plate voltage in- grid voltage. 


stantaneously to zero, or to such a value 
that ionization will cease, and to maintain it below this value for a time 
sufficient for the ions to diffuse to the electrodes. This may be accom- 
plished by means of a condenser, as shown in figure 8. In order to 
make the process clear, let us assume a line voltage of 250. When the 
tube is carrying current the voltage at its anode, and hence at the conden- 
ser terminal, is about 15. The other condenser terminal is charged through 
the resistance R, to line voltage, 250. If now this terminal is grounded 
by the switch its voltage will fall to zero, a drop of 250 volts. The terminal 
connected to the plate must also suffer instantaneously the same fall, 
which will bring it to —235 volts, from which it will rise gradually to +250 
as the condenser is charged through the resistance R2. If the rate of charg- 
ing is such that there is time for complete de-ionization of the space before 
the anode voltage reaches zero, we have the same condition as with alter- 








oS a nea 





224 PHYSICS: HULL AND LANGMUIR Proc. N. A. S. 


nating anode voltage, and the grid, being now negative, will prevent re- 
starting of the current. 

The time required for de-ionization is thus an important factor. This 
time depends upon the distance between electrodes, the pressure, the 
current density, and the voltage. Figures 9-11 show the values of this time 
for a typical tube. The time is directly proportional to the pressure, in- 
versely proportional to the */. power of the voltage of the grid with re- 
spect to the space, and directly proportional to the 0.7 power of the cur- 
t = 0.0012 pr? 

_" 
servations. This equation is plotted in figure 11 as full lines, with the 
experimental data as points. The grid voltages in figure 11 are measured 


rent. The semi-empirical equation, , appears to fit the ob- 
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FIGURE 11 
De-ionization time as function of current. Solid lines repre- 
sent empirical equation. 


from the cathode, which is 23 volts negative with respect to the space 
around the grid. 

Instead of using a switch to ground the condenser and produce instan- 
taneous negative plate voltage, one may use a second tube, as shown in 
figure 12. The anodes of the tube are connected by a condenser, and the 
current to each is limited by a series resistance. Assuming tube 1 to be 
carrying current and 2 idle, the voltage at the anode of 1 is about 15, while 
that at anode 2 is the full line voltage, 250, for example. If now the grid 
voltage of 2 is raised above the critical value its current will start abruptly 
and its anode voltage will fall suddenly from 250 to 15, a fall of 235 volts. 
The voltage at anode 1 must suffer instantaneously a similar fall, since the 
two are tied together by a condenser which cannot instantaneously dis- 
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charge. This brings the voltage of anode 1 to —220, which momentarily 
stops the current and the production of ions. If the condenser is large 
enough so that the ions in the space between grid and anode have time to 
diffuse to the walls before the anode voltage again becomes positive, the 
grid, now negative, can prevent the current from restarting. The current 
is thus shifted from tube 1 to tube 2. Similarly, by making grid 2 negative 
and grid 1 positive, it may be shifted back to tube 1. By impressing volt- 
age surges, of any kind, first upon one grid and then upon the other, the 
current can be shifted from one tube to the other at will. The grid volt- 
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FIGURE 12 FIGURE 13 
Method of stopping current through Method of shifting current from one 


tube with d. c. supply, by means of tubetoanother. This method may be used 
another tube. to change direct current into alternating. 


age may be impressed by a transformer, as shown in figure 12, either in 
the form of individual surges or continuous sine waves, or in any other 
manner. 

It is obvious that an inductance may be used in place of resistance to 
limit the rate of discharge of the condenser, as shown in figure 13, provided 
the grid impulses are properly timed. And a secondary winding coupled 
to this inductance will receive alternating current. 

This method of controlling the arc current through a tube has many 
applications. An obvious one is the transformation of direct current into 
alternating current, to which the circuit arrangement shown in figure 13 
is directly applicable. There are many scientific problems for which it 
may be used, such as the timing of events and the measurement of short 
time intervals. 
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THE HEAT OF DISSOCIATION OF NITROGEN 


By JOSEPH KAPLAN 
UNIVERSITY OF CALIFORNIA AT Los ANGELES 


Communicated December 16, 1928 


In a recent letter in Nature, Gaviola' presented evidence for a value 
of the heat of dissociation of the nitrogen molecule considerably smaller 
than the heretofore accepted value of about 11.4 volts.2 The purpose 
of this paper is to present arguments which show that the heat of dissocia- 
tion is equal to about 9.0 volts. Recently, R. T. Birge and R. S. Mulli- 
ken,* independently of each other, arrived at the conclusion that the 
heat of dissociation had a value of about 9.5 volts. These arguments 
have been presented in letters to Nature and before discussing the reasons 
for the 9.0 volt value, a brief résumé will be given of the methods by 
which Birge and Mulliken arrived at their conclusions. 

Birge’s arguments are based on the work of Herzberg‘ on the negative 
nitrogen bands. Herzberg, by an extension of the negative bands was 
able to obtain a better curve for the variation of the frequency of vibra- 
tion with the vibrational quantum number than that obtained by Birge 
and Sponer.’ From Herzberg’s data Birge obtained a value of 7.1 volts 
for the heat of dissociation of the normal nitrogen molecule ion. This 
value was obtained by Birge in such a way that it was an upper limit, 
the most probable value being about 0.2 volt less. Using this upper 
limit for D’ and a value of 16.9 volts for the ionization potential of the 
molecule, Birge calculated D, the heat of dissociation of Ne, and obtained 
a value of 9.5 volts. This was done by solving the equation 


In — D’ = D—I,, 


where J, is the ionization potential of the atom. Since the 7.1 value 
is an upper limit, the value of D found above is also an upper limit. It 
is certainly not surprisifig therefore to find that D is less than 9.5 volts. 

Mulliken® obtained his value of 9.5 volts for D from a consideration of 
the products of adiabatic dissociation from the X’ and A’ levels of the 
molecule ion. The best experimental values for the total energies in 
these levels are 26.0 and 24.0 volts, respectively. Assuming, as Herz- 
berg suggested, that the products of dissociation from the X’ level are 
a nitrogen atom in the 7D state and an unexcited atomic ion, Mulliken 
has shown that the value of D is 9.5 volts. This value agrees quite well 
with the 24.0 volt total energy for the A’ level. It is not possible to 
reconcile the old value of D with the above facts relating to the A’ and 
X’ levels. It seems certain therefore that the heat of dissociation of 
N¢z is equal to about 9.5 volts. 
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In the October 6 number of Nature, Herzberg calculated D, using D’ 
= 6.9, lm = 16.7 and J, = 14.5 volts and obtained a value of 9.1 volts. 

The present method for determining D, which is based on ideas entirely 
different from those presented by the above authors, will now be dis- 
cussed. Sponer’s value for the energy of the nitrogen molecule in the D 
level is 14.8 volts. The only bands that have been observed in transitions 
between the D level and other levels of the molecule are the fourth posi- 
tive bands. These bands all correspond to transitions between the Do 
level and the first seven B levels. The absence of transitions from higher 
D levels can be explained if it is assumed that the higher D levels do not 
exist, i.e., that the total energy of the nitrogen molecule in the D level 
is almost the same as the energy in the Dp level. It would be difficult 
to explain the absence of the higher transitions in any other way. The 
value of wo for the Do level is only slightly greater than that in the Bo 
level. This statement is qualitative and it is made on the basis of the dis- 
tribution of intensity among the seven fourth-positive bands.’ This 
intensity is such as to indicate a small relative shift in the two potential 
energy curves of the D and B level. 

It is possible to obtain a fairly good estimate of the total energy in the 
D level if the truth of the above arguments is granted. The vibrational 
energy in centimeters~! of the level having as its vibrational quantum 
number is 

Vy = won(1 — xn). 


If the number of the vibrational level which corresponds to w, = 0 is 
called 0, then the limiting vibrational energy D” is equal to 


D’ = wono(1 — No). 


In the last paragraph the assumption was made that the total energy of 
the D level had a value very nearly equal to 14.8 volts because of the 
non-existence of the Do and higher D levels. Therefore if mo is taken to 
be equal to unity the value of D” so calculated will be a fair estimate of the 
true value of D’. Differentiating v, with respect to n we obtain 


dv, /dn = wo — 2wonx. 
dv,/dn = 0, n = No 


@o — Qwonor = 0 
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Consequently 


D’ = Wolo (1 — St) = = 
2no 2 


“= 


For the D level, 


Ny = Land D’ = > 


Reasons have already been given for thinking that the wo for the D 
level is only slightly greater than wo for the B level. It will be assumed 
here that the value of wo is about 2000 cm.~!. Obviously, the value of 
wo won't affect the discussion very greatly. If the value of 2000 cm.—! 
is used then D’ = 1000 cm.~! = 0.1 volts, and the total energy in the 
D level is equal to 14.8 + 0.1 = 14.9 volts. 

Since the identity of the D level is not certain, Mulliken* has not dis- 
cussed the products of adiabatic dissociation from it. According to the 
arguments of both Birge and Mulliken the heat of dissociation of Ne is 
certainly not much different from 9.5 volts and it is therefore possible 
to predict definitely the products of dissociation from the D level. There 
can be only two possible products of dissociation from this level and 
they are an atom in the 7D state and one in the *P level. The proof 
of this statement follows. The energy of the atom in the *P state is 
given by Compton and Boyce® as 3:56 volts and that in the 7D level has 
been shown to be 2.39 volts. If it were assumed that the products of 
dissociation are two 2D atoms then the total energy in the D level would 
have to equal 9.5 + 4.8 = 14.3 volts and this is, of course, impossible. 
On the other hand, if the products are assumed to be two *P atoms, then 
the total energy based on D = 9.5 volts would be 9.5 + 7.1 = 16.6 volts. 
This, in view of what has been said before, is obviously too high. It 
seems therefore absolutely necessary to say that the products are a *D 
and a *P atom. 

It is possible therefore now to calculate the heat of dissociation D from 
this total energy in the D level, using the fact that adiabatic dissociation 
results in an atom in the ?D state and one in the *P level. The value of 
D so found is 14.9 — 5.9 = 9.0 volts. The total energies in the other 
levels of the molecule will now be discussed in the light of this new value 
for D. 

The total energy in the A level has been found by Birge and Sponer to 
be 11.9 volts.!° If we assume that the products of dissociation are a 
4S and a 2D atom, the total energy should be 11.37 volts. The difference 
between this and the observed value is about the same as the error in most 
of the other cases that are to be discussed. On most of the levels the 
observed value for T. E. is larger than the predicted value. Only in the 
case of the C level, where the extrapolation is non-linear, does the experi- 
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mental value turn out to be lower than the theoretical one. Mulliken" 
notes that in this case the extrapolation should yield a low value. 

The total energy in the ‘‘a’’ level is 14.24 volts and the products of dis- 
sociation according to Mulliken are two *D atoms. The predicted total 
energy is 9.0 + 4.8 = 13.8 volts. It should be noted that the energy in 
the ‘‘a” level is known absolutely from the study of the spectrum of nitro- 
gen in the far ultra-violet’? and is therefore not comparable with the energy 
in the other levels. The error in this case is 0.4 volt and measures the 
error of the linear extrapolation. 

In the B level the total energy is 14.6 volts and the products of dissocia- 
tion are two 7D atoms. The total energy should again be equal to 13.8 
volts, the difference between the experimental and observed values being 
0.8 volt. 

In the C level the total energy was calculated by non-linear extrapolation 
and has been found to be equal to 14.6 volts. The products of dissociation 
are a *D and a *P atom and the calculated total energy is therefore 14.92 
volts. The error in this case is 0.3 volt and is in the expected direction 
since the extrapolation was non-linear, the extrapolated value being 
probably too small. 

In conclusion it may be mentioned that the data on which the new 
value of 9.0 volts has been calculated has been, as far as possible that 
determined by a single experimenter. It seems therefore that it is justi- 
fiable to draw the conclusion that the heat of dissociation of nitrogen is 
9.0 volts within the same accuracy with which Miss Sponer’s results have 
been observed. 

1 Gaviola, Nature, 122, 313, 1928. 

? Birge and Sponer, Phys. Rev., 28, 259, 1926. 

*R. T. Birge, Nature, in press; R. S. Mulliken, Phys. Rev., 32, 761, 1928. 

‘ Herzberg, Ann. Physik, 86, 189, 1928. 

5 Birge and Sponer, loc. cit. 

6 Mulliken, loc. cit. 

7 Strutt, Proc. Roy. Soc., A85, 377, 1917. 

8 Mulliken, Phys. Rev., 32, 186, 1928. 

® Compton and Boyce, Bull. Am. Phys. Soc., Nov. 30, 1928. 

10 Birge and Sponer, loc. cit. (2). 

11 Mulliken, loc. cit. (8). 

12 Sponer, Proc. Nat. Acad, Sci., 13, 100, 1927; Birge and Hopfield, Phys. Rev., 29, 
356A, 1927. 
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EXPONENTIAL YIELD OF POSITIVE IONS IN ARGON 
By KATHARINE B. BLODGETT 
RESEARCH LABORATORY, GENERAL ELectrric Co., SCHENECTADY, N. Y. 


Read before the Academy, November 21, 1928 


When an electric discharge is passed through a gas-filled tube it is often 
a problem to know how much of the current is carried by electrons and 
how much by positive ions. ‘This paper will describe experiments in which 
these currents were measured, and will suggest a mechanism to account 
for the observed relationships. 

The tube used in the experiments had a pure tungsten filament 0.5 mm. 
in diameter as hot cathode, a molybdenum disk-shaped anode 6 mm. 
distant from the cathode, and the tube could be filled with argon at any 
desired pressure. A side bulb sealed to the tube had a glass surface that 
was out of the range of any particles sputtered from the cathode, so that 
it always afforded a clean window through which the filament could be 
observed with a pyrometer. 

Figure 1 shows a plot of 5 volt-ampere characteristics all made with argon 
in the tube at a pressure of 0.26 mm., but made at a series of different 
filament temperatures. A short 0.5 mm. filament has a temperature 
gradient at practically every point along its length so that the temperature 
marked on each curve refers only to the hottest section near the mid- 
point. The plot in each case shows only that part of the characteristic 
which represents the discharge after the arc has struck, that is, after a 
sufficient number of positive ions are flowing to the cathode to neutralize 
the electron space charge at the cathode and permit the full saturation 
electron current to flow to the anode. It is clear from the curves in figure 1 
that with this pressure of gas in the tube the current does not reach a 
limiting value, as it does in a vacuum tube, but on the contrary appears 
to be rising without limit as the voltage is increased. 

The electron current from the filament was calculated for each curve 
by measuring the temperature of the filament at intervals along its length 
with a pyrometer, and calculating the saturation electron emission from 
these temperature data.' It was assumed that the emission from a fila- 
ment in argon can be calculated from the same data that hold for a fila- 
ment in vacuum. A horizontal line on each curve in figure 1 marks the 
electron current, determined in this way. Each mark occurs at a voltage 
between 15 and 19 volts, so that the characteristic in each case passes 
through the current that corresponds:to the electron current at a poten- 
tial close to the ionization potential of argon (15 volts). A large number 
of characteristics, taken under widely different conditions of arc current 
and gas pressure, have all shown this feature. It sometimes happens at 
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very low currents that the arc will not strike until higher voltages are 
applied, and cannot be maintained at voltages less than 30 or 50 volts, 
but in every case in which the arc will run at 15 volts, or less, the current 
in the neighborhood of 15 volts is found to be the same as the electron 
current. 

The curves in figure 1 have nearly the same slope. This is not the case, 
however, when the data are taken at currents above 10 milliamps., for 
then the number of positive ions bombarding the filament at high voltage 
becomes sufficient to heat the filament. Since the filament-heating cur- 
rent is always held constant while the anode voltage is raised, this addi- 
tional heating causes the filament to emit a higher electron current, so 
that the curves bend upward at higher voltages and will eventually rise 
perpendicularly. With a filament at 2000°K., a change of temperature 
of only 10° will raise the electron emission 15%, so that even a small in- 
put of energy into the filament from positive ion bombardment will seri- 
ously affect the slope of the curves. It is necessary, therefore, when ex- 
ploring the region of currents higher than those in figure 1 to measure the 
temperature distribution along the filament for each point on the curve, 
and allow for the changes in electron current when plotting the data. 
When this correction is applied the resulting curves have the same slope 
as those in figure 1. For the purpose of comparing data taken at different 
temperatures it is useful to convert curves of observed current, such as 
those in figure 1, to curves giving the ratio of the observed current to the 
electron current. Since the curves in figure 1 are plotted on semi-logarith- 
mic paper, this amounts in effect to displacing these curves vertically in 
such a way that the lines marking the electron currents all coincide, and 
when this is done the curves nearly coincide throughout. 

Figure 2 shows the ratio of total current to electron current plotted for 
three different pressures of argon. At the higher pressures the current 
rises so rapidly with voltage that at 5.0 cm. pressure and 90 volts the total 
arc current is 13 times the saturation electron emission from the filament. 
Since this means that the arc current in this case is composed of one part 
electron current and 12 parts positive ion current, the question arises of 
how one 90-volt electron can produce 12 positive ions in a gas that has 
an ionization potential of 15 volts. 

In order to see how this rapid multiplication of current may take place, 
and also why it takes place only at higher pressures, let us consider the 
course which the electron follows as it leaves the cathode. It is well- 
known that in a gas tube nearly all the fall of potential takes place close 
to the cathode, and Langmuir and Mott-Smith? have demonstrated that 
this is because the cathode becomes surrounded by a positive ion sheath, 
that is, a region from which electrons are repelled by the negative field of 
the cathode, while at the same time positive ions are attracted in such 
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numbers that they produce a positive space charge which is able to neutral- 
ize the negative field of the filament at a short distance. These authors 
have shown that the thickness of the sheath can be calculated from space- 
charge equations of the positive ion current flowing to the cathode, so that 
in the experiments with this argon tube the sheath thickness was known 
to be usually of the order of magnitude of the filament diameter. The 
remainder of the space between the cathode and anode is a region of almost 
no electric field, that is, a plasma in which there are practically the same 
number of electrons as of positive ions per unit volume. Consider what 
will happen when an electron is emitted from a cathode in a gas in which 
‘the electron mean free path is sufficiently short compared to the sheath 
thickness so that the electron can make several ionizing collisions before 
reaching the boundary of the sheath. The conditions will then be right 
for the Townsend type of ionization by collision which occurs when an 
ionizing agent produces electrons in a region of electric field so that the 
new electrons, being accelerated by the field, become themselves new ion- 
izing agents. If the field is sufficient to allow this process to be repeated 
n times, the yield of positive ions will be (2"—1). Thus if an electron, in 
traversing a sheath across which there is a 90-volt potential drop, is able 
to make 6 ionizing collisions, the theoretical yield of positive ions may be 
as high as 63. Consider on the other hand what will happen if the sheath 
is thin compared to the mean free path, so that the electron reaches the 
sheath edge, or even passes beyond it, before making its first ionizing colli- 
sion. In this case the new electrons produced by collision will be in a 
region of no field so that they will not acquire the energy that they must 
have to ionize, and the only positive ions formed will be those due to the 
original electron which, after falling through 90 volts, will carry with it 
into this region sufficient energy to produce, at the most, 6 positive ions. 

The dotted line (a) in figure 3 represents the maximum theoretical ratio 
of total current to electron current in the first case when ionization takes 
place inside the sheath, and the dotted line (b) represents the case when 
the electron ionizes beyond the sheath. The data for the three current- 
voltage curves in figure 3 were all taken at 2.0 cm. pressure of argon. The 
positive ion currents for the curve at 2088°K. were 10 times those at 1922°, 
and the currents at 2232° were 5 times those at 2088° Now at any given 
voltage the sheath thickness varies inversely with the square root of the 
positive ion current. Consequently the curve at 1922° corresponds to 
the thickest sheath, which would give the greatest opportunity for colli- 
sions within the sheath, and it is seen that this curve has the steepest slope, 
whereas the curve at 50 times higher current has so far departed from this 
steep slope that in the upper region it is approaching the slope of line (b) 
for ionization beyond the sheath. These data, therefore, bear out the 
theory that the relation of the sheath thickness to the distance which an 
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electron must travel from the cathode before it ionizes is a factor of funda- 
mental importance in fixing the magnitude of the positive ion current. 
1 Jones and Langmuir, G. E. Rev., 30;310 (1917). 
? Langmuir and Mott-Smith, G. E. Rev., 27, 449, 538, 616, 762, 810 (1924). 


ON THE RAMAN EFFECT IN DIATOMIC GASES 
By F. RASETTI 


NorMAN BRIDGE LABORATORY OF Puysics, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated February 17, 1929 


The only observation on Raman effect in gases, as far as I know, is 
due to Ramdas.' The substance investigated was ethyl ether vapor, 
and the Raman spectrum was found to be identical with that of the liquid. 

The aim of the present paper is to give some information about the 
Raman spectra of molecules of simple and well-known structure. It 
would have been perhaps more interesting. to observe the effect in atoms, 
but it seems bard to get appropriate metallic vapors of sufficient density, 
and the rare gases, on the other hand, have no energy levels above the 
normal state low enough to give a Raman effect in the optical region. 
But, in the last years, our knowledge of molecular spectra has progressed 
far enough to enable us to verify if the Raman spectrum of a given mole- 
cule is connected with its band spectrum as is to be expected from the 
theory. 

The theory of the Raman effect was given, long before the experimental 
discovery of the phenomenon, in Kramers and Heisenberg’s treatment of 
dispersion? from the point of view of the correspondence principle. In 
this theory some of the fundamental assumptions of the later-developed 
matrix mechanics were already implied, so that the formula given by 
Kramers and Heisenberg for the electric moment of a system perturbed 
by the field of a light wave is exactly the same as it is when derived with 
the methods of wave mechanics. 

It has been shown’ that if we take an atom in the quantum state dis- 
tinguished by the index k (and which we will assume to be the normal 
state), the perturbation due to a light wave of frequency »v gives rise to 
harmonic components in the electric moment of the system having fre- 
quencies given by 


where / represents any other state of the system satisfying only the con- 
dition that the expression above must be positive. 
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Considering the simplest case, in which the electric moment associated 
with the transitions involved in the process has the same direction, x 
say, of the electric vector of the light wave, then the amplitude in the 
direction x of the harmonic component vj, in the electric moment is pro- 
portional to 





‘ 1 1 
ms Xis Xps Vs +y + 4 — ? ’ 
where the quantities 


Xs = i xp.dr 


give the transition probabilities between the states / and s (in particular, 
if x,, = 0, the corresponding transition is ‘‘forbidden’’). 

From this it appears clearly that the condition for observing a transition 
(kl) as a Raman shift is not that this transition should occur in emission 
or absorption (which means x, + 0), but that both the k and / states 
should combine with a third state s (or with more than one), in order 
to give at least one of the products x,,x;, different from zero. 

This means that in an atom with a simple spectrum, such as that of 
any element in the first three columns of the periodic system, practically 
only forbidden transitions can take place in Raman effect; for instance, 
if the normal state is multiple (such as Al, Ga, In, Tl) the transitions 
between the terms of the multiplet. 

All this applies also to the case of a diatomic molecule. Only, if the 
molecule is not polar; such as O2 or Ne, transitions between the different 
vibrational states, without change in the electronic levels, do not give 
rise to radiation because of the absence of an electric moment, and there- 
fore do not occur; though they are not properly ‘‘forbidden” transitions, 
as the integral x, = .f xy,dr does not vanish, the eigenfunctions 
being the same (Hermite’s orthogonal functions in the case of the harmonic 
oscillator) for both the polar and the non-polar molecule. Such transitions 
are expected to take place in the Raman effect, exactly as the ‘‘forbidden” 
transitions in atoms, because of the fact that different vibrational states 
belonging to the same electronic level combine with other states with a 
change in the electronic configuration. 

I have undertaken experiments on diatomic gases to test this conse- 
quence of the theory. 

The gases were used at atmospheric pressure, and contained in a glass 


tube, about 4 cm. in diameter. A quartz mercury arc 12 cm. long was ~ 


placed close to the tube, and a cylindrical aluminum reflector was used to 
concentrate as much light as possible. The light scattered by the gas 
was collected by a lens at the end of the tube and projected on the slit 
of a spectrograph. 
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Considerable care had to be taken to prevent any stray light from falling 
on the slit. Therefore the glass tube was cemented at one end of a wide 
iron tube, about one meter in length, accurately blackened on the inside, 
and provided with many diaphragms with circular holes, in such a way 
as to ensure a perfectly black background; several diaphragms had to be 
placed also between the lens and the illuminated part of the tube. 

The spectrograph employed was very luminous, having a Tessar Zeiss 
1:2.7 aperture lens for the camera. The length of the spectrum on the 
plate was 16 mm. from \ 3650 to 5460. An exposure from 48 to 60 hours 
proved necessary to obtain Raman lines from gases at atmospheric pres- 
sure. 

The gases investigated were Nz, O2 and CO; all of them gave Raman 
lines, excited by \ 4047 or \ 4358 of mercury. 

The wave-lengths were measured by comparison with a copper arc spec- 
trum, a linear interpolation on a frequency scale between neighboring 
copper lines being accurate enough. 

The error in the measurements, evaluated from the constancy of the 
shift for a given transition excited by the two intense mercury lines, was 
always of a few cm.~', so that the data given here may be considered as 
accurate within ten cm.~. 

The results are reported in the following table: 

’R Aug Av “osc. 

Oz 21384 4358 1553 1552 

Ne 22374 4047 2330 2332 
20598 4358 2339 


Co 20781 4358 2156 2133 
22550 4047 2154 


The first column gives the frequencies of the observed Raman lines, the 
second, the mercury exciting lines, the third, the shift in cm.~', the last, 
the energy of the first excited vibrational state of the molecule (in the 
lowest electronic state), as given by R. T. Birge.‘ 

The agreement is very good in the case of Nz and Os, something less 
for CO; but the difference is not much larger than the error in the measure- 
ments, and can be easily accounted for by the structure of the bands due 
to rotational motion. With the dispersion available, the Raman line is 
not resolved, and its position depends on the relative intensities of the 
rotational components. 

In the case of CO, the Raman line corresponds to the known infra-red 
band at 4.66uy. 

The Raman line excited in O2 by \ 4047 could not be observed, as it 
happened to fall very near the over-exposed line \ 4358 of mercury. 

The harmonics of the fundamental frequency were not observed; 
probably their intensity is much smaller, and the exposure was just long 
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enough to bring out the Raman lines given above. Theoretically, they 
should be expected to occur. 

All the observed facts can be accounted for by the quantum theory of 
dispersion, the fact that oscillation frequencies of molecules without an 
electric moment appear in the Raman effect being closely connected with 
the existence of electronic transitions. This appears very clearly if we 
consider how, when the frequency of the exciting light coincides with an 
absorption frequency, the Raman effect becomes the fluorescence observed, 
for instance, by R. W. Wood on iodine vapor. 

Between the two phenomena exists the identical relationship as between 
the Rayleigh non-selective scattering of light by gases, and the resonance 
of mercury or sodium vapor. 

1 Ramdas, Indian J. Phys., 3, 1928. 

? Kramers, H. A., and Heisenberg, W., Z. S. Phys., 31, 681, 1925. 

3 Klein, O., Z. S. Phys., 41, 407, 1927; Dirac, P. A. M., Proc. Roy. Soc., 114, 710, 


1927. 
4 Birge, R. T., Bull. Natl. Res. Council, 11, 132, 232, 1926. 


ADIABATIC EXPANSION IN CASE OF VANISHING 
INCREMENTS. II. 


By Cari BARus 
DEPARTMENT OF PuHysics, BROWN UNIVERSITY 


Communicated January 31, 1929 


1. Introductory Apparatus.—The apparatus is the same as that used 
(figure 3a) in the preceding experiments and consists of an interferometer 
U-gauge, mm’, with one cell (v, capacity 80 cm.*) closed and in com- 
munication with a relatively capacious (555 cm.*) Dewar flask D. The 
other cell, v’, is open to the atmosphere and hence windy days are to be 
avoided. ‘The half-inch clear way stopcock was used at C. It was opened 
and closed suddenly, by reversing the handle, and the number of such 
successive releases of the tension in the Dewar flask is given in the ab- 
scissas of the following diagrams (figures 1-3). Nearly ten of these, 
curiously enough, were admissible without fresh initial exhaustion or 
compression of the air content of the Dewar flask. 

The degree of exhaustion is measured in terms of the parts s of an 
arbitrary scale of the ocular micrometer in the telescope of the U-gauge 
interferometer. One division of s corresponded to 48.1 x 10~* cm. of 
difference of level of the mercury pools of the U-gauge and thus to about 
6.3 X 10-7 of an atmosphere. The initial increment of pressure or ex- 
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haustion was as a rule 120s or 5.8 X 10-* cm. of mercury thus about 
7.6 X 10-* atm. 

2. Data.—Inasmuch as ratios (As/s) only are in question, the ordinates 
of the diagrams are given in s, as observed. In figure 1, the initial ex- 
haustion or compression was of the orders of 110s to 120s; in figure 2 
(drawn to a coarser scale) of the order of 225s; in figure 3 of the order 
of 400s to 450s, thus even exceeding 2 X 10-* cm. of mercury. The 
graphs are numbered (only examples being given) in accordance with 
the data of table 1 in which details the experiments made will be found. 
Expansions of imprisoned air are distinguished from compressions of like 
value by primed numbers. Here so gives the initial exhaustion or com- 
pression; 7 = As/An the initial fall of s per exhaustion or compression 
on release (single sudden opening of the cock), the first or the first 
two releases being taken. Ry = 17/s» is then the corresponding relative 
value, or the fall of s for a single release, per unit of mean air pressure 
under which the release was made. Naturally s refers only to the initial 
increment or decrement of air pressure, the barometric pressure being 
enormous by comparison, as explained in §1. 

If the Dewar flask is initially exhausted; the graphs correspond to the 
full initial pressure S — so where S is the pressure of the atmosphere in 
the same scale. Here s is properly recorded as a decreasing quantity. 

When the gas is initially compressed (S + 5o) the graphs merely record 
the s beginning with zero and ending with so, the pressure being S + so — s. 
The aim was to make the constant s» (not given in the graphs) of about 
the same value for initial compression as for initial exhaustion. 


TABLE 1 
RATES OF APPARENT INSTANTANEOUS EXPANSION AND COMPRESSION 
Sudden Compression Sudden Exhaustion 
NO. So A ro Ro NO. So An ro’ Ro’ 
1 126 
2 114 ; = ie ; sane : Ae oe 
3 117 2 13 0.13 3’ 120 2 25 0.25 
os 103 2 13 0.14 4’ 120 2 22 0.22 
5 120 2 15 0.16 5! 120 2 22 0.22 
6 113 2 15 0.14 6’ 120 2 21 0.21 
7 224 2 26 0.15 “ fe 240 2 40 0.20 
8 224 2 26 0.15 3’ 240 2 45 0.22 
9 405 | 7 0.19 9’ 450 1 90 0.22 
10 410 = oa Are ‘ pas : ae 
Mean Ry = 0.152 0.221 
Mean (Ro’ nase Ro)/Ro = 0.45 


3. Inferences.—As the graphs all depend on the time during which the 
cock is open, and as this time of release is dependent on the quickness 
with which the handle of the cock is reversed between its opposed (closed) 
positions, it is remarkable that the curves come out as smoothly as they 
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do. For the time element here introduced is of very crude constancy. A 
mean time element is thus in question; but the trustworthiness of results 
is enhanced by the generally good agreement of the ratios Ry and others. 
The effective constancy of the time element is improved by the reciprocal 
action of the gauge as explained in the preceding paper,* by which the 
time element becomes of less importance. 

The data of marked interest in the table are the values of Ro=(As/s)/ An 
or the fall of s per unit of s for a single release. When the air in the re- 
ceiver D is originally exhausted, the mean value is Ry = 0.15. When the 
air is originally compressed Ry = 0.22, small increments As being in 
question in both cases. Whether the original s (positive or negative) is 
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increased from 1 to 2 to 4 times makes~no essential difference in the 
result certainly to be recognized, though Ry, may be increasing toward 
R,’ as So increases from 1 to 4 fold. In other words since mean (R,’ 
— R,)/Ro = 0.45 at least, or under otherwise like conditions, the com- 
pressed air escapes from the vessel about half as fast again as the atmos- 
pheric air enters the exhausted vessel. Efflux tends to show a marked 
excess time rate over influx, again ‘under symmetrical conditions. Ex- 
pansion of imprisoned air seems to occur more rapidly than its compression 
throughout. 

These results, which are repeated consistently in the ten series of experi- 
ments made, are not easily explained. Thermodynamically the two 
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processes are essentially symmetrical. It is true that on expansion, the 
gas density and therefore the mass in D diminishes whereas on com- 
pression it augments. But these quantities (10-° to 6 X 107°) are too 
small for consideration in relation to Ry and R,’. Similarly, though the 
elastic changes of the volume of the vessel are in step with the results for 
R, they are also symmetric while the pressure increments to be expected 
are inadequate (2 X 10°, about). The U-gauge consists of two parts 
which again move symmetrically. 

4. Numerical Results —The curves figure 3 are sufficiently definite to 
be tested by the equation s = so 10~™ or again s = so(1 — 10~%), 
from which we get the rates ds/dn = — sclog 10. If this is done, it is found 
that the equation is inadequate in the simple form and must be modified 
tos = so 10° + *”" where cy and a are constants. The rates ds/dn = 
—(co + 2am) log 10 are thus linearly variable with n. The results for 
the ten releases in each case were thus computed (figure 4). 

Considering the method of opening the cock, the mean differences of 
observed and computed s were smaller than would be expected. Only 
in case of the expansion after » = 6, there must have been some un- 
detected irregularity, as is shown by the crosses in figure 3, which give 
the computed curve. This makes the values of the initial compression 
(so — s = 450) too large, implying that the observed rates are too small. 
The correction of the error would therefore accentuate the remarks more 
as to the efflux excess. 

The rates computed as (ds/s)/dn = — (co + 2an)log 10 = —(0.08 + 
0.0035)2.3 and (d(so — s)/(so — s))/dn = —(0.09 + 0.0060n)2.3, if 
n = 1, come out as —0.192 for compression, and as —0.221 for expan- 
sion, respectively, which agree with the mean values Ry and R)’ given in 
the table. 

The same computations made with regard to figure 1 where a single 
original exhaustion or compression is in question deal with an interval 
rather too small for smooth results, unless the fringes are very steady. 
In the case of compression the simple equation s = s)10~™” reproduces 
the results adequately and the replacement of c by co — am was not 
necessary. It indicates merely that c does not increase. 

The lines for — (ds/s)/dn and —(d(so — s)/(so — s))/dn have been 
drawn in figure 4, a, b, c corresponding, respectively, to figures 1, 2, 3. 
Throughout the coefficient for expansion exceeds that for compression, 
but the relation of the lines a, b, c to each other does not follow any ob- 
servable uniform development. The discrepancy between expansion 
and compression lines is largest when the original exhaustion or com- 
pression So is smallest and seems to fade away when sy exceeds 400. Some 
undetected phenomenon probably associated with the failure of the equa- 
tion s = s910°™ and probably of greater menace in the expansions is 
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running in parallel with the apparent simple changes of density of the 
imprisoned air. I have not succeeded in fully searching it out; but the 
evidence so far obtained points to the heat produced on turning the stop- 
cock as the source of the discrepancy. It is improbable, of course, that 
small increments of the expansions and the compressions of a gas should 
not be symmetric. 

* These PROCEEDINGS, 15, 1929, 22-25. 


ON ELECTRONS THAT ARE “PULLED OUT” FROM METALS 
By Epwin H. HA. 
JEFFERSON PuHysiIcAL LABORATORY, HARVARD UNIVERSITY 


Communicated February 7, 1929 


The work! done by Millikan and his co-laborers, Eyring and Lauritsen, 
on the subject of ‘pulling’ electrons out of metals has attracted much 
well-merited attention, but in the present paper I shall undertake to 
discuss certain aspects of this matter which I have not seen presented in 
print. 

It is well known, to the editor of these PROCEEDINGS at least, that I 
am the somewhat persistent advocate of the proposition that we can 
best explain the phenomena of conduction, of thermoelectric action and 
of the Volta effect in metals by supposing that the electrons taking part 
in these phenomena are of two classes, the classes differing from each 
other in respect to the mean kinetic energy and mean potential energy of 
the individual electrons composing them. 

The electrons of one class are supposed to be in a state of higher potential 
energy than those of the other class and to share the energy of thermal 
agitation proper to the temperature of the metal. These I usually call 
the ‘free’ electrons, but I shall in this paper, adopting the name used 
by Millikan, call them ‘‘thermions.”’ Conduction electrons of the other 
class are not supposed to share the energy of thermal agitation, being 
less free from atomic union than the others, and I have usually called 
them the “‘associated”’ electrons. I shall in this paper call them ‘‘valence 
electrons.” 

I have been gradually coming to the conclusion that the thermions 
are very few® compared with the atoms, a conclusion that receives some 
support from the experimental evidence offered by F. Evelyn Colpitts 
in a recently published paper.’ I shall refer to this matter again further 
on. 

Years ago I inferred from conduction and thermo-electric data that in 
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most metals the thermions carry a much smaller fraction of the current 
than do the other class of conduction electrons. For example, I estimated‘ 
that in tungsten at 0°C., the thermions are to be credited with only about 
4% of the current. 

The conclusions drawn from their experiments by Millikan and Eyring 
were in some respects strikingly similar to those indicated above. In the 
abstract of their paper they say, ‘“The lack of dependence of field currents 
upon temperature furnishes strong evidence that most of the conduction 
electrons do not share in the energy of thermal agitation. The thermions, 
however, do share in this energy; they are presumably responsible for 
the Peltier and thermo-electric effects.” 

I called attention® to this similarity of views soon after the paper just 
referred to appeared, but in one respect, apparently, similarity was lack- 
ing. I have not detected in the paper of Millikan and Eyring or in that 
of Millikan and Lauritsen any indication of an opinion that the two 
classes of conduction electrons are in different states of potential energy, 
implying different amounts of work against opposing forces required to 
remove them from the metal. My own theory, on the other hand, as- 
sumes such a difference of potential energy. 

I write as the amount of energy required to bring a conduction electron 
from the non-thermal state to the thermion state 


’ = N+ SRT, (1) 


where X, and s are constants and R is the gas-constant for a single atom. 
At room temperature and thereabout \; is for most metals much smaller 
than sRT. Thus when X’ is expressed in ergs, \, is for tungsten taken 
as 40 X 10~-", while s is taken as 11.4 and R = 1.37 X 10-". I have 
no great confidence in the accuracy of the values assigned to and 
s, and there is risk in extending the application of equation (1) to high 
temperatures, but taking this risk we get for tungsten the following 
values of \’: 


TABLE 1 
r »’ (IN ERGS) »’ (IN ELECTRON-VOLTS) 
300°K. 4.7 X 107" 0.30 
600 °K. 9.4 X 10716 0.59 
900 °K. 14.1 K 1076 0.89 
1000 °K. 15.7 X 107% 0.99 
1100°K. 17.2 X 10-"* 1.08 
1200°K. 18.7 X 10-6 1.18 


From my point of view the value of \’ measures the advantage which 
a thermion has, compared with a conduction electron of the other class 
(a ‘“‘valence’’ electron), in the operation of escape from the metal. If 
thermions and valence electrons were equally numerous within the 
metal, the thermions would furnish more than half of the field current, 
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and the field current, for a given strength of field, would be a function 
of temperature. The same general proposition would hold, of course, if 
d’ represented merely the energy inherent in the gaseous state of the 
thermions. The fact that below 1100°K., and with fields not above a 
certain strength, no such dependence upon temperature was indicated, 
in the experiments of Millikan and Eyring, was doubtless their reason 
for concluding that the thermions are few compared with the other con- 
duction electrons. 

The main purpose of the following discussion is to find, as definitely as 
possible, how few the thermions must be, compared with the valence 
electrons or with the atoms, in order with my values of \’ to contribute 
so little to the field current as the experiments of Millikan and Eyring 
would indicate. 

I shall not be disturbed if I find the thermions relatively very few. 
In a paper already referred to (in footnote 2), applying my ideas to the 
emission data furnished by Dushman and his co-workers, I have obtained 
the following tentative expression for », the number of thermions per 
cu.-cm., in thoriated tungsten: 


n = 3.86 X 10" x T!4, (2) 


Taking N, the number of atoms per cu.-cm. in tungsten, as 6.3 X 107%, 
I get by use of (2) the following estimates: 


TABLE 2 
i n N+ 
300°K. 1.8 X 10% 35 X 108 
900 °K. 5.0 X 10% 13 X 10 
1100°K. 12.2 X 1016 5 X 108 


Miss Colpitts, in the paper already referred to, made use of a method 
having no relation to thermionic emission. Working at room tempera- 
ture she found that her experiments indicated for one tungsten film » = 
1.4 X 10!° and for the other tungsten film » > 9.2 X 10'°, and so of course 
the conclusion to be drawn is somewhat doubtful; but, taken as a whole, 
the evidence she presents is not unfavorable to the view that the estimates 
which I have given above are of the right order of magnitude, and this is 
all I can claim for them at the best. 

In dealing with the experimental data furnished by the papers of Milli- 
kan and his co-workers I shall make much use of certain conclusions ar- 
rived at by Dr. J. R. Oppenheimer, who‘ has applied the ideas and methods 
of wave-mechanics to the phenomena in question. He remarks, ‘‘Any 
field, no matter how weak, will in time dissociate an atom. This is es- 
sentially a consequence of the fact that the motion of the electron is no 
longer absolutely restricted to a region of the dimensions of the Bohr 
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orbit; it will now occasionally, though not very often, be found at points 
much further from the nucleus; and the further it is, the smaller will be 
the field required to insure that it does not return to the nucleus. Since 
the probability that an electron be at a distance R from the nucleus falls 
off exponentially with increasing R, the rate at which the field ionizes the 
atom may be expected to decrease rapidly when the field strength is 
decreased.’”’ Presently he gives an expression, numbered (1), for the 
rate at which a ‘“‘hydrogenic atom with ionizing potential W’’ dissociates 
under the influence of a field F. This expression begins with a numerical 
factor, and this observation is made: ‘‘If the wave-functions of the atom 
are not quite hydrogenic in character, . . . the numerical factor may vary 
by a factor, perhaps, of two.”’ 

In Oppenheimer’s expression (1) W is taken in ergs, while F and the 
electron charge e are taken in e.s.u.; but later W is expressed in electron- 
volts and F in volts percm. With this understanding we have the follow- 
ing statement: “The formula (1) gives in amperes for the current per 
series electron: 


i=3xX 10°F” W” exp. (-10° W” F-).” (3) 


By ‘“‘series electron’ I understand, after consultation with younger 
colleagues, “‘valence electron.” For tungsten, the metal used in the 
experiments under consideration, there are six valence electrons per 
atom. 

Going over carefully the steps needed to proceed from expression (1) 
to equation (3) I have been unable to verify the two constants 3 X 107% 
and 10* given by Oppenheimer in this equation. I get this 


i = 72 X 10°F“ W™* exp. (—6.85 X 107W"? F-), —(3’) 


and my values of the constants have been verified by my colleagues, 
Professor Kemble and Dr. Guillemin. I have made no attempt to allow 
for the departure of the tungsten atom from the hydrogenic type, under- 
standing from Dr. Oppenheimer’s statement, quoted above, that this 
consideration would account for no great variation of the initial numerical 
factor in (3). 

It is to be noted that W means the work done against opposing forces 
in taking a valence electron from the surface layer of atoms out into the 
space surrounding the metal, beyond the reach of the metal’s attraction. 
Oppenheimer, whose paper is dated from the Norman Bridge Laboratory 
at Pasadena, and who doubtless had all the data there available at his 
command, gives 4.7 as an estimate of the value of W, in electron-volts. 
The process by which this estimate was made is not fully shown, and 
probably no great accuracy would be claimed for this particular value 
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of W. But Warner’ gives 4.79 volts as the ‘‘photo-electric work-function”’ 
for tungsten, and Dushman® gives to the constant by in Richardson’s 
thermionic emission equation a value which indicates a work-function of 
4.51 volts for tungsten. There is a strong tendency’ now to identify 
the volt-equivalent of bo with the photo-electric work-function, and so we 
can take Warner’s 4.79 and Dushman’s 4.51 as estimates of the same 
quantity as that denoted by Oppenheimer’s W. The value 4.7 is near 
enough to the mean of the three estimates, and I shall use this value, taking 
it as a constant independent of temperature. 

The F of Millikan and Eyring’s paper is calculated in a purely con- 
ventional way, by means of the rule 


F = 228 


“where ¢ is the difference in potential between the electrodes’; that is, 
between the emitting thoriated tungsten wire and the receiving copper 
cylinder which surrounded it. The diameter of the wire was 0.00123 
cm., that of the copper cylinder was 1.625 cm. F was intended to be the 
field strength at the surface of the wire regarded as uniform, and the 
factor 228 was calculated from the dimensions of the wire and the cylinder 
as given above. ‘The value of F, as thus calculated, which corresponded 
with an emission current of 9.7 X 10-'? ampere was called by Millikan 
and Eyring the ‘‘critical’’ potential gradient. This current of 9.7 X 10~" 
ampere, which I shall call the “‘critical’’ current, is, I believe, the smallest 
one recorded by the authors, and I suppose that it was regarded as the 
smallest current that could be measured in their experiments with a satis- 
factory degree of accuracy; in no other sense, apparently, was it “critical.” 
I shall assume for the purposes of this paper that this current could be 
measured with an accuracy of 1%. If the accuracy of this measurement 
was less than I have taken it to be, the accuracy of any numerical conclu- 
sions that I may arrive at in the following discussion will be reduced 
accordingly. 

The copper-collecting cylinder was about 6 cm. long, and I shall assume 
that all the measured field current came from the 6-cm. length of wire 
contained within the cylinder. But the current came very unequally 
from different points on the surface of this wire. Millikan and Eyring 
say, ‘The ‘field-currents’ in general have their origin in a few minute 
surfacé spots.’’ Microscopic examination of the wire shows, according 
to Oppenheimer, “that at appropriate positions there are small craters, 
and that these are surrounded by protuberances with very small radius 
of curvature; and the dimensions of these points, which should be respon- 
sible for the current, vary from a tenth to a fortieth of those of the wire.” 
“Diffraction, however, makes the precise estimation of the radii of curva- 
ture of the points impossible.’ Oppenheimer gives an estimate of the 
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total area of the effective points as ‘about 10~° that of the wire.” I 
shall presently give a reason, drawn entirely from the data furnished by 
Millikan and Eyring, for believing it to have been very much smaller 
than this. 

The F of Oppenheimer’s equation (3) and of my equation (3’) is the 
potential gradient at these sharp effective points of emission and is there- 
fore very different from, many times greater than, the conventionally cal- 
culated F of Millikan and Eyring’s expression, F = 2289. 

I shall in what follows use first Oppenheimer’s equation (3) and later 
equation (3’). 

With Use of Equation (3).—Case 1. The total area of the emission spots 
1s taken to be approximately 10-5 times the surface area of the 6-cm. length 
of wire: Taking F “critical’’ = 29.3 X 10° volts per cm., I get from 
equation (3) 7 = 4.61 X 10~*! ampere as the contribution per valence 
electron of the surface atoms subject to the potential gradient F. 

The number of atoms per cu. cm. of tungsten being about 6.3 x 107, 
the number of atoms per sq. cm. of the surface will be about 1.58 X 10%. 
The wire being 0.00123 cm. in diameter and, within the receiving cylinder, 
6 cm. long, has a surface area of 0.0232 sq. cm., on which are 3.67 X 10" 
atoms, each with 6 valence electrons, a total of 2.2 X 10'* such electrons. 
If every one of these electrons were subject to the potential gradient F, 
the total field current would be found by multiplying the value of 7 by 
2.2 X 10'4. The product is 1.02 X 10-* ampere. Now this is about 
105,000 times the ‘‘critical’’ field current, 9.7 X 10~!* ampere, which the 
“critical’’ potential gradient, by definition, would produce. We there- 
fore conclude that only 1 in 105,000 of the surface valence electrons is 
active in the emission—that is, only this small fraction of the surface 
atoms are in such a position as to get the full force of the assumed F. 
The active atoms are supposed to be on small protuberances subject to 
a higher value of the field force than the other surface atoms. The ratio 
1 in 105,000 is near enough to the ratio 1 in 100,000 which Oppenheimer 
appears to accept, and so the value F “‘critical’’ = 29.3 X 10° volts per 
cm. may be taken as the one that results from his equation (3) and his 
estimate of the active emitting area. 

Let us now consider how the thermions might be expected to act under 
this same force F. I shall assume that equation (3) can be applied to the 
thermions by merely reducing the value of W. For the valence electrons 
I have, following Oppenheimer, taken W as 4.7 electron volts. I must 
now reduce W by subtracting the value of \’ as given in table 1. With 
T = 300°K., W = 4.7 — 0.3 = 4.4 electron-volts. So for each thermion 
where F = 29.3 X 10° volts per cm. we have by equation (3) 7 = 1.16 
X 10-'* ampere. This is about 25.2 times as great as the 7 found for a 
valence electron. If there were 1 thermion for every 2520 valence elec- 





























Vor. 15, 1929 PHYSICS: E. H. HALL 247 





trons, or for every 420 atoms, the thermions would furnish 1% of the 
“critical’’ current, and if they had done this their effect would perhaps 
have been perceptible even at 300°K. 

But Millikan and Eyring declare that “the ‘critical’ gradients and the 
‘field-currents’ are completely independent of temperature between 300° K. 
and 1000°K.”” Let us try at 900°K. With T = 900°K., W = 4.7 — 
0.89 = 3.81 electron-volts. Then with F “critical” still = 29.3 x 10, 
7 = 4.86 X 10- ampere. This is about 10‘ times the value of 7 found 
for a valence electron with the same value of F. Accordingly, if there 
were 1 thermion for every 10° valence electrons, or for every 167,000 
atoms, the thermions would produce 1% of the field current. Apparently 
they did not make so large a contribution as this. 

At 1100°K., according to table 6 of the -Millikan and Eyring paper, 
the effect of the thermions was on the verge of showing at the “‘critical’’ 
potential gradient and did show in some cases, though not in all. Let 
us, then, try T = 1100°K.., still keeping F ‘“‘critical” = 29.3 X 10° volts 
per cm. We now have W = 4.7 — 1.08 = 3.62 electron-volts andz = 
3.04 X 10~ ampere, which is about 66,000 times the value of z for a 
valence electron. Accordingly, if there were 1 thermion for every 6,600,000 
valence electrons, or for every 1,100,000 atoms, the contribution of the 
thermions to the field current might begin to show. Comparing the 
number 1,100,000 with the value of (V + ) at 1100°K., as given in table 
2, which is 5 X 10®, we see that the two quantities are of the same order 
of magnitude. 

But I am not satisfied with the value of F ‘critical’ that I have been 
using. It was so chosen as to make the active emission area equal to 
about one 10-5 part of the wire surface, microscopic evidence indicating 
some such ratio. But such evidence was, according to Oppenheimer, 
dubious, because diffraction ‘‘makes the precise estimation of the radii 
of curvature of the points impossible.’ A different criterion for the 
estimation of F “‘critical,” and apparently a better one, is furnished by 
the rate of increase of 7 when values of F are increased. Millikan and 
Eyring, in figure 2 and the accompanying text on p. 55 of their paper, 
describe a case in which, starting with critical F and increasing the field 
to a final value about 2.83 times the initial, they observed an increase of 
field current in the ratio of 1 to about 4 X 10’. 

Now if, going back to the valence electrons with W = 4.7 electron 
volts, we increase F from 29.3 X 10° to 2.83 X this value, 82.9 X 10%, 
we get from equation (3) 1 = 3.44 X 10-!! ampere. But this value of 
i is much more than 4 X 10? times the value given by F = 29.3 X 10° 
volts per cm. It is about 7.5 X 10° times that value of 7. 

So I shall now start anew with equation (3). Case 2. “‘Critical” F to 
be determined by the criterion that increasing F to 2.83 times F “‘critical” shall 
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multiply the value of 1 about 4 X 10’ times: With W = 4.7 electron-volts 
and F “critical” = 38.2 xX 10° volts per cm., I get from equation (3) 
4 = 1.63 X 10~" ampere, and with the same W but making F = 2.83 
X 38.2 X 10° = 108.1 X 10° volts per cm.; I get ¢ = 6.54 X 10-!° ampere. 
The second value of 7 is 4.01 X 10" times the first, thus satisfying the 
ratio-of-increase criterion. 

The moderate relative increase of critical F, from 29.3 * 10® to 38.2 
X 10°, that has taken us from the criterion of case 1 to that of case 2, 
changes greatly the estimated total area of the emission spots. In case 
1 this area is about 1 part in 10° of the surface of the wire; in case 2 it 
is about 1 part in 3.7 X 10° of the total surface. This makes the number 
of emitting atoms about 10°. I shall presently consider a much smaller 
estimate of their number. 

Let us now discuss the emission possibilities of the thermions under 
case 2. At 1100°K. we have, with F “critical’” = 38.2 X 10° volts per 
cm., 7 for a thermion = 7.80 X 10~'* amp., which is 4780 times as great 
as the corresponding valence electron 7. If the number of thermions 
were to the number of valence electrons as 1 to 478,000, or to the number 
of atoms as 1 to 79,700, they should furnish 1% of the field current at 
1100°K. with F “critical’’ = 38.2 X 10° volts percm. Apparently they 
did not do so much as this in all cases at 1100° K., though in some cases 
with critical F their effect was visible at this temperature. 

I shall now try equation (3’) instead of equation (3). 

With Use of Equation (3').—Case 1, defined as in use of equation (3): 
Taking W = 4.7 electron-volts and F “‘critical’’ = 18.45 X 10° volts per 
cm., I get = 4.61 X 10-7! ampere. Multiplying this quantity by 2.2 
X 10'4 and dividing by 9.7 X 107'*, as before, I get 105,000, approximately, 
as the ratio of the total wire surface to the emitting area. This comes 
near enough to satisfying the condition of case 1. At 1100°K., with 
F unchanged, the 7 for a thermion would be 8.18 X 107! ampere, and here 
1 thermion to about 17.8 10° valence electrons, or to 2.96 10° atoms, 
would furnish 1% of the field current. 

Case 2, defined as before: ‘Taking W = 4.7 electron volts and F “critical”’ 
as 26.22 X 10° volts per cm., I get 7 = 3.65 X 10~-% ampere.* With 
W unchanged and F multiplied by 2.83 I get 7 = 1.44 X 10%. This 
second i is 3.95 X 10” times the first, which satisfies well enough the 
criterion of case 2. 

With T = 1100°K., F remaining unchanged, the value of 7 per thermion 
is 1.76 X 10-!2 ampere. This is 4800 times as great as the corresponding 
value of 7 for a valence electron. Accordingly the thermions, if there 
were 1 of them for every 480,000 valence electrons, or for every 80,000 
atoms, should furnish 1% of the field current at 1100°K. 

In table 3, I have assembled various estimates of n, the number of 
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thermions per cu. cm., and of (N + n), the number of atoms to a thermion, 
in thoriated tungsten. Of the estimates which I have in the preceding 
pages made from the data of Millikan and Eyring with the help of Oppen- 
heimer’s equation, I shall here retain only those for 1100°K., since in 
dealing with lower temperatures I have merely found lower limits, per- 
haps very remote lower limits, for (V + m). I have no assurance that 
the estimates for 1100°K. are at all close in any one of the several cases 
considered—there are too many elements of uncertainty for any such 
assurance. But I think it is fairly evident that we can harmonize the 
main features of the dual theory of conduction with the experimental 
data of Millikan and Eyring by supposing the thermions to be very few 
compared with the number of atoms, not fewer, however, than other 
evidence would make them. 


TABLE 3 

T SOURCE OF DATA n (N + n) 
300°K. Dushman 1.8 X 10% 35 X 108 
900 °K. Dushman 5.0 X 10% 13 X 106 
1100°K. Dushman 12 xX 10% 5 X 106 
290 °K. (?) Colpitts, film 1 1.4 x 10% 45 X 108 
290°K, Colpitts, film 2 >9 xX10% <ST KIS 
1100°K. M. & E., with eq. (3), case 1 6 xX 10% 1 X 108 
1100°K. M. & E., with eq. (3), case 2 8 xX 10" 8 X 104 
1100°K. M. & E., with eq. (3’), case 1 2 X 1016 3 X 106 
1100°K. M. & E., with eq. (3’), case 2 8 X10" 8 X 104 


General Conception of Conduction.—It may at first glance appear as if 
thermions so few as 1 in 10° or 10° of the atoms could have only a negligible 
effect in conduction. It is to be remembered, however, that, if the ‘‘free 
path,” or rather the free time, of a thermion within a metal is terminated 
only by the capture of the thermion by a metal ion, the conducting power 
of the thermions is, other things being equal, independent of their number. '° 

If the thermions are, let us say, 1 to every million atoms, each thermion 
will be, on the average, about 100 atomic diameters distant from its 
nearest neighbor thermions. If, as I assume, the positive ions within the 
metal are equally numerous, or equally few, with the thermions, we 
naturally ask whether they have any other function in conduction than 
that of terminating the careers of thermions by capture. Have they 
anything important to do with the conduction due to the valence elec- 
trons, which I have in preceding papers called the ‘‘associated” electrons? 
For many years I have been insisting upon the probability that they do 
have such a réle. I am of the same opinion still, but my conception of 
the manner in which they serve has somewhat changed. 

I used to think of ‘‘associated-electron’’ conduction as effected by the 
passage of individual electrons from atoms to adjacent ions, and I was 
much concerned with the kinematics of the operation. The conceptions 
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of wave-mechanics have happily loosened up the relations of valence 
electrons to their nuclei, and with my growing conviction of the relative 
fewness of the thermions, and so of the ions, I have come to think of 
“associated-electron,’’ or valence-electron, conduction as effected by the 
intermittent movements of trains of electrons, each train perhaps extend- 
ing over many atoms. The essential conditions for the movement of a 
train are a potential gradient and a vacant space—that is, an ion—at 
the forward end of the train." 

The conception of electrons as moving in longer or shorter trains through 
a metal is no novelty, but usually, I think, those who have adopted this 
conception have made no provision of terminal facilities for such trains. 
I have made this suggestion before, but I make it now with new emphasis. 

Naturally, heat movements of the atoms would tend to break up the 
trains (Bridgman’s idea) and so, probably, would the presence of atomic 
inequalities such as would exist in alloys, each influence lessening the 
conductivity of the metal. 


* With this value of 7, only 1 in 8.2 X 10° valence electrons on the surface of the 
wire need be active in order to emit the “‘critical’”’ current, 9.7 X 10-12 ampere. The 
number of actively emitting surface atoms on the 6-cm. length of wire is thus reduced 
to about 4400. It appears that in one case (see figure 2, p. 55 of Millikan and Eyring’s 
paper) the total field-current was about 4 X 10~‘ ampere. The question naturally 
arises whether a current of this strength passing out through a total surface area con- 
taining only 4400 atoms would not cause a considerable local heating of the emitting 
surface. 

Assuming as an extreme case that all of the current passes out through a single area, 
I find by calculation that, if the ordinary laws of electric and thermal conduction hold, 
with the ordinary values of the conduction coefficients, the emitting surface, if con- 
tinuous with the general surface of the wire, would be heated only a fraction of a degree 
above the general temperature of the wire, supposed to be at 300°K. If, on the other 
hand, this emitting surface were the outer end of a cylinder projecting a distance of 5 
times its own diameter straight out from the general wire surface, it might be heated 
5 or 6 degrees above the temperature of the wire. 

If the total emission current of 4 X 10~‘ ampere is supposed divided into ten equal 
parts, each going out through an emitting surface containing 440 atoms, this surface 
being the outer end of a cylinder projecting 5 times its own diameter, the surface in 
question will be heated about one-tenth as much as the larger surface previously con- 
sidered. 

It appears, then, that no considerable local heating of the emitting spots would occur 
under any of the conditions contemplated in this paper. On the other hand, spots on 
the receiving copper cylinder opposite the emitting spots on the wire were in some 
cases heated to luminosity by the electrons projected from the wire. 

1 Phys. Rev., 27 (Jan., 1926), pp. 51-67; Proc. Nat. Acad. Sci., 14 (Jan., 1928), pp. 
45-49. 

2See Proc. Nat. Acad. Sci., 13 (1927), p. 323: 

3 J. Franklin Institute, October, 1928. 

4 Proc. Nat. Acad. Sci., 7 (1921), p. 108. 

5 Ibid., 12 (1926), p. 327. 

6 Jbid., 14 (1928), pp. 363-365. 
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7 Ibid., 13 (1927), p. 56. 

8 Phys. Rev., 25 (1925), p. 357. 

® See the paper by Warner, referred to above, and a paper by L. A. Du Bridge in 
the Phys. Rev. for Feb., 1928. See also a paper by the present author in the Proc, 
Nat. Acad. Sci., 15 (1929), p. 126. 

10 FE. H. Hall, Proc. Nat. Acad, Sci., 14 (1928), p. 378. 

11 It seems probable that such a movement would begin at the forward end of the 
“‘train’’ and extend itself backward, somewhat as a railway train starts. 


A PRELIMINARY REPORT ON THE MEASUREMENT OF THE 
Ka LINE OF CARBON 


By Cari E. Howe 
RYERSON PuysIcAL LABORATORY, UNIVERSITY OF CHICAGO 


Communicated February 8, 1929 


In a recent paper by Weatherby! the wave-length of the Ka line of 
carbon is given as 45.4 A, the weighted mean of a number of individual 
measurements that agree very well with each other. This value is quite 
far above that of 44.9 A as determined by Thibaud? and rather close 
to that of 45.5 A, as measured by Dauvillier? using a fatty acid crystal. 
Thibaud? estimates a maximum possible error in his work of 0.2 to 0.3 A. 
On the other hand it is generally recognized that the index of refraction 
of the crystal used by Dauvillier would give a result that is too high. 
Consequently it was thought advisable to re-measure the K a line of carbon. 
This paper is a preliminary report on this measurement. 

The measurements were made by reflection at grazing incidence from a 
ruled grating in a vacuum spectrograph. The spectrograph used was 
that described by Osgood,‘ modified for use with a plane grating. The 
radiation, after passing through a series of slits to collimate the x-rays 
and to reduce the intensity of the optical light from the tungsten filament, 
falls on the grating at a small grazing angle. The size of the ruled sur- 
face (3 cm. X 3 cm.), the width of the x-ray beam (approx. 0.1 mm.) 
and the magnitude of the grazing angle (more than 2° 30’) make it possible 
to have the entire beam “‘reflected’”’ from the ruled surface of the grating. 
The lower half of the radiation is recorded on a photographic plate placed 
fairly close to the grating, the upper half on a plate some distance farther 
back. The direct beam is permitted to fall on the plates by the complete 
removal of the grating from its path. This arrangement permits all 
measurements to be made from the centers of the recorded lines. Fur- 
thermore, the distance between the two plates together with the separa- 
tion of the lines on them permits an accurate calculation of the distance 
from the grating (see the accompanying figure). 
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The two plates are separated about 38 cm. and made parallel to within | 
0.02 mm. The sharpness of the lines together with a large number of 
readings permits in most causes the determination of the position with a 
probable error of less than 0.005 mm. 


_| Plate A 





Plate B 














FIGURE 1 


a’ 
eo 4, tan 20 = “0 oa , tan(2@ + a) = > 


b 
r ins 
The ordinary formula m\ = D(sin 7 — sin 7’) becomes »\ = D[cos 0 
—cos (9 + a)] where @ is the angle of grazing incidence and 6 + a is 
the grazing angle of diffraction. By a simple trigonometric transforma- 
tion this expression reduces to 
6@+a.a 


oe <6 
~~. 





2 
= 2D sin 


The angles 26 + a and 26 are determined from their tangents. By sub- 
traction a is found and substitution is made in the formula. 

Three gratings were used in these measurements: no. 1, 600 lines/mm. ; 
no. 2, 1010.53 lines/mm.; no. 3, 600 lines/mm. Following is a tabulated 
record of the measurements. 


PLATE GRATING 20 DIS. ORDER N 
3 1 ‘ae 36.33 cm. 1-Outside 44.51A 
3 2-0 44.53 
4 1 : td 74.32 1-0 44.46 
4 2-0 44.59 

11 2 8°59’ 25.32 1-0 44.73 
11 4-0 44.50 
11 5-O 44.64 
12 2 8°59’ 63.30 1-0 44.70 
23 3 8°55’ 25.25 1-0 44.59 
23 1-Inside 44.67 
31 3 7°35’ 25.29 1-0 44.70 
35 3 5°35’ 25.34 1-0 44.63 
35 2-0 44.64 
36 3 5°35’ 63.32 1-0 44.56 
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An unweighted mean of the measurements gives a value of 44.60 with 
a probable error not greater than 0.04 A. This value differs from that 
of Thibaud by an amount which is approximately equal to his maximum 
possible error. It is much lower than that of Weatherby. 

In conclusion I wish to thank my friends in Ryerson for the interest 
taken in this work, and especially Prof. A. H. Compton under whose 
direction the work has been carried on. I am also deeply indebted to 
Dr. T. H. Osgood who began this work and who initiated me into it. 

Note Added in Proof. By comparison with the Ka line of aluminum, 
Séderman’ has obtained a value of 44.70 + 0.09 for the Ka line of carbon. 


1 Weatherby, Phys. Rev., 32, 707 (1928). 
2 Thibaud, J. Op. Soc., 17, 145 (1928). 
3 Dauvillier, J. Physique, 8, 1 (1927). 

4 Osgood, Phys. Rev., 30, 567 (1927). 

5 Séderman, Zeit. f. Phys., 52, 795 (1929). 


ON THE IMPERFECTIONS OF CRYSTALS 


By F. Zwickxy 
NorMAN BRIDGE LABORATORY OF Puysics, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated February 13, 1929 


(a) Introductory Remarks.—Imperfections of crystals may be either of 
mechanical or chemical nature. As will appear in the course of this 
investigation, the phenomena related to the chemical imperfections can 
hardly be understood without a preliminary knowledge of the mechanical 
defects of the crystals considered. ‘The aim of the following considerations 
is to give theoretical reasons for the necessity of microscopic cracks even 
in the most perfect crystals. It will also be shown how the characteristics 
of these cracks, such as their size, shape, etc., may in principle be de- 
termined from the properties which characterize the atoms or ions which 
constitute the lattice. The author was led-to his conclusions presented 
here in the course of an attempt to reconcile the discrepancy between the 
theoretical breaking strength of crystals with the observed values for 
this quantity. It is proposed therefore to consider this special problem 
first. 

(b) The Breaking Strength of Crystals—As Polanyi' has shown an 
estimate of the breaking strength Z can be obtained in this simple way. 
Let o be the surface tension or the free energy per unit surface. Then 
the work done in forming two surfaces of 1 cm.? each by breaking a crystal 
is equal to 2c. The forces acting between the two surfaces during the 
process of rupture become negligible after a separation 6 corresponding to 
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a few atomic instances. We have therefore this approximate equation: 
2c¢ = Z.5. Now it has been calculated from the theory of heteropolar 
crystals? that for lattices of the NaCl type we have: o = 0.12e?/R3, 
e = 4.77 es. u. For NaCl itself the lattice constant is R = 5.6 X 107% 
em. Therefore* ¢ = 156 erg/cm.? and Z ~ 2 X 156 X 108 dynes/cm.? 
~ 30,000 kg./cm.? The experimental value is Z = 50 kg./cm.? This 
enormous discrepancy made a rigorous treatment of the problem involved 
a necessity. The author* therefore gave an exact deduction of the break- 
ing strength of rocksalt, making use of the theory of heteropolar crystals 
as developed by Madelung, Kossel, Born and others. The value obtained 
for Z at T = Ois Z = 20,000 kg./cm.?, in qualitative agreement with the 
above estimate. Now the theory of the ionic lattices mentioned before 
is of a very general nature. All the conclusions drawn from it concerning 
the elastic constants, the frequency of the residual rays, the heat of sub- 
limation, etc., have been checked quantitatively by the experiment. The 
total failure to account for the breaking strength therefore confronted 
the theory with a serious difficulty. An attempt to account for the low 
“technical’’ strength by the influence of temperature did not lead to the 
desired result. 

There then followed several very thorough experimental investigations 
of the strength of NaCl crystals. It was proved by different authors‘ 
that under certain special conditions the breaking strength could indeed 
be increased enormously. An NaCl crystal, for instance, whose surface 
is wetted will stand stresses as high as 16,000 kg./cm.? before breaking, 
a value very near to the theoretically deduced strength. On application 
of great thermal gradients to a NaCl sphere stresses were produced in 
the center of the sphere as high as 7000 kg./cm.? without causing rupture. 
The experimental results then seemed to prove that, proper conditions 
being observed, an NaCl crystal will in fact exhibit the strength required 
by the theory. There arose now the question as to the cause of the low 
breaking strength of the crystals in the dry state. 

Griffith® suggested that the existence of fine cracks on the surface of 
crystals might provide a satisfactory explanation. Application of a rela- 
tively low average stress would then indeed produce much higher stresses 
at the edges of the cracks and there cause the rupture. Smekal*‘ elaborated 
this theory. He thinks that actual crystals are never of an ideal geo- 
metrical structure but are made up of a great many microscopic blocks 
leaving cracks and other imperfections between them. We might call 
such a structure a mosaic crystal according to a term used by Ewald.’ 
Smekal correlated this conception with a great many facts such as elec- 
tric conductivity, electric strength, photoelectric absorption, etc. He 
succeeded in determining the size of the elementary blocks in a half em- 
pirical way and he found that they contain in the case of rocksalt and in 
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many other cases about 10,000 atoms each. As Smekals considerations 
enable us to understand a great number of facts it has to be concluded 
that the formation of cracks and blocks cannot be of a purely accidental 
nature. Otherwise, the surface cracks should occur in all kinds of different 
sizes on different samples. But then, it is impossible to understand why 
the technical breaking strength is a relatively well-defined quantity, as 
the different samples of a crystal would break on application of widely 
varied specific stresses. It is possible, indeed, to justify the above views 
on mosaic crystals on a more fundamental basis. 

(c) On the Theory of Mosaic Crystals——In a paper already mentioned, 
the author showed that the maximum lateral contraction y on extenson 
is 5.4% for rocksalt. Calculations of the same kind have recently been 
published by Lennard-Jones* who also found that for the alkalihalides y 
is 6% approximately. This means that considering the ions of one lattice 
plane only, they will try to arrange themselves in a lattice characterized 
by spacings which are 6% smaller than the spacing in the solid crystal. 
This same tendency to contract is also exhibited by the surface planes, 
although in a lesser degree. 

Now the equilibrium configuration of the crystal is one of minimum 
energy.{ But this configuration according to the above is characterized 
by a lattice spacing which is constant in the interior of the crystal and 
decreases toward the surface. Now it is easily seen that the surface 
cannot contract as a whole, but it will do so in patches of a definite average 
size between which cracks will be left open. The linear extension L of 
such a patch corresponds roughly to: L = 100R/y. This is approxi- 
mately 100 A for rocksalt. A much larger extension is excluded because 
of the fact that otherwise ions of a certain sign in the top layer would come 
to lie above ions of the same sign in the lower less contracted layers. 
This then would tend again to increase the energy. Now let ¢.S be the 
decrease in energy corresponding to the contraction of a patch of surface 
S. Notice, also, that the largest cracks will be those which are formed 
adiabatically. Then the new surface S’ (normal to the old) which can 
be formed at the expense of the energy gained by the contraction of S 
will be S’ = S.¢/o. Now « can easily be calculated from the theory of 
the heteropolar crystals. For NaCl the value 220 ergs/cm.? is obtained. 
Although this value is undoubtedly somewhat too large when applied to 
the surface of a crystal, it is nevertheless evident that S’ is approximately 
equal to S. It should be stated also that not the entire surfaces of the 
crack will consume an energy o erg/cm.? for its formation, as at the point 
P» (see figure) the surface energy evidently must become zero. Assum- 
ing then, for instance, an average value o/2 for the surface energy per cm.” 
along PP) we have according to the above this approximate relation for 
the depth of the crack PsP) = 2.L/4 ~ 50 A. The surface of a crystal 
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will then exhibit an appearance as approximately represented in figure 1. 
The dimensions given are those which we would have to expect for the 
largest cracks. 

One further conclusion which can be drawn from the numerical relations 
given above is this: The crystal will not only crack on the surface, but also 
exhibit a great tendency for block formation all through its interior when 
not controlled properly during its growth. This follows from the fact that 
the energy gained from the contraction of a surface is of the same magnitude 
as the surface energy. It should be possible at least partially to avoid the 
formation of a mosaic structure in the crystal by growing it under very 
well-controlled conditions. But it is then to be expected from the theory 
that such crystals will be very unstable against mechanical and thermal 
disturbances. There are a number of facts known which seem to justify 
the above conclusions. A very instructive experiment is this. Certain 
organic crystals grown in the saturated solution are stable as long as they 
are entirely covered by the liquid. They break off immediately, however, 


BFR . 


io 
FIGURE 1 


P,P, ~8A P,P; ~ 100A PoP, ~ 50A. 





as soon as the surface becomes dry. The fact that very thin glass and 
quartz fibers show a high strength just after formation, followed by a 
gradual decrease, might be interpreted in the same way. It is probable, 
also, that the instability of the properties known as cold hardening falls 
into this category. 

It can now be understood that the technical breaking strength will be 
a well-defined quantity. It will be determined by the size and shape 
of the largest cracks, which in turn may be obtained in principle by the 
above considerations. To account for the ratio of the technical to the 
ideal strength will be the aim of a more rigorous treatment of the problem 
involved. 

It is also comprehensible now why an apparently trivial conditioning 
of the surface like wetting rocksalt with water produces the enormous 
strengthening mentioned before. We know from the interionic theory 
of strong electrolytes, that the ions in the saturated solution will influence 
the ions in the solid surface in nearly the same way as the ions arranged 
in lattice planes below the surface. The minimum of the free energy 
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will be obtained with an almost smooth surface. ‘The crystal, therefore, 
is much strengthened as far as rupture is concerned. It should be said, 
however, that the application of great stresses will produce slipping of 
the different blocks in the inner parts of the crystal. This process leads 
to the cold hardening, a phenomenon for which a proper theory is not 
yet given. 

As to the number of particles constituting one of the elementary blocks 
it follows from the dimensions given in figure 1 that this number is ap- 
proximately 10,000. ‘This is exactly the number which Smekal has de- 
duced in a more phenomenological way. 

(d) Generalization of the Theory.—In the case of heteropolar crystals a 
rigorous mathematical treatment of the block formation seems possible. 
We have no quantitative knowledge, however, about the forces of cohesion 
in homopolar crystals. A qualitative generalization of the considerations 
given is nevertheless possible. 

The formation of cracks and blocks in crystals must be a very general 
phenomenon as it is essentially related to the existence of a lateral contrac- 
tion on extension. These qualitative statements can be made on the size 
of the elementary cracks and blocks: 

(1) The surface of the elementary blocks will increase with increasing 
lattice constant. This is due to the fact that a relative shift of half a 
lattice constant (or eventually another simple fraction of R) represents 
a configuration of a relative maximum of the potential energy. 

(2) ‘The surface of the blocks will decrease with an increasing tendency 
to a lateral contraction of the lattice plane forming this surface. -An 
empirical estimate of this tendency can be obtained by studying the geo- 
metrical and energetic aspects of the lateral contraction of the crystal 
on extension. 

(3) The depth of the-largest cracks will essentially be determined by 
the ratio of the energy of contraction ¢ to the surface of energy o. 

(e) Additional Remarks.—In general the cracks will be distributed 
irregularly over the surface of a crystal. It might be expected, however, 
that a surface grown under properly controlled conditions would show 
some regularity in the distribution of the blocks. This then suggests 
the possibility of using such a surface as a grating for soft x-rays, the 
grating constant being for rocksalt approximately 100 A. In certain 
cases like Bi we should expect blocks of considerable size, so that inter- 
ference fringes might eventually be obtained with visible light. 

A very nice proof of the existence of cracks in single crystals of Bi was 
given by Kapitza.* He showed that on application of external pressure, 
the cracks could be closed so that the electric resistance was diminished 
correspondingly. 

Some very interesting conclusions are also suggested in regard to the 
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chemical imperfections of crystals. Suppose that we have two kinds of 
atoms which do not form a mixed crystal. If a crystal which is built 
up by atoms of the first type actually contains atoms of the other kind it 
is evident that these atoms must be absorbed along the cracks. This 
fact combined with our considerations provides a satisfactory explana- 
tion for the following very interesting phenomenon. It is known that 
the amount of absorbed gas in a metal varies in this way with the tem- 
perature. The amount of absorbed gas increases when the temperature 
is raised. At the melting point a sudden jump occurs inasmuch as the 
melt absorbs much more gas than the solid. Sometimes the absorption 
curve shows a minimum considerably below the melting point. From our 
theory it follows that there are two effects which have to be taken into 
account. Firstly, the amount of gas p absorbed per unit surface will 
decrease when the temperature is raised. But the total surface S avail- 
able for absorption on the other hand will increase with 7. In section d;, 
it was shown that S « ¢/o. Now, e will depend very little on tempera- 
ture, as it is essentially related to the lateral contraction. The surface 
tension o, on the other hand, decreases with increasing T so that S increases 
with 7. The total amount of absorbed gas S.p varies with T precisely 
as required by the experiments. To get the number of atoms per ele- 
mentary block just below the melting point, we propose this approximate 
estimate. In the melt, the absorbed gas atom can be associated with 
any one of the N metal atoms present. The surface S,, available for the 
absorption, therefore, is: S,, = 4xr?.N, where r is of the order of the 
diameter of an atom. ‘The total surface S of the elementary blocks in 


4 2/3 
the solid just below the melting point is of the order S = 6N/n G nw) : 


The ratio of the amount of absorbed gas in the solid (G) and in the melt 
(Gi_) is: 
Rin See: S n) * 0.80" 
2\3 

For the absorption of oxygen in silver (10) it is G,,/G ~ 20. ‘Therefore, 
n ~ 15,000, which is of the same order of magnitude as the value required 
by the theory. This satisfactory agreement supports the statement made 
in the beginning, that the chemical impurities in crystals are very essen- 
tially related to the mechanical imperfections of the crystal structure. 

I am indebted to Dr. A. Goetz of this laboratory for having called my 
attention to many experimental facts related to the problem under dis- 
cussion. 


* NaCl in the molten state has at 800°C. a surface tension ¢~100 erg/cm?. 

t See, however, the objections raised by Polanyi in Naturwissenschaften, 1928, p. 1043. 

t Actually the free energy has to be a minimum. It is easy to show, however, that 
the difference between the energy and the free energy is negligible in this case. 
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THE MECHANISM OF SPARK DISCHARGE 


By LEONARD J. NEUMAN 
PHYSICAL LABORATORY, UNIVERSITY OF CALIFORNIA 


Communicated February 11, 1929 


In attempting to account for the ionization required for the initiation 
of a self-sustained electrical discharge, J. S. Townsend! developed a theory 
of the electric spark. While the cumulative ionization of the gas at high 
fields by electrons constitutes part of the mechanism of the spark, Town- 
send has shown that for a self-sustained discharge it is essential that 
electrons be liberated in the neighborhood of the cathode. From his in- 
vestigations on the ionization current between parallel plates produced 
by high electrical fields Townsend concluded that for sufficiently high 
values (X/p), (X = electric field strength, p = gas pressure) the positive 
ions could ionize the gas molecules by impact. As a result of these re- 
searches he developed his well-known equation for spark discharges, At 
that time experiment seemed to indicate that the sparking potential was 
a function of the gas and the electrical field strength and was practically 
independent of the electrode material, consequently Townsend’s equation 
was quite generally accepted. Later Townsend! extended his theory to 
include ionization by positive ions at the surface of the cathode, an effect 
that could conceivably take place at low pressure when the positive ions 
had little chance to strike gas molecules near the cathode. In 1922 
Holst and Oosterhuis? found that the sparking potential in neon (p of 
order of a few mm.) varied in the ratio of 1-3 when the cathode material 
was changed from carbon to rubidium or caesium. Assuming a uniform 
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field between the anode and cathode they showed that the probability 
of positive helium ions acquiring sufficient energy to ionize molecules by 
impact was less than 10~!’. They therefore concluded that the amount 
of ionization produced in the gas by positive ions was negligible and that 
the electrons were liberated from the cathode by some process which was 
not dependent on the energy of impact of the positive ions with the cathode; 
for liberation of electrons by positive ion bombardment of a surface re- 
quires energies of the same order of magnitude as does the ionization of the 
gas molecules. They accordingly set forth the hypothesis that these 
electrons were pulled out of the cathode surface by some mechanism de- 
pendent on the electrical image forces which exist between the approaching 
positive ions and the cathode surface. J. Taylor* repeated, extended and 
confirmed these experiments qualitatively and on the basis of his results 
postulated that electrons were liberated from the cathode photoelectrically, 
according to the suggestion of J. J. Thomson,‘ by the neutralization radia- 
tion produced at the cathode. 

The results of Holst and Oosterhuis and J. Taylor, together with the 
apparently low ionizing efficiency of the positive ions in a gas for low 
fields assumed uniform and the low electron emission from an out-gassed 
metallic surface due to bombardment by positive ions, have led many 
workers in this field to question the validity of Townsend’s theory. Grant- 
ing that a source of electrons exists at or near the cathode surface we may 
classify the possible mechanisms of this source as follows: 

(A) Ionization of the gas by impacts between positive ions and gas 
molecules. (Townsend.') 

(B) Liberation of electrons from the cathode by positive ion bom- 
bardment. This effect depends on the energy of impact and the work 
function of the cathode. (Townsend! and J. J. Thomson.°) 

(C) Liberation of electrons from the cathode or the gas by methods 
which are independent of the energy of motion of the positive ions. 

1. By the photoelectric action of the neutralization radiation on the 
cathode. This radiation is produced by the neutralization of positive 
ions at the cathode or in the body of the gas. (Taylor.*) 

2. By the electrical image forces exerted by the approaching positive 
ions on the cathode surface. (Holst and Oosterhuis.?) 

3. By the thermionic emission produced by the local heat of neutral- 
ization of the positive ions neutralized on the cathode surface. (von 
Hippel.*) 

4. Ionization of the gas by inelastic collisions of the second class or 
by radiation falling on the activated metastable atoms near the cathode 
as a result of electron impacts. (Brode.’) 

5. Ionization of the gas by the photoelectric action of the neutraliza- 
tion radiation on the gas. (Thomson.*) 
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Now these types of mechanisms, (A), (B) and (C) are clearly distinguish- 
able in terms of experimental facts. Type (A) depends on the value of 
X/p and the nature of the gas and is completely independent of the cathode 
surface, Class (B) is a function of the cathode material, the value of 
X/p and the gas used. This type of mechanism is dependent on the gas 
because both the work function and the formation of positive ions are 
functions of the chemical nature of the gas. Type (C) is dependent 
upon the cathode material and the ionization potential of the gas but is 
independent of the velocity of the positive ions. (C) depends upon the 
properties of the cathode surface with regards to its photoelectric char- 
acteristics, thermionic emissivity or work function. The mechanism is 
also dependent upon the gas in that electron emission from the cathode 
is influenced by the gas layers formed on the cathode surface, and by the 
ionization mechanism of the gas by 
electrons. ‘These differences permit us 
clearly to distinguish between mecha- 
nisms (A), (B) and (C) as follows: (A) 
is independent of the cathode material 
and depends on the kinetic energy of 4 
the positive ions whereas (B) depends 
upon both the cathode material and 
the energy of the positive ions. (C) however is independent of the velocity 
or kinetic energy of the positive ions, but is a function of the gas and of 
the cathode material. 

It Was suggested to the writer by A. Joffé and L. B. Loeb that the 
dependence of the mechanisms (A) and (B) on the velocity or energy of 
the positive ions together with the independence of (C) of this factor 
should lead to a critical experiment which could definitely decide which 
mechanism was the essential one in the spark discharge. The discharge 
tube shown in figure 1 was used to distinguish between the two mechanisms 
of ion production. A and G are nickel wire grid electrodes so constructed 
that they can be heated by means of an electrical current. G is placed 
a few mm. in front of the flat end of a glass tube C which can be cooled 
by a stream of water. A tungsten lead is sealed to C. The discharge 
tube T was placed in a furnace while being evacuated and baked at 600°C. 
Carefully prepared sodium was distilled into the hot discharge tube and 
condensed on the cooled tube C while the filaments were kept at dull 
red heat. A clean bright sodium cathode was thus secured. Carefully 
purified argon was then admitted at various pressures. When the spark- 
ing potential is placed across the two clean nickel electrodes A and G 
(A positive) and an auxiliary field of about 10% of the sparking potential 
is placed between the nickel electrode G and the sodium cathode C (G 
positive) the following results should be observed: positive ions generated 
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between A and G are accelerated toward G by a high electric field. Now 
most of these positive ions should pass through the nickel grid and be 
carried toward C by the weak field between G and C. However, in travel- 
ing this distance G-C they make many impacts with neutral molecules and 
thus loose most of their energy of motion gained in the field A—-G before 
reaching C. Any electrons created between G and C or at or on C by 
mechanisms of the (C) class will, however, be carried toward G by the 
weak field between G and C. Now in pure argon (p = 0.1 to 0.3 mm.) 
the sparking potential for a sodium cathode is from 30 to 50 volts less 
than the sparking potential for a clean nickel cathode (see Fig. II). If 
we assume that processes (A) or (B) determine the mechanism of sparking 
we would expect to find the high value of the sparking potential character- 
istic of nickel, since the positive ions of high energy are only found near 
G, and as there is no sodium on G. If, however, processes (C) constitute 
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the mechanism of sparking, then, since it is probable that a large portion 
of the positive ions reach C with only a low energy, and are neutralized 
there, electrons liberated by (C) class mechanism will be carried by the 
weak field, between G and C, into A-G with relatively little loss. We 
should thus expect the sparking potential to be characteristic of a sodium 
cathode, for it is immaterial where the electrons come from as long as they 
are present in the neighborhood of G. 

It was found when the grid G was cooled soon after the sodium had been 
deposited on C that the sparking potential was characteristic of sodium. 
(Curve F, Fig. II). However, if the grid G was heated to dull red heat 
the sparking potential jumped up to the characteristic value for nickel. 
(Curve E, Fig. II). These readings were reversible and could be repro- 
duced at will by condensing sodium on the grid and then vaporizing the 
sodium by heating the nickel electrode. 











VoL. 15, 1929 PHYSICS: L. J. NEUMAN 263 


Upon extending these experiments to higher pressures (about 10 mm.) 
it was observed that the sparking potentials for nickel and sodium cathodes 
were practically the same (i.e., they differed by less than 5 volts which is 
the limit of experimental accuracy) and that the change could be observed 
in a reversible fashion indefinitely by changing the pressure. This fact is 
exceedingly instructive in that it clearly shows that mechanism (A) plays 
a more and more important part in the spark discharge as the pressure 
increases, a fact in complete accord with the observations of previous 
experimenters.'*® Figure III gives a schematic picture of the sparking 
potentials for nickel (curve M) and sodium (curve N) cathodes as a func- 
tion of the pressure in argon.* Whether the gas is ionized by positive 
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ions or whether electrons are liberated from the cathode depends upon the 
relative probability of ionization by these two processes under any given 
set of conditions. The probabilities are a function of the nature of the 
gas and the cathode surface as well as the field strength and pressure of 
the gas. As the pressure increases even with a low probability of gaseous 
ionization by positive ion impact, the increase in the number of molecules 
struck increases the number of electrons from this source relative to the 
number from the cathode which must be decreased by a greater pressure 
and a consequent decrease of high energy positive ions striking the cathode. 
Consequently the ionization of the gas by positive ions becomes more 
important than the liberation of electrons from the cathode. This is 
clearly shown by the curves in figure ITI: 
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From these experimental results one is led to the following conclusions: 

(1) In argon at low pressures the cathode material plays an important 
part in the mechanism of the spark discharge as was previously observed 
by Holst and Oosterhuis, and J. Taylor. 

(2) Under these conditions the principal mechanism by which positive 
ions liberate electrons from the cathode depends upon the velocity or 
energy of impact of the positive ions with the cathode. ~ 

(3) The liberation of electrons from the cathode by the photoelectric 
action of the radiation produced by the neutralization of slowly moving 
positive ions near, at or on the cathode is not the primary mechanism of 
electron production in the spark discharge in argon. ‘This photoelectric 
action may play a minor part. 

(4) Very thin sodium films which distill over from the sodium cathode 
to the grid at room temperature (p = 0.1-0.3 mm.) suffice to reduce the 
sparking potential to the value characteristic of sodium. The formation 
of such alkali films has already been studied by Ives and Johnsrud.!° 
Dull red heat is sufficient to drive the sodium from the nickel grid and 
change the sparking potential to the higher value characteristic of nickel. 

(5) As the pressure increases the predominating mechanisms of elec- 
tron production in the spark discharge change from type (B) to type (A), 
i.e., those in which electrons are liberated from the cathode by the bom- 
bardment of swiftly moving positive ions to those in which electrons are 
generated in the gas by collisions between swiftly moving positive ions and 
neutral molecules. This is in accord with experimental fact, observed by 
many experimenters, that at atmospheric pressure the sparking potential is 
practically independent of the cathode material. 

These results are in conformity with the recent advances in our knowl- 
edge of positive ions. It is true that Jackson! and others have shown that 
the probability of secondary emission from a gas denuded surface due to 
bombardment by positive ions is very low; however, nickel and sodium 
surfaces in argon are not gas denuded. Baerwald!? has found that even 
20-volt positive rays can liberate electrons from metallic surfaces which 
are not out-gassed. Under these conditions the secondary emission is a 
function of the energy of the positive ions and the metallic surface. 

Ionization of gas molecules by positive ions has been a point of much 
controversy. The recent experiments of R. M. Sutton,!* however, in- 
dicate very definitely that 100-volt potassium ions can ionize helium and 
argon in appreciable amounts. In hydrogen this ionization is much less. 

The serious criticisms raised about the mechanisms of the (A) and (B) 
type which these experiments seem to uphold merit some discussion. 
They arise from the fact that to gain enough energy from the field the 
positive ions must fall through potentials of from ten to fifty or more volts 
over relatively short mean free paths (of the order of 10—' to 10~* cm.). 
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If we make the usual assumption that the fields existing at spark over are 
uniform, the acquisition is highly improbable as Holst, Oosterhuis, Taylor 
‘ and Loeb'* have shown. ‘This together with the facts detailed above in 
which it is shown that the mechanisms (A) and (B) are the probable mech- 
anisms lead one to conclude that the fields existing under conditions of 
breakdown are far from uniform before the spark passes. If, then, the 
gratuitous assumption as to uniformity of fields be discarded, the question 
presents no serious discrepancy. The recent measurements of time lag 
in spark discharge observed by Rogowski,!® Torok'® and Beams,!’ where 
intervals of the order of 10~7 seconds are observed with considerable over- 
voltages lead one to the conclusion that their fields are built up in very 
short periods of time. Since the mobilities of the electrons and ions will 
under these conditions barely permit space charges to build up in intervals 
of this order of magnitude, the problem still requires study from the point 
5 of view of space-charge conditions at the time of passage of the spark, and 
the mechanism of their formation. 

The writer wishes to acknowledge his indebtedness to Professors A. 
Joffé and L. B. Loeb for valuable suggestions and to especially express 
his appreciation to Professor Loeb for assistance and his constant en- 
couragement and inspiration. 


* The curves indicated schematically have been observed by the writer. A more 
complete set of data is needed before giving actual results. These are now being ob- 
tained. This work was delayed by lack of argon but it was felt advisabie to publish 
the results to date in the interim. The curve is represented in sections because of the 
enormous range of pressures covered and is typical of curves passing through the mini- 
mum sparking potential upwards, except for the fact previously unknown that the 
differences in sparking potential for different cathodes vanished above 20 mm. pressure. 
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THE CELL DIVISION AT WHICH CROSSING-OVER 
TAKES PLACE 


By JoHn W. GOWEN 


DEPARTMENT OF ANIMAL PATHOLOGY, THE ROCKEFELLER INSTITUTE FOR MEDICAL 
RESEARCH, PRINCETON, N. J. 


Communicated January 29, 1929 


Of cardinal importance to the understanding of segregation and the 
interchange of genes between chromosomes is a knowledge of the par- 
ticular cell division or divisions where crossing-over takes place. 

Two schools of thought have tried to approach this problem. The 
first attempted to explain crossing-over on the basis of segregation of the 
different types of gametes in a cell division shortly after the germ cells 
had separated from the soma cells, and the subsequent differential multi- 
plication of these potentially different cells to the proportions necessary 
to account for the gametic ratios in the progeny. The other school ac- 
counts for the crossing-over ratios on the basis of the direct interchange 
of genes between the homologous chromosomes, the rate of this inter- 
change furnishing the gametic ratios visible in the offspring. The pro- 
ponents of this view early criticized the other or reduplication theory on 
the grounds that, in theory at least, the numbers of cells required in the 
reduplication was immense—beyond the bound of all probability. 

Plough! furnished the first objective contribution to the solution of 
this problem. In this work data were presented to show the effect of 
heat on the rate of crossing-over in Drosphila. Heat was shown to cause 
a marked rise in the crossing-over rate for the factors black, purple, curved 
in the second chromosome. ‘This rise in crossing-over occurred abruptly. 
It had a delayed period, however, before the increase was noted. Thus 
heat-treated females would lay about 150 eggs before any appeared which 
showed the effect of the heat by the increase in the proportion of the eggs 
which showed crossing-over between these factors. The significance of 
the delay became clear when the eggs were counted in the ovaries of two 
flies. About 140 eggs, past the last odgonial division, were believed to 
be present in each fly. From this it was concluded that crossing-over 
probably takes place in the very earliest odcytes. 

The writer’s experiments have now thrown further light on this problem. 
In normal cell division one cell on dividing gives rise to 2 viable cells. 
These 2 cells on dividing give rise to 4 cells—16 cells give rise to 32. This 
process is believed to be carried on in the whole germ-cell tract until the 
maturation divisions. In these maturation divisions in place of two equal 
cells—a larger viable egg cell and a much smaller, nonviable cell—the polar 
body is formed. In the second division of the egg cell the process is re- 











VoL. 15, 1929 GENETICS: J. W. GOWEN 267 


peated—a viable egg cell and a nonviable polar body are formed. Con- 
sideration of these facts presents an hypothesis for the cell division in 
which crossing-over takes place. If crossing-over takes place earlier than 
these last two maturation divisions the expectation would be that both 
cross-overs would be found in the progeny since they would then be in 
separate cells capable of developing eggs. The number of progeny in 
each cross-over class expected from such a single cross-over would depend 
on how much earlier than the last odgonial division the crossing-over of 
the homologous chromosomes occurred. Should it be the last cell di- 
vision, each crossed-over chromosome would be present in each of the 
oégonial cells and therefore would be expected to appear in the progeny 
as one each of the complementary cross-over classes. Were the cross-over 
two-cell divisions back from the maturation division, then two of each 
complementary class of crossed-over chromosomes would be expected in 
the progeny, four such cross-over individuals in all. Crossing-over in the 
oégonial cells thus entails on any view the production of both complemen- 
tary cross-over offspring in the progeny. A female which gives one cross- 
over individual can be expected to give the other complementary cross- 
over offspring in her progeny. The relative degree of crossing-over be- 
tween the given factors does not alter this expectation. It is a constant 
depending simply on the chance of fertilization and subsequent develop- 
ment of the egg. 

If, however, crossing-over takes place in either of the maturation di- 
visions the expectation is quite different. The complimentary cross- 
over goes into the polar body and does not survive. There will, therefore, 
be only one of the cross-over classes present in the progeny from any such 
single cross-over. ‘The balancing-up of the complementary classes of cross- 
overs is brought about by another separate crossing-over when the comple- 
mentary cross-over chances to remain in the egg cell instead of going to 
the polar body. In contradistinction to crossing-over in the odgonial cells 
there is no more likelihood of obtaining the complementary cross-over indi- 
vidual in the progeny of a female showing the first cross-over than the 
chance of getting another independent cross-over. 

The crucial data necessary to determine which of these conditions occurs 
in odgenesis are had when the crossing-over rates are reduced to a mini- 
mum, say one cross-over in a thousand flies with no disturbing lethal or 
other unwanted variables present. Data meeting these conditions are 
available in certain of my stocks. In many thousands of flies involving 
more than 10,000 matings, 16 cross-overs have been noted. These all 
occurred in the progeny of separate females and although the offspring 
of these females were examined as long as possible there was only one 
cross-over progeny present in each culture; the complementary cross-over 
progeny was always notable by its absence. This result is easily explained 
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by the view that crossing-over occurs in the chromosomes as they prepare 
for the first maturation division, for then the expectation would be that 
the complementary class formed by the crossing-over would be absent. 
It must be in the early stages preparatory to this first maturation division 
rather than the second division as crossing-over is known to take place 
in the four-strand stage. The result is accountable on the possibility of 
crossing-over in the oégonial divisions only by a wide stretch of imagina- 
tion encompassing improbable chances. 

These facts are therefore entirely in agreement with Plough’s conclusion 
that the stage when crossing-over takes place is in the very earliest odcytes. 
The proof is furthermore independent of that of Plough and therefore 
freed of certain possible although weak criticisms of technical difficulties 
involved in his work. A further proof of the validity of this approach 
and reasoning would lie in a comparison of the cross-over classes in forms 
which have crossing-over in the male as well as the female. Here, if the 
rate of crossing-over is low enough, it should be possible to show that the 
male had both complementary classes of cross-over in the progeny, whereas 
the female had but one for each cross-over which took place. The results 
of this study bring out one more point which might be mentioned, i.e., 
the equality in number found in each complementary class of a culture 
is always dependent upon the chance occurrence of two separate odcytes 
in which crossing-over has taken place and not, as some have conceived 
it, of a cross-over in one separate cell to form the two reduced gametes. 

1 Plough, Harold H. “The Effect of Temperature on Crossing-Over in Drosophila.”’ 
J. Exp. Zobl., 1917, 24, 147-208. 


AN UNEXPECTED ASSOCIATION OF FACTORS BELONGING TO 
THREE LINKAGE GROUPS IN OENOTHERA AND ITS 
EXPLANATION! 


By GrEorceE H. SHULL 
PRINCETON UNIVERSITY, PRINCETON, NEW JERSEY 


Communicated, January 14, 1929 


It has been shown by extensive genetical studies with the Oenotheras 
that independent inheritance is a relatively rare phenomenon in this 
group of plants. The occurrence of 14 chromosomes in all the known 
species of Oenothera would naturally be assumed to provide a basis for 
the occurrence of seven linkage groups, with independent inheritance 
whenever two pairs of factors belonging to different linkage groups should 
be simultaneously under consideration. Although a considerable number 
of traits have been studied in these species, only three linkage groups have 
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been thus far clearly demonstrated: Besides one large linkage group 
(group I) there are two small groups, II containing as yet only B,b, for 
brevistylis vs. long styles, with another probable pair of this group still 
under investigation, while group III contains the three factors Vv for 
old-gold vs. yellow flower color, S,s, for double vs. single flowers, and 
B,b, for bullata vs. Lamarckiana type of foliage. 

The discovery that in most of the known wild species of Oenothera 
all or nearly all of the chromosomes are connected with one another, end 
to end, to form a multisomal chain or circle which persists right up to the 
early anaphase of the reduction division has made it clear that such an 
arrangement provides an alternative basis of linkage, provided it be 
assumed that the chromosomes of the multisomal circles are alternately 
of paternal and maternal origin, so that when adjacent chromosomes pro- 
ceed to opposite poles, all of the chromosomes of maternal origin go to 
one pole and all the paternal chromosomes collect at the other pole. 

I have recently pointed out” that the arrangement of the chromosomes 
in the multisomal circles is probably not that of an alternation of maternal 
and paternal chromosomes, but one in which the chromosomes keep their 
association in homologous pairs, the orientation of maternal and paternal 
homologs being atrandom. The former arrangement (determinate alterna- 
tion of maternal and paternal chromosomes) would make the multisomal 
circle the basis of a single linkage group, while the latter (indeterminate) ar- 
rangement preserves the chromosome pair as the basis of a linkage group. 

The existence of these two hypotheses as to the bases of linkage groups 
has added interest to the study of linkages in the Oenotheras and makes 
the search for crucial evidence for either of these two hypotheses as against 
the other, one of the most fundamental problems of Oenothera genetics 
and cytology. 

I have already demonstrated that two of the linkage groups are based 
on a single chromosome pair each, since three groups of linked genes are 
found in forms having only one circle of chromosomes. The second and 
third groups are probably the ones thus related to single chromosome pairs, 
as they are independent of group I which is most reasonably assumed to be 
related to the multisomal circle, if any group be so related.* 

Cleland has found that the mutants and hybrid derivatives from species 
which have most or all of their chromosomes included in a single multisomal 
circle, very commonly present varying degrees of resolution of these 
circles, into circles containing smaller numbers of chromosomes, the re- 
maining chromosomes either forming a second multisomal circle in the 
same nucleus, or else (and more commonly) they form a series of separate 
pairs. He has also found cases in which two parents with only paired 
chromosomes have produced offspring in which a portion of the chromo- 
somes form a new multisomal circle. . 
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Such resolution and composition of multisomal circles cannot be without 
fundamental significance for the permanence and consistency of the linkage 
groups if such circles are the cause of the linkages, and it would seem 
reasonable to expect that factors which occur in a single linkage group in a 
biotype with all or nearly all of its chromosomes in a circle will be found 
separated into several independent linkage groups in the derivatives which 
have smaller multisomal circles, and, conversely, groups of genes which are 
independent in forms with small circles or no circles should be fused into a 
single linkage group when multisomal circles are produced de novo, or when 
the number of chromosomes involved in the circles is increased on crossing. 

The discovery of such resolution and composition of linkage groups 
becomes, therefore, a major objective in Oenothera genetics. Until they 
are found, genetical evidence is lacking that the multisomal circles con- 
stitute the physical bases of linkage groups. 

These considerations have kept me on the alert for cases in point for 
several years. I reported to the Joint Genetics Sections last year a case 
which seemed suggestive at first, but which proved on examination to 
involve no new linkage relations, but only a recognition of the fact that 
a single plurigenic trait like old-gold flower color may show linkage with 
as many linkage groups as are represented among its plural determiners.‘ 

During the past summer a case of wholly different nature presented 
itself, in which traits belonging to three different linkage groups were 
associated together much more frequently than was expected on the 
assumption that they were independently inherited. The family in 
question (27152) was segregating for (1) rubricalyx buds versus green- 
hypanthium buds, (2) long styles of Lamarckiana versus short styles of 
brevistyls, and (3) yellow versus old-gold flower color, representing respec- 
tively the first, second and third linkage groups. All of the eight expected 
recombination classes were present, but the ratio deviated markedly in 
several particulars from the expected ratio of equality, as will be shown 
below. 

The history of this family was as follows: The mother of the family 
27152 was an F,; hybrid (26198 (5)) rubricalyx long-styled sulfur, pro- 
duced by mating a long-styled rubricalyx sulfur as seed parent with a 
green-hypanthium outside-in’ gold-center® pollinator. The factorial com- 


position of the egg and sperm which united to produce this mother were 
Oe a 


ot RS a acne “~~ 

R's B, S,V X 1's b, syv, so that her constitution was R's rs B,b, SyV sy. 
The father of family 27152 was 26100(41) with green hypanthia, short 

styles, single flowers and old-gold flower color. His three known linkage 


ss ie. 
groups carried the factors r’s b, S,v. The mating which gave rise to 
family 27152 was, therefore, a trihybrid, R"r'B,b,Vv backcrossed with 
the corresponding triple recessive r'r’b,b,vwv. There was present also 
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the possibility of producing by the process of crossing-over, sulfur-colored 
flowers, double flowers and outside-in flowers, but none of these last three 
types actually appeared, although one plant which developed so late that 
only one well-developed bud was produced may have been a sulfur, as 
this bud, when opened artificially, was found to be a long-styled rubricalyx, 
and seemingly with sulfur flowers; but in buds younger than one day 
previous to anthesis the flower-color determination is too precarious to 
justify any confidence in this apparent instance of crossing-over. In 
addition to the above-mentioned segregating factors, which were pre- 
dictable in advance, there was present also a heterozygous condition for 
green versus chloralbinic foliage which showed segregation in the family 
under discussion (27152) into 114 green, 152 more or less chloralbinic. 
Of the chloralbinic offspring 91 died in the seed pan, 17 were lost in the 
pots and 44 were finally set in the field. Of the latter 7 died and 18 failed 
to bloom. All of the albinos which bloomed had rubricalyx buds and all 
had yellow flowers except the one mentioned above which had rubricalyx 
buds but may have had sulfur-colored flowers. The segregates in family 
27152 may be tabulated as follows: 


ALBINIC SERIES GREEN SERIES 
11 Rubricalyx long-styled yellow 3 (+15) Rubricalyx long-styled yellow 
8 Rubricalyx brevistylis yellow 3 Rubricalyx brevistylis yellow 
16 Rosettes only 2 Rubricalyx long-styled old-gold 
2 Failed to bloom (of which one 1 Rubricalyx brevistylis ,old-gold 
may have been a rubricalyx 13 Green-hypanthium long-styled yellow 
long-styled sulfur) 22 Green-hypanthium brevistylis yellow 
7 Missing (died in field) 16 Green-hypanthium Iong-styled old-gold 


23 Green-hypanthium brevistylis old-gold 
6 Rosettes (4 probably triploid) 


Several features of these results merit special consideration: (a) the 
absence of old-gold-flowered plants in the albinic series, (b) the small 
proportion of the green series showing rubricalyx buds except in the long- 
styled yellow class, and (c) the considerable excess of rubricalyx plants with 
long-styled yellow flowers in comparison with the other classes. Only 
the last of these phenomena has been chosen for special consideration 
here, but the other two points may be briefly considered first. 

The absence of old-gold among the albinos is undoubtedly due solely to 
the retarding effect of old-gold compared with yellow on the process of 
development, as shown in a previous paper.’ The low development of 
chlorophyll made it difficult to get any of this class to come to bloom. 
Those that did bloom began late in the season, and it is reasonably certain 
that owing to the additional handicap of the old-gold factor those albinos 
which possessed this factor failed to bloom. 

The small number of rubricalyx plants in the green series, and the fact 
that all of the albinic series have rubricalyx buds is reasonably explained 
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as indicating a linkage between rubricalyx and albinism, the rubricalyx 
among the green plants being thus of the nature of cross-overs. The fact 
that rubricalyx old-golds occur among the cross-over plants leaves no 
doubt as to their potential presence in the albinic series. 

The condition which attracted special attention to this family, however, 
was the fact that of the 24 green plants which had rubricalyx buds 18 had 
long-styled yellow flowers, while there was no such excess in the green- 
hypanthium group and there was no ground for assuming a differential 
elimination among the green plants, such as is so obvious in the albinic 
group. Only six plants (one-fourth of the total) should have been rubri- 
calyx long-styled yellows on the assumption of independence of these 
three traits. This suggestion that factors which had been found in- 
dependently assorted in other matings were remaining together in this 
family in an undue proportion of the offspring led me to plan at once 
a genetical analysis of a goodly number of the plants which showed the 
unexpected association in question, but on preparing to make the necessary 
backcrosses to carry out such analysis, I discovered at once another 
solution, which placed a wholly different face on the problem, for on 
careful examination it was seen that fifteen of the eighteen rubricalyx long- 
styled yellows in the green class were triploids, as indicated by large, rather 
heavy foliage, large buds, short, relatively heavy capsules and scanty 
pollen, of which always a greater or less proportion of the grains were 
four-sided and four-angled. The fact that all of these triploids had the 
known dominant factors of the mother R’B,VS,, together with the two 
dominant factors SS, of the pollen-parent, shows that they resulted from 
eggs which had undergone no reduction division, and which, therefore, 
had 14 chromosomes carrying all of the factors brought in by the con- 
tributing egg and sperm of the preceding generation. These unreduced 
14-chromosome eggs fertilized by normal 7-chromosome sperms produced 
these triploid plants which must therefore have the genotypic formula: 


a nr 


SN ge 
R's r's rs+B,bb,+SyV Sv syv 


Besides these 15 plants which were demonstrably triploid, four of the six 
green plants which remained rosettes were almost certainly triploid. 
This makes more than 18 per cent of the green plants triploid. 

Triploid mutants, usually referred to as ‘‘mut. semi-gigas’’ or “‘heroes,”’ 
are of frequent occurrence in the Oenotheras as also in many other plants, 
and have been reported many times by Oenothera geneticists. They have 
occurred every year in my cultures, but until this year they have never 
been present in numbers sufficient to constitute a considerable percentage 
of any family or to interfere in any essential way with the inheritance 
ratios. Since these plants arose from eggs in whose production no segre- 
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gation had been possible, because no reduction division had taken place, 
it is clear that they should be simply set aside as having no relevancy for 
the determination of problems of linkage. When these are deducted 
from the table on page 271, giving the segregations in family 27152, it 
is seen that there remains only a slight suggestion of linkage between 
brevistylis and green-hypanthium buds, and this of so low a grade that 
it lies within the possibility of a simple chance deviation from the ex- 
pected 1:1:1:1 ratio. 

In most cases of triploidy it is not so easy as in the present case to see 
just what process has been involved in their production. It seems prob- 
able that in most cases the origin of triploids is the same as here demon- 
strated, namely, the fertilization of unreduced eggs by reduced sperms, 
for when extra chromosomes are present in the pollen they seem to present 
obstacles to the functioning of the pollen which possesses them. On the 
other hand, it should not be overlooked that 14-chromosome pollen does 
effect fertilization in Oenothera gigas, as also when gigas pollen is applied 
to the stigmas of Oe. Lamarckiana and other diploid species. In these 
cases there is no competition between 14-chromosome and 7-chromosome 
pollen, and it is probable that this lack of competition is the crucial differ- 
ence which gives 14-chromosome sperms an opportunity to effect pollina- 
tion in the one case and not in the other. 

The only suggestion yet available as to what may have prevented 
reduction in so large a proportion of the eggs of the mother of family 27152, 
is the discovery of Blakeslee and Cartledge,* that low temperature may 
prevent or modify the reduction division in Datura. It seems a little 
doubtful that this explanation can have a bearing on the present case, 
however, since family 27152 resulted from a series of pollinations which 
were made between August 23 and 30, when there were no low tempera- 
tures obtaining, comparable with those which were found to be effective 
in Datura. Also, it may be recalled that Blakeslee and Cartledge’s 
studies on Datura related to microsporogenesis and conditions producing 
non-reduction in megasporogenesis may be entirely different. 

The clearing up of this case of unexpected association of factors belong- 
ing to three linkage groups, together with the simple explanation found 
last year for an apparent linkage between factors belonging to two different 
linkage groups, may suggest to some reader that I deem it demonstrated 
that the linkage groups in Oenothera are not subject to resolution and 
composition, in the manner pointed out above as one of the expected 
results if the delayed separation of telosynaptic non-homologous chromo- 
somes is the causal mechanism of linkage in the Oenotheras. It is im- 
possible to prove a universal negative, and although two cases have been 
explained without resort to a change in the inclusiveness of the known 
linkage groups, the next case which arises may have no other reasonable 
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explanation than on the resolution and composition of the chromosome 
circles. Not much work has yet been done along this line, and it is hoped 
that more workers will take up the study of Oenothera linkages in the hope 
of an early decision as to whether multisomal circles have any significance 
for the linkage of factors in this group of plants. 

1 The experiments on which this paper is based have been supported in part by 
grants from the American Association for the Advancement of Science and the BACHE 
Funp of the NATIONAL ACADEMY OF SCIENCES. 

? Shull, G. H., ““Oenothera Cytology in Relation to Genetics,’’ Amer. Nat., 62, 97-114, 
1928. 

3 Shull, G. H., “Crossing-Over in the Third Linkage Group in Oenothera,”’ Proc. 
Nat. Acad. Sci., 13, 21-24, 1927. 

4 Shull, G. H., “Linkage with Crossing-Over between Rubricalyx Buds and Old- 
Gold Flower Color in Oenothera,” Ibid., 14, 147-149, 1928. 

5 Outside-in flowers are the joint product of two recessive factors, b-b; for short 
styles and spsp for double flowers. See Shull, G. H., “The Outside-in Oenothera 
Flower, a New Morphological Type Produced by the Interaction of Two Recessive 
Mendelian Factors,’”’ Ibid., 14, 142-146, 1928. 

6 Gold-center flowers represent the double recessive of the two dominant comple- 
mentary factors whose coéperation produces the typical yellow color of the wild Oeno- 
thera flowers. Yellow has the genic composition SSVV and gold-center is ssvv. See 
footnote 7 for reference. 

7 Shull, G. H., “ ‘Old-Gold’ Flower Color, the Second Case of Independent Inheri- 
tance in Oenothera,’’ Genetics, 11, 201-234, 1926. 

8 Blakeslee, A. F., and Cartledge, J. L., “Sterility of Pollen in Datura,’ Mem. Hort. 
Soc., New York, 3, 305-312, 1927. 


AN INHERENT DEFECT IN THE THEORY THAT GROWTH RATE 
IS CONTROLLED BY AN AUTOCATALYTIC PROCESS 


By GEorGE D. SNELL 
Bussey INSTITUTION, HARVARD UNIVERSITY 


Communicated January 30, 1929 


Since the increase in size of most plants and animals during development 
is due principally to the synthesis of new protoplasm and hence is essen- 
tially a chemical process, it would seem that some inference as to the nature 
of this process might be made by a comparison of the rate of growth with 
the rate of the various types of chemical reactions. Such a comparison 
was first made by Jacques Loeb in 1906, and enlarged upon by him in later 
papers. From a study of the early cleavage stages of fertilized eggs he 
concluded that not only the number of nuclei, but also the total mass of 
nuclear material, is doubled with each successive cell division.! Start- 
ing with one nucleus following fertilization, there are two nuclei as large 
as the first after the first cleavage, four after the second cleavage, eight 





VoL. 15, 1929 PHYSIOLOGY: G. D. SNELL 275 


after the third, and soon. Moreover, the length of time elapsing between 
successive cleavages is practically constant provided the temperature be 
constant. Hence ‘‘the mass of the nuclear material during successive 
cell divisions (at least during the early stages) increases in geometric pro- 
gression.” 

Since among chemical processes the only ones which show this same 
geometric increase are those in which the end product is itself a catalyst 
for the reaction, Loeb concluded that the growth of nuclei must be a proc- 
ess of this type. He says: 

“For this peculiar condition there is, so it seems to me, only one expla- 
nation: namely, that the nuclear material, the end product, acts itself 
as catalyst for the synthesis; and therefore that the synthesis of nuclear 
material is an autocatalytic reaction.”’ 

This type of reaction alone shows the geometric increase in rate char- 
acteristic of the growth of nuclei, and for the reason that it is the only one 
in which the amount of catalyst increases as the reaction progresses. 
Moreover, as Loeb points out, the continuity of the germ plasm is most 
easily explained by assuming that the nuclei (or better perhaps the indi- 
vidual genes as suggested by Muller, 1922) are catalysts for their own 
synthesis. It therefore seems probable that autocatalytic processes play 
an important part in the accelerating nuclear growth found in the early 
stages of cleavage of many organisms. Hence, too, it is probable that 
they are involved in the steady increase in slope of the growth curve which 
occurs during the first half of every growth cycle. Such an increase in 
rate as is there observed is strong evidence that if catalytic processes are 
involved at all, and they surely are, they must be autocatalytic. 

The question then naturally arises, if the upward bending of the growth 
curve during the first half of the growth cycle is due to the action of auto- 
catalytic processes, may not also the gradual decrease in slope and ulti- 
mate flattening which occur in the second half of the cycle be due to the 
same cause? Cannot the retardation and ultimate cessation of growth 
with the approach of maturity, a retardation and ultimate cessation which 
is also to be observed in autocatalytic processes, be explained on the ground 
that growth is regulated throughout by such processes? 

This question has been answered in the affirmative by both Robertson 
and Crozier. They hold that growth throughout the whole life cycle, 
both in the period of acceleration and in the period of retardation, is regu- 
lated by one, or at most a few, monomolecular autocatalytic reactions. 
The theories which they advance rest primarily on the similarity which 
they (and also Ostwald, Reed, and others) have demonstrated to exist 
between the rate of organic growth and the rate of autocatalytic reactions 
as carried out in the laboratory. Both processes start slowly, gradually 
accelerate to a maximum ‘velocity, and finally slow down again until they 
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cease altogether. The similarity is so close that the equation which 
describes the one will accurately describe the other.? This fact they take 
to be evidence that the rate of growth of most plants and animals is con- 
trolled by a single or ‘‘master’’ autocatalytic reaction. 

Robertson has worked out the theory in detail along the following lines. 
The master reaction results, so it is assumed, in the production of an end 
product, let us call it X, which is essential to growth. The reaction goes 
very slowly, more slowly in fact than any other involved in the synthesis 
of new protoplasm, with the result that all the others are kept in check 
by it, and the growth of the whole organism waits upon this one process. 
Hence the size of the organism is proportional to the amount of X that 
has been produced, and the equation for an autocatalyzed monomolecular 
reaction which gives the amount of X at any given time also gives the 
size of the organism at that time. 

Robertson develops this equation in the following manner. He assumes 
that the end product X is derived from a substrate whose concentration 
in the organism remains constant throughout growth. He assumes further 
that the reaction is a reversible one and that it is finally brought to a stop by 
a steady increase in the concentration of X, which accelerates the reverse 
reaction to a point where it equals in velocity the forward one. Being 
autocatalytic, both the forward and reverse reactions are catalyzed by X. 
LetA = the concentration of the substrate (a constant by hypothesis). 

x = the concentration of X at the time ?. 

Ky the velocity constant of the forward reaction. 

Ke the velocity constant of the reverse reaction. 

Now by the law of mass action the velocity of a chemical process is di- 
rectly proportional to the concentration of the reacting substances. ‘There- 
fore the velocity of the forward reaction is proportional to A, the concen- 
tration of the substrate, and to x, the concentration of the catalyst which 
in this case is also the end product. Likewise the velocity of the reverse 
reaction is proportional to x’, the substance X acting both as substrate 
and as catalyst.* Hence 


dx/dt = K,Ax — Kox’. (1) 


Crozier’s autocatalytic theory differs from Robertson’s in certain im- 
portant details, though the principles involved and the equation which he 
derives from them are approximately the same. He states the case as 
follows: 

“Let it be supposed that at the beginning of a developmental cycle 
there is available an unrenewable quantity, A, of a substance giving rise 
to another, x, which determines the velocity of growth. We are especially 
interested in ‘velocities of growth’ as measured by the reciprocals of the 
times required to attain a given stage of development. We will suppose 
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that the material A gives rise to x by a first order reaction, and that x 
serves as catalyst for this change. The reaction A —~> x will therefore 
be governed by a velocity constant (K,) proper to it in the absence of the 
influence of x, and also by the velocity constant due to catalysis by x. 
The decomposition of A must therefore be conceived as made up of two 
parallel reactions, and its differential equation is then 


dx/dt = (K, + Kx) (A — x) (2) 


where K¢ is the velocity constant associated with x as catalyst.”’ 

“This formulation seems simpler, and leads to more suggestive conse- 
quences so far as concerns the planning of experiments, than does Robert- 
son’s (1923) suggestion of the pseudo-reversible character of the ‘master 
process’ in growth; moreover, it would appear to be the obviously correct 
equation for an autocatalytic system.” 

Both equations (1) and (2) above are typical equations for monomolecu- 
lar autocatalytic reactions, and both describe with accuracy the appro- 
priate chemical reactions of this type as carried out in the laboratory. 
Both of them, moreover, can be made to fit the growth curve of a large 
variety of plants and animals, though Crozier’s equation is the more flexible 
in this respect and can be applied to the greater number of cases. 

There are, however, serious objections to their application to the chemi- 
cal processes involved in growth. Thus T. H. Morgan has pointed out 
that an organism, since it is constantly taking in food and giving out waste 
products, is not a closed system, and hence that to account for the retarda- 
tion and ultimate cessation of growth requires special postulates. He says: 

‘The symmetry of the upper and lower halves of the curve of growth.... 
and that of the autocatalytic curve..... cannot be used as direct evidence 
of the relationship of the two curves, for this would hold only under special 
conditions, i.e., in a closed system. In an organism that is growing and 
receiving food from outside there is no reason for supposing that the upper 
half of the curve should be the reverse of the lower half; or even that there 
should be an inflexion in the curve. Only on the assumption that the re- 
action is reversible and unimolecular in both directions can symmetry 
be expected in this case.” 

An even more fundamental objection, one that unfortunately renders 
invalid the simple theories of growth proposed by Robertson and by Crozier, 
is that the equations they use do not apply to growing organisms because 
of the effect of increasing size on the concentration of the reagents involved 
in growth. The equations hold true only on the condition that the volume 
occupied by the reacting substances remain constant, and since a growing 
organism is constantly increasing in volume this condition, of course, is 
not met. Even the condition stated by Morgan to be sufficient, namely, 
that the reaction be ‘‘reversible and unimolecular in both directions,’’ 
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is not adequate to make the customary equations for autocatalysis fit the 
processes involved in the synthesis of new protoplasm. When a chemical 
process is carried out in a laboratory the reagents are ordinarily dissolved 
in water or some similar solvent, and the volume of this solvent is kept 
constant throughout the whole process. It is to reactions carried out 
under such conditions that the usual formulas based on the mass action 
law are intended to apply. To make the conditions of a laboratory proc- 
ess comparable to those involved in the synthesis of new protoplasm the 
volume of the solvent would have to be increased as fast as the amount 
of the end products is increased. As the solids of new protoplasm are 
formed they do not stay in the same little parcel of liquid occupied by the 
old; rather they cause the liquids to expand with them. Hence the vol- 
ume occupied by the end products of growth is proportional to their 
amount, and the concentration of these products, instead of constantly in- 
creasing, remains constant. This is a very important difference, for it 
is on the concentration and not the amount of the reagents that the ve- 
locity of a reaction depends. ‘“The rate at which any reaction proceeds 
is directly proportional to the amount, or rather the concentration, of the 
reacting substances.”” (Bayliss, p. 184. Italics mine.) 

Because of this difference the customary equations for the rate of 
chemical processes cannot be applied to processes involved in growth, 
and new equations will have to be derived. These equations will involve 
a term for the volume‘ of the organism. We shall first consider the form 
of autocatalytic process which Robertson holds to be the master reaction. 

Let A = the concentration of the substrate (a constant by hypothesis). 

x the amount (not concentration) of the end product and 
catalyst at the time ¢. 

V = the size (volume) of the organism at the time ?. 

V, = the size of the organism at the beginning of the growth cycle. 

Now it is assumed by Robertson and Crozier that at the beginning of a 
new growth cycle the governing substance, X, is entirely absent. As 
the cycle progresses, however, the formation of new living tissue is limited 
by the formation of X, and hence the amount of new tissue (V — V,) 
increases in proportion to the amount of X. We have therefore 


ax=V-—YV, where ¢ is a constant, or 
V = cx +t V,. 


The concentration of X is defined as the amount of X divided by the 
volume which it occupies. Therefore 


x/V = concentration X. 


By the law of mass action the velocity of the forward reaction is pro- 
portional to the concentration, A, of the substrate and to the concentration, 
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x/V, of the catalyst. In like manner the velocity of the reverse reaction 
is proportional to x?/V?, the substance X acting both as substrate and as 
catalyst. There is, moreover, one other factor to be considered which is 
not present under ordinary conditions, namely, the effect of growth on 
the amount of the reacting substances (see quotation from Morgan). 
This will be made clear by the following consideration. If a reaction 
such as that with which we are dealing were going on in a single beaker 
we would have a given amount of X produced in a given time. If, now, 
instead of a single beaker we have two just as full as the first the amount 
of X produced in both of these together in the given time will obviously 
be doubled. In like manner the amount of X produced in the growing 
organism is doubled when the size is doubled, other things being equal, 
and the rate of production is thus proportional to V. We therefore get 
as the differential equation for the production of X 





2 
ds = K,A~ V- i. V. 
dt V Vy? 
But since V = cx + V, 
2 
os te Be bem 
dt cx + V, 


The integrated form of this equation, for certain values of the constants, 
gives a sigmoid curve just as does Robertson’s but one of decidedly different 
form. It is doubtful if it can be made to fit the majority of growth curves. 

The same principle may be applied to Crozier’s proposed type of master 
reaction. 


Let A = the initial amount of substrate. 
x = the amount of end product at the time ?. 
V, = the volume of the organism at the beginning of the growth 
cycle. 
V = cx + V, = the volume of the organism at the time ¢. 
Then 
: A —x = the amount of the substrate at the time #. 
(A — x)/V = the concentration of the substrate at the 


+ time ?¢. 
x/V = x/(cx + V,) the concentration of X at the time f. 
Hence the differential equation for the formation of X becomes 


dx x A-x\+3 
— mi K K. V 
dt ( ae +)( V ) 


x 
* A 
(m+m—4—)( ) 








This is the correct form for the equation of an autocatalyzed mono- 
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molecular reaction in which the conditions are similar to those in a growing 
organism; that is to say, in which the volume occupied by the reagents 
increases in proportion to the amount of the end product. In the in- 
tegrated form it cannot be made to fit the growth curves of the rat and 
the sunflower given by Crozier as illustrations of development controlled 
by an autocatalytic process. Moreover, in cases where the size of the 
organism at the beginning of the growth cycle is small compared to the 
size at the end, as in the embryonic cycle of mammals, the expression 
x/(cx + Vo) is almost constant and the curve approximates that of a mono- 
molecular reaction not autocatalyzed. This curve of course is altogether 
different from any growth curve. Hence the theory proposed by Crozier 
when developed so as to be in consonance with the law of mass action 
will not explain the rate at which organisms grow. 

Conclusions.—We therefore find that owing to the expanding volume 
of the growing organism the usual equations derived from the law of mass 
action cannot be applied to the processes involved in the synthesis of new 
protoplasm. Equations suited to describe the rate of such processes are 
markedly different from those cutomarily used. This fact renders invalid, 
at least in their present form, the hypotheses put forward by Robertson 
and by Crozier, according to which growth rate is controlled by a single or 
master autocatalyzed monomolecular reaction. It does not, however, 


affect Loeb’s conclusions that the production by geometric progression of 
new nuclear material in the fertilized egg is evidence for the autocatalytic 
nature of the processes involved. 
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1 Conklin did not find this to be true in Crepidula plana; rather, he found compara- 
tively little increase in nuclear material with each cell division. This, however, does not 
necessarily invalidate Loeb’s conclusions. 

? Pearl, however, asserts that the fit of Robertson’s theoretical growth curves to the 
observed data is too poor “‘to afford evidence of any particular value in favor of Robert- 
son’s ingenious...... hypothesis.” 

MacDowell, moreover, finds that when new-born mice are fed to the limit, their growth 
curve during the first fourteen days is parabolic rather than exponential. This one fact 
alone would compel a modification of the theories which Robertson and Crozier propose. 

3 As both Morgan and Crozier have pointed out, however, Robertson’s assumption 
that the master reaction is monomolecular in both directions, an assumption to which he 
is forced in order to make the autocatalytic curve coincide with that of growth, is a rather 
improbable one. Ina later paper (Robertson, 1926) a different equation is used to per- 
mit a slightly asymmetrical curve, but it does not remedy this defect, or any of the 
other defects that are here discussed. 

4 It is assumed that volume is proportional to mass throughout growth. While this 
is not strictly true, the deviation from proportionality is not sufficient to affect the 
present discussion. 


CHLOROPHYLL CONTENT AND RATE OF PHOTOSYNTHESIS 


By Ropert EMERSON* 
LABORATORY OF GENERAL PHYSIOLOGY, HARVARD UNIVERSITY 


Communicated January 25, 1929 


Two sets of factors are generally recognized as affecting the rate of photo- 
synthesis in green plants: external factors, including light, temperature, 
and carbon dioxide; and internal factors. The latter is a vague group, 
the only obvious one being chlorophyll. Owing probably to the difficulty 
of experimentally controlling internal factors, little is known about them. 
The only work leading to definite information about an internal factor 
is that of Willstatter and Stoll (1918). They investigated the photo- 
synthetic rates of a variety of leaves differing widely in chlorophyll content. 
They were able to determine chlorophyll accurately, but they could not 
control the amount of it present in the leaves. Thus in order to obtain 
leaves differing in chlorophyll content they were obliged to use leaves 
which were not at all comparable in other respects. This may be the 
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reason why they failed to find any correlation between chlorophyll content 
and rate of photosynthesis. 

The chlorophyll content of Chlorella vulgaris may be controlled over 
a fairly wide range without materially disturbing its rate of growth or res- 
piration. This is done by growing pure cultures of cells in a medium con- 
taining the appropriate salts, with about 1.5% glucose, and with the con- 
centration of iron kept below that normally used. 

The rate of photosynthesis of cells differing in chlorophyll content was 
determined at high light intensity in glass vessels attached to manometers, 
after a method described by Warburg (1924). Relative chlorophyll con- 
tent per unit volume of cells was determined by extracting the cell pig- 
__ ments (green and yellow together) 

with methyl alcohol, the extract 

| being made up toa definite volume 

| in a volumetric flask. ‘The ex- 
tinction coefficient ¢« for the wave- 
length \ = 670 yy was determined 
for each extract thus prepared. 
This wave-length was chosen be- 
cause it is outside the region of 
absorption of the yellow cell 
pigments, and near the center of 
| strong absorption by chlorophyll 

; ; | in the red. 

0005 0050 Q075 0100 Figure 1 shows rates of photo- 
£, TEMAT WE, RENAL. CREM RET synthesis plotted as ordinates 
ceieaarth against corresponding values of ¢ 


The rate of photosynthesis plotted as a : 
function of the chlorophyll content per unit as abscissae. Qos represents the 


volume of cells. number of cubic millimeters of 
oxygen produced by 10 cubic mil- 
limeters of cells in one hour at 20°C. € represents the extinction coeffi- 
cient of the extract of 10 mm.* cells made up to 25 cc. with methyl alcohol. 
The points, like those of several other similar series, lie on a smooth 
curve. The rate of photosynthesis at high light intensities is clearly a 
function of the chlorophyll content of the cells. This seems to be good 
experimental evidence that chlorophyll is really the pigment concerned 
in photosynthesis—a fact doubted by no one, but supported so far by in- 
ference only. 
_ It is a familiar fact that cutting down the rate of photosynthesis by de- 
creasing the light intensity causes a fundamental change in the character 
of the process. At low intensities of light the rate scarcely varies with 
temperature, while at high intensities it varies rapidly with temperature. 
Dilute prussic acid has practically no effect on photosynthesis at low in- 
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tensities of light, but inhibits the process considerably at high intensities. 
For these observations we are indebted to F. F. Blackman (1905) and 
O. Warburg (1919). The so-called ‘‘Blackman reaction,” a reaction with 
a high temperature coefficient and sensitive to prussic acid, is supposed 
by Warburg to be limiting the rate of the photosynthetic process at high 
light intensities, while at low light intensities a photochemical reaction, 
relatively insensitive to changes in temperature and to prussic acid, governs 
the rate. The characteristics of the two reactions are well summarized 
by Warburg (1925). If this interpretation is correct, then it should be 
possible to separate the Blackman reaction from the photochemical reac- 
tion either by reducing the light 
intensity, or by cutting down the 
chlorophyll content. The result A 
in either case should be a decrease 
in the rate of the photochemical 
reaction until it is slower than the 
Blackman reaction, and becomes 
the “‘pacemaker” of the process. 
Figure 2 shows two curves rep- 
resenting rate of photosynthesis 
as a function of temperature at 
widely different chlorophyll con- 
centrations. For curve A, ¢€ = i 
0.101, and for curve B, e = 0.016. 
The data for most processes, when 
plotted in this way (logarithm of rt an 
the rate or frequency against re- aes IBSOLUTE 
ciprocal of the absolute temper- FIGURE 2 
ature), usually show a straight The rate of photosynthesis plotted as a 


line relationship. ‘The culiar function of temperature, at two different 
P ” chlorophyll concentrations. For curve A, 


shape of the curves presented will ~ _ 9 101. for curve B, e = 0.016. 

be discussed in a later paper. I 

wish here to bring out the point that these curves are similar in form, 
showing that photosynthesis is the same function of temperature whether 
the chlorophyll concentration be high or low. 

Table 1 shows the inhibition of photosynthesis by prussic acid at three 
different concentrations of chlorophyll. There is no tendency for the 
inhibition to be less at lower chlorophyll concentrations. On the con- 
trary, at the lowest chlorophyll concentration the inhibition is 15 per cent 
more than at the highest. 

Although either lowering the light intensity or decreasing the chlorophyll 
content decreases the rate of photosynthesis, the characteristics of the 
process are very differently affected in each case. It is unreasonable to 
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suppose that the light intensity affects anything but the photochemical 
part of the process. ‘The amount of chlorophyll must also affect this part, 
and in the same way. However, since the characteristics of the process 
(sensitivity to temperature and to prussic acid) cannot be similarly altered 
by changing chlorophyll content and light intensity, it must be supposed 
that chlorophyll plays some other part in the process, besides its rdle in 


the photochemical reaction. 
TABLE 1 


Tue Errecr or Prussic Acip ON PHOTOSYNTHESIS AT THREE DIFFERENT 
CONCENTRATIONS OF CHLOROPHYLL 


MM.? O2 INHIBITION 
© HCN CONCENTRATION PRODUCED BY HCN 
O 
0.037 9X 10-°N 93.1 40% 
56.2 : 
O 
0.060 9X 10-*N 122.0 32% 
82.7 
O 
0.083 9X 10-*N 171.5 25% 
129.5 


A more detailed presentation of this work is in preparation, and will 
appear in the Journal of General Physiology. 
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FUNCTIONAL DISTURBANCES OF HEARING IN GUINEA PIGS 
AFTER LONG EXPOSURE TO AN INTENSE TONE 


By Morcan Upton* 
PsYCHOLOGICAL LABORATORY, HARVARD UNIVERSITY 


Communicated January 17, 1929 


I.—Present controversy regarding theories of the perception of sound 
seems to be mainly concerned with two sharply defined views. One, the 
“frequency” theory, regards the basilar membrane as a vibrating unit, 
which responds in whole or in part in a way depending upon the energy 
of the stimulating force. A stimulus of great energy would disturb a 
relatively large part of the basilar membrane and a stimulus of less energy 
would disturb a relatively shorter section beginning at the base of the 
cochlea and extending upward. Under this theory pitch is correlated with 
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frequency of nerve impulse and intensity with the number of basilar units 
set into motion by a stimulus. The other, the ‘‘place’’ theory, regards 
the transverse fibers of the basilar membrane as resonators which are set 
into motion by the sound waves which correspond to their natural periods 
of vibration, thus correlating pitch with the place of disturbance on the 
basilar membrane. The “‘place’’ theory in its modern form is derived 
from the theory of Helmholtz, and includes certain considerations made 
necessary by the establishment of the ‘‘all or none” law of nervous con- 
duction. Helmholtz assumed that the experience of intensity was corre- 
lated with the energy of the nerve impulse, but it now appears that such 
an assumption is untenable. 

The results of an experiment by Yoshii,' published in 1909, were widely 
accepted as clear evidence for the Helmholtz theory. In this experiment 
guinea pigs were’ exposed to tones of different pitch for extended periods. 
Histological examination of the cochleas of the animals so exposed revealed 
definite structural changes which were limited to transverse strips in the 
organ of Corti and occurring in the order of the frequencies of the tones 
to which the animals were exposed, with the highest frequency nearest 
the base of the cochlea. It has been shown clearly (Boring, 1926),? that 
the data obtained in Yoshii’s experiment stand just as firmly in support 
of a “frequency” theory of hearing. Both the amplitude and the fre- 
quency of the stimuli used by Yoshii varied and either one might have 
determined the locus of degeneration. 

Although Yoshii’s work demonstrates that certain definite structural 
changes in the auditory mechanism accompany long exposure to intense 
tones, without functional differentiation his results must occupy an equiv- 
ocal position as regards the two types of theories of hearing. The 
present experiment was designed with a view to the determination of the 
functional changes correlated with extended exposure to intense tonal 
stimuli. 

II.—The effectiveness of relatively pure tones as stimuli in the environ- 
ment of guinea pigs was demonstrated in a preliminary experiment by the 
establishment of a specific change in the breathing rhythm as a condi- 
tioned response to the sounding of a tone. The unconditional stimulus 
was an electric shock applied to the rear legs of the animals. A tuning- 
fork oscillator of 600 cycles driving a loud speaker was the source of the 
tone. The tone presentation lasted through six respiration cycles and the 
shock was applied at the cessation of the tone. After 250 presentations 
of the tone and shock in combination, evidence of conditioning was ob- 
servable in the breathing curve as recorded on a kymograph. ‘The respira- 
tion cycles coincidental with the tone and preceding the shock were equal 
in amplitude, whereas the breathing curve not coincidental with the tone 
was characterized by wide and constant variation of amplitude. After 
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approximately 500 presentations of the tone-shock combination a respira- 
tion cycle of very wide amplitude occurred in the breathing curve coinci- 
dental with or slightly later than the end of the tone when the tone was 
presented without application of the shock. ‘These changes were regarded 
as clear evidence of the existence in the guinea pig of a mechanism which 
mediates the reception of pure tones. ‘Tonal differentiation in the guinea 
pig was demonstrated for tones whose vibration rates differed by 400 
cycles. Finer differentiation is very likely possible, but this was close 
enough for the purpose of the present experiment. 

Upon completion of the preliminary experiment two groups of animals 
were exposed to a very intense tone with a frequency of 600 cycles, for 
a period of seventy days. ‘Tuning-fork oscillators and a loud speaker 
served as the source of the tones employed in this experiment. Distances 
of 24 inches and 96 inches, respectively, separated the two groups from 
the source of the tone so that each group would be affected by a different 
intensity. 

At the termination of the exposure period the animals were tested for 
direct responses (i.e., change in respiratory rhythm) to three intensities 
(low, intermediate and high) of two frequencies (600 cycles and 1000 
cycles). Corresponding intensities of the two frequencies were equal. 

The low test tone with a frequency of 600 cycles had been used as the 
original conditioning tone while the high test tone of the same frequency 
was identical with the tone to which the animals had been exposed. A 
group of normal animals was introduced as a control. The responses 
of the members of this group to corresponding intensities of both frequen- 
cies were uniformly similar. 

A comparison of the responses of the exposed groups with the responses 
of the control group revealed significant deviations from the normal. 
Group I (exposed with a 96-inch separation from the source of the tone) 
gave definite responses to all three intensities of the exposure frequency. 
Group II (exposed with a 24-inch separation from the source of the tone) 
gave no responses to any of the three intensities of the exposure frequency. 
Two of the animals of this group had previously been conditioned to re- 
spond to the low intensity of the 600 frequency. An attempt to recon- 
dition them to the same tone was unsuccessful. 

The sensitivity of the animals in group I to tones of the exposure fre- 
quency had evidently been increased, while the sensitivity of the animals 
in group II to the exposure frequency had been practically eliminated. 
The control group gave no responses to the lowest intensity of both fre- 
quencies. The responses of both groups I and II to the three intensities 
of the 1000 frequency were uniformly similar to the responses of the con- 
trol group to the corresponding intensities of both frequencies. 

These results indicate that exposure to a particular frequency is accom- 
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panied by certain specific changes: first, by an increase of sensitivity to 
the exposure frequency, and then by a practical desensitization to all in- 
tensities of the exposure frequency within a very wide range. ‘The simi- 
larity of the responses of all groups to the frequency of 1000 cycles indicates 
that no general change in the auditory mechanism accompanied exposure 
to a particular frequency. 

III.—From the results of the experiment it may be said that: (1) ex- 
posure to an intense tone over an extended period of time is accompanied 
by a total loss of sensitivity to tones of the exposure frequency within a 
very wide range of intensities; (2) exposure to an intense tone is not, ap- 
parently, accompanied by any change in the sensitivity of the exposed 
animals to other frequencies. 

The results stand as evidence for a ‘‘place-pitch” or ‘‘resonance’’ theory 
of hearing. Such a theory would regard the cochlea as an organ of pitch 
discrimination. The discrimination would depend upon the vibration of 
specifically localized transverse strips of the organ of Corti which would 
be set into vibration by the sound waves corresponding to their natural 
periods of vibration. Yoshii (1909) found that the width of the strips 
in the organ of Corti which were damaged by exposure to intense tones 
varied with the amplitude of the tones. This linear spread over the basilar 
membrane, with its center at the fibre whose natural period of vibration 
corresponds to the stimulating frequency, offers an explanation of the 
perception of intensity. Each increment of energy would result in an 
increase of amplitude of the stimulus and therefore in a wider spread along 
the basilar membrane, thus including additional nerve fibres. Upon this 
assumption the perceived intensity would be a direct correlate of the 
number of fibres involved in the disturbed section of the basilar mem- 
brane. 


* NATIONAL RESEARCH COUNCIL FELLOW, Laboratory of General Physiology, Har- 
vard University, Cambridge, Massachusetts. 

1 Yoshii, U. “Experimentelle Untersuchungen iiber die Schadigung des Gehérorgans 
durch Schalleinwirkung,” Zeit. Ohrenheilkunde, 1909, 58, 201-250. 

2 Boring, E.G. ‘‘Auditory Theory with Special Reference to Intensity, Volume and 
Localization,’’ Amer. J. Psychology, 1926, 37, 157-188. 

















288 PALEONTOLOGY: T. W. VAUGHAN Proc. N. A. S. 


STUDIES OF ORBITOIDAL FORAMINIFERA: THE SUBGENUS 
POLYLEPIDINA OF LEPIDOCYCLINA AND ORBITOCYCLINA, 
A NEW GENUS 


By THomAs WAYLAND VAUGHAN 


Scripps INSTITUTION OF OCEANOGRAPHY, UNIVERSITY OF CALIFORNIA, 
LA JOLLA, CALIFORNIA 


Communicated February 18, 1929 


Introduction.—In 1924 I proposed (Vaughan, 1, p. 807) the name Polylep- 
idina for a subgenus of Lepidocyclina with Lepidocyclina chiapasensis 
Vaughan as its holotype. The diagnostic features of the subgenus were 
based on the embryonic chambers of the megalospheric form and on the 
equatorial chambers. It was stated that the megalospheric embryonic 
apparatus is composed of four or five chambers, one or two of which are 
larger than the others, and that in some specimens there are two subequal 
chambers somewhat larger than two other subequal chambers, the four 
chambers so arranged as to form a cross. The equatorial chambers were 
described as having arcuate outer walls, pointed or truncate inner ends, 
and greater transverse than radial diameters. 

Two other species, referred to Polylepidina, L. adkinsi and L. protet- 
formis, were described at the same time, but it was pointed out in dis- 
cussing L. adkinsi that in some specimens of it the embryonic apparatus 
is very similar to that of specimens of Lepidocyclina (Lepidocyclina) in 
which accessory embryonic chambers lie adjacent to the two larger cham- 
bers. For this and other reasons, which I shall not elaborate here, it is 
my opinion that Polylepidina is not valid as a genus. Therefore, I do not 
concur in a suggestion in a paper by Prof. J. J. Galloway mentioned below. 
The species that I am placing in Polylepidina represent an interesting 
group of forms in that they appear to throw light on the possible phylogeny 
of the genus Lepidocyclina and perhaps even on that of the entire family 
of the Orbitoididae. 

Mrs. Helen K. Hodson in 1926 (pp. 23-28, pl. 7, Figs. 1-8) described 
two species of Lepidocyclina, L. zuliana and L. einiealianes; which she re- 
ferred to Polylepidina, and she doubtfully referred to it another species, 
which she designated L. churuguaritana. 

Prof. J. J. Galloway in 1928 proposed to elevate Polylepidina to generic 
rank and added to my original three species L. antillea Cushman and L. 
mortont Cushman and described as new a species from the Upper Cretaceous 
on the east side of Cardenas, San Luis Potosi, Mexico, assigning to it the 
name Polylepidina cardenasensis. My doubt as to the generic validity 
of Polylepidina has been stated above. Lepidocyclina antillea is referable 
to Polylepidina, but L. mortoni, in my opinion, is not. Professor Gallo- 














VoL. 15, 1929 PALEONTOLOGY: T. W. VAUGHAN 289 


way’s Polylepidina cardenasensis is a synonym of Lepidorbitoides minima 
H. Douvillé, which appears to represent a new genus to which I am apply- 
ing the name Orbitocyclina. 

In 1928 I published (Vaughan 2, p. 350, also p. 349) additional notes 
on Polylepidina. 

Supplementary Notes on the Original Species of Polylepidina.—Recently 
I have restudied the embryonic and equatorial chambers of each of the 
original species of Polylepidina. No significant new information on the 
embryonic chambers was obtained, but with better thin sections it is possi- 
ble to present more satisfactory illustrations (pl. , Figs. 3, 4, 5). 

The equatorial chambers need more detailed consideration. There is 
a pronounced tendency for them to be arranged in radial rows. As the 
sides of the walls of chambers in the same circle usually do not meet each 
other, rather regularly the walls of chambers in adjacent rows alternately 
overlap. These features are shown on plate, figures 7 and 8. 

The structure of the walls of the equatorial chambers is complex, that 
of the roofs and floors being markedly different from that of the parts of 
the walls between adjacent chambers of the equatorial layer. The roofs 
and floors will be described first. 

The roofs and floors in a section cut vertically to their upper and lower 
surfaces exhibit a columnar perforate structure. The perforations are the 
means for protoplasmic communication between the equatorial and lateral 
chambers, and similar perforations in the roofs of the lateral chambers 
afford means for communication between the successive lateral chambers. 
These perforations may be more advantageously examined in sections 
parallel to the equatorial plane and are illustrated in the upper part of 
figure 8 on plate. Their diameter ranges between 5 and 7y. 

The distal part of the wall of each equatorial chamber consists of three 
parts. Exteriorly there are pectinations which produce an appearance 
of perforations, especially as this part of the wall is or appears to be con- 
tinuous with the roofs and floors, but in many instances at least the spaces 
between the tooth-like processes do not reach the dark line next to be 
described. 

Immediately interior to the pectinate part of the wall is a dark zone as 
seen in both vertical and horizontal section (see pl., Figs. 7 and 8). The 
nature of this part of the wall has not been adequately studied. Its color 
may be due to the carbonaceous residue of organic matter. A point of 
interest in connection with it is that no indication of any canal system 
could be found. 

Interior to the dark zone above described is a very minutely fibrous 
zone, the position of which is indicated by plate, figures 7 and 8, but the mag- 
nification is insufficient to show its structure. No perforations could be 
detected in this layer which lines the distal and lateral parts of the cham- 
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DESCRIPTION OF PLATE 


FIGURES 1, 2 


Lepidocyclina (Lepidocyclina) mortoni Cushman, x about 33. Figure 1, embryonic 
chambers of a megalospheric specimen. Figure 2, vertical section, the equatorial layer 
extends vertically with the distal part upward, to show the apertures through which 
..Jlons connected the protoplasmic content of adjacent chambers. Two such apertures 
show as circular white dots to the right of the letter ‘‘a’’ on the figure. The specimens 
figured are topotypes, collected by T. W. Vaughan in the Jackson formation at Mont- 
gomery, Louisiana. Coll. Scripps Inst. of Oceanography. 


FIGURE 3 


Lepidocyclina (Polylepidina) chiapasensis Vaughan, x about 28. To illustrate the 
embryonic chambers of a megalospheric specimen cut in nearly the equatorial plane. 
Three larger chambers of nearly equal size, two of which are crossed by cracks, and 
five smaller chambers are shown in the figure. The specimen is a topotype collected 
by M. F. Nesbit, 3.75 km. south and 6.25 east of Chil6én, Chiapas, Mexico, Coll. 
U. S. Nat. Museum. 


FIGURE 4 


Lepidocyclina (Polylepidina) adkinsi Vaughan, x about 27. Part of an equatorial 
section of a megalospheric specimen. To show the embryonic apparatus which con- 
sists of two larger subequal chambers and smaller accessory chambers. Specimen, 
a topotype, collected by M. F. Nesbit, 2.75 km. south of Hacienda ‘El Triunfo,”’ 
Chiapas, Mexico. Coll. U. S. Nat. Museum. 


FIGURES 5, 6,7 


Lepidocyclina (Polylepidina) antillea Cushman, x about 85. Figure 5, section, nearly 
in the equatorial plane, to show the embryonic chambers of a megalospheric specimen. 
Figure 6, section, approximately vertical, to show the embryonic chambers of another 
megalospheric specimen. Specimens illustrated by figures 5, 6, collected by T. W. 
Vaughan, U. S. G. S. locality no. 6902, near Négre Point, Island of St. Bartholomew. 
Coll. U. S. Nat. Museum. Figure 7, equatorial section of a microspheric specimen. 
Topotype, collected by T. W. Vaughan, U. S. G. S. locality no. 6897, point between 
Anse Ecaille and Anse Lézard, Island of St. Bartholomew. Coll. Scripps Institution. 


FIGURE 8 


Lepidocyclina (Polylepidina) proteiformis Vaughan, x about 80. Section in the 
equatorial plane of a megalospheric specimen to show the embryonic chambers and the 
surrounding equatorial chambers. The embryonic apparatus consists of a larger 
chamber entirely surrounded by smaller chambers, the uppermost chamber of the ring 
as shown in the figure is about twice the size of those that extend from its ends to com- 
plete the ring. The arrangement of the chambers is similar to that in L. chiapasensis 
Vaughan, figure 3 of this plate. Specimen illustrated, a topotype, collected by T. W. 
Vaughan, M. 118 V., about 0.5 km. southwest of Palma Sola, Mecapala Hills, Vera 
Cruz, Mexico. Coll. Scripps Institution. 

(The photographs for the illustrations were made by Mr. E. M. Thorp at the Scripps 
Institution.) 
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ber. It is possible to examine the inside of the walls of many chambers, 
especially of those which show the stoloniferous apertures described be- 
low. There are no perforations in the sides of the equatorial chambers 
such as there are in their roofs and floors. Under a magnification of 450 
diameter the texture is minutely granular in appearance, perhaps, but 
not necessarily, because of recrystallization. The inner surfaces of the 
chamber walls were examined in each of the three original species of Polylep- 
idina but the material of L. (Polylepidina) proteiformis was the best. L. 
mortont furnished excellent material for such a study and gave results 
accordant with what is said above. 

Communication between the equatorial chambers of each of the species 
now under consideration is by means of stolons for the passage of which 
there are apertures in the chamber walls. I published notes on these aper- 
tures in L. chiapasensis and L. proteiformis in 1928 (Vaughan, 2, p. 349). 
In the former species the range in diameter is from 30 to 38u and there 
may be three in a row in a chamber 210, tall. In the latter species the 
range in diameter is from 23 to 30u, and there is one in a chamber 110y 
tall and 2 in a chamber 160, tall. These apertures are well shown in one 
of the original illustrations of L. proteiformis (Vaughan, 1, pl. 32, Fig. 7) 
where they appear as circular white dots in the chambers at and near 
each end of the figure. In L.-adkinsi the diameter of the apertures is be- 
tween 20 and 25u. Three are shown in a chamber 120y tall. Each 
equatorial chamber seems to communicate with four adjacent chambers, 
two of which are in the next inner and two in the next outer circle of cham- 
bers. 

Notes on Lepidocyclina (Lepidocyclina) mortont Cushman. (P1., Figs. 
1, 2).—Professor Galloway in his paper cited above (p. 301) states after 
the name of this species ‘“Upper Eocene, Ocala of Alabama and Jackson 
of Mississippi.”’ ‘The type locality of L. mortoni is Montgomery, Louisiana, 
in the exposure there of the Jackson formation, where the type specimens 
were collected by me. Professor Galloway on p. 300 of the same paper 
figures specimens from 5 mi. north of Geneva, Alabama, and he refers 
them to P. mortoni. I am confident that he made an error in ‘identifica- 
tion, because the embryonic chambers of the megalospheric form and the 
equatorial chambers as figured by him do not accord with those of L. 
mortoni. (Cushman, 2, pl. 27, Figs. 14; Vaughan, 1, pl. 36, Figs. 1, 2.) 
In order to make this clear a figure of the embryonic apparatus of a topo- 
type, collected by me, is here given (pl., Fig. 1). Comparison of it with 
Professor Galloway’s figure shows that he was dealing with a different 
species. L. mortoni belongs to the subgenus Lepidocyclina. A second 
figure of L. mortoni (pl., Fig. 2) represents a vertical section and illustrates 
the stoloniferous apertures between chambers in the equatorial layer. 

Orbitocyclina Vaughan, ‘New Genus.—A few words will be said about the 
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species to which Prof. H. Douvillé (2, 1927) has applied the name Lepidor- 
bitoides minima and Professor Galloway (1928) the name Polylepidina 
cardenasensis, from the Upper Cretaceous at Cardenas, San Luis Potosi, 
Mexico. Although after much study I concurred with Professor Douvillé 
in his assignment of this species to Lepidorbitoides and have so expressed 
myself in an article submitted for publication in the Journal of Paleontology; 
still further study has resulted in the discovery of stoloniferous apertures 
in the lateral walls of the equatorial chambers (see H. Douvillé, 1, 1923). 
The species, therefore, does not belong to the genus Lepidorbitoides. Since 
the embryonic apparatus of the megalospheric form differs from that of 
Polylepidina and the other subgenera of Lepidocyclina it is not referable 
to Lepidocyclina. The equatorial layer differs from that of Pseudorbi- 
toides, the genus to which Professor Douvillé (1923) originally referred 
the species. Under the circumstances the only course to take is to pro- 
pose a new generic designation. Therefore, I propose the generic name 
Orbitocyclina with Lepidorbitoides minima H. Douvillé as defined by me 
in my article above mentioned as the genoholotype. The first two cham- 
bers of the embryonic apparatus of the megalospheric form are similar 
to those of Lepidorbitoides, that is reniform, except that they are rather 
smaller. The second chamber is followed by three or four other chambers 
which extend around the initial chamber and entirely or nearly envelop 
it. The equatorial chambers are similar to those of Polylepidina. ‘Two 
apertures, ranging from 7 to 8.5u in diameter, one near the roof and one 
near the floor, were observed on the inner sides of several chambers. 
Larger apertures, as much as 17y in diameter, were noticed in the anterior 
parts of at least two chambers. The inner surface of the walls in those 
chambers in which they could be examined exhibited no perforations ex- 
cept those for the stoloniferous apertures. 

Lepidocychina (Polyepidina) antillea Cushman (Plate, Figs. 5, 6, 7).— 
Lepidocyclina antillea Cushman, 1, p. 63, pl. 3, figure 3, and 2, p. 78, pl. 35, 
figures 4, 5, 1920; Vaughan, 1, p. 800, pl. 35, figures 1, 2, 3, Polylepidina 
antillea Galloway, p. 300, 1928. 

Cushman’s original description of the general features of Lepidocyclina 
antillea is good and as it is easily accessible it will not be repeated. In 
my paper, 1924, the equatorial chambers were described and figured in 
more detail. I pointed out that in its initial stages the microspheric form 
is spiral and figured a specimen; but I was not able to describe and illus- 
trate the embryonic chambers of the megalospheric form. Additional 
sections recently cut for me at the United States Geological Survey ren- 
der it possible to ascertain the general features of the embryonic chambers. 

The embryonic apparatus of the megalospheric form as shown in two 
sections consists of at least four chambers, one of which is somewhat larger 
than those that follow it; and the chambers are arranged in an indistinct 
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partial spiral. The clearer of the two sections is illustrated by plate, figure 
5. It is probable that a complete picture of the embryonic apparatus 
has not been procured, but enough is now known to show that it conforms 
to the Polylepidina type. 

The first and largest chamber in the best section is somewhat irregular 
in shape, flattened on the side adjacent to the next formed chamber. Its 
greatest diameter is about 2254. The next chamber is subelliptical, 
flattened on the side adjacent to the first chamber, and 150 by 110y in 
dimensions. The third chamber has an arcuate outer wall and on its 
inner side overlaps the dividing wall between the first two chambers. Its 
dimensions are 190 by 1304. The fourth chamber has an arcuate outer 
wall which overlaps from the third to the first chamber, and its dimensions 
are 210 by 1504. There is in this specimen a small fifth chamber which 
overlaps the boundary between the third and fourth chambers. Its outer 
wall is arcuate and its dimensions are 170 by 1004. The chambers are 
not all in the same plane. The fourth chamber in the illustrated embryonic 
apparatus occupies a somewhat higher position, as the section was viewed, 
than the preceding chambers. An obliquely vertical section of another 
specimen shows a larger chamber, on one side of which there are two ad- 
jacent apparently smaller chambers, one above the other (pl., Fig. 6). The 
embryonic apparatus is composed of a so-called ‘mulberry mass.”’ 

The walls of the embryonic chambers are perforated in a complicated 
way. ‘They are composed of elements approximately perpendicular to 
the outer surface of the chamber wall with interspaces which are crossed 
by elements approximately parallel to the outer surface. The structure 
is irregularly trabecular. The thickness of the walls is variable, ranging 
from about 16 to 35yu. The thickness of the vertical elements ranges 
between about 2 and 4y (measurements only approximate); the width 
of the interspaces is about the same or somewhat less. The structure 
of the outer wall of chamber no. 4 of the figure is distinct and is somewhat 
coarser than that indicated by the measurements just given. 

That the microspheric form of L. antillea has a spirally coiled nucleoconch 
was clearly shown in an illustration I published in 1924, but the specimen 
figured does not reveal the detail of the early chambers. A subsequent 
preparation made by Mrs. E. G. Moberg renders their study possible. 

The boundary of the actual primary chamber is cloudy on one side, 
but its diameter is about 34y. It appears to be followed by a somewhat 
smaller chamber, about 17y in diameter, from which it seems to be di- 
vided by a curved wall, the convexity of which is toward the center of the 
primary chamber. Since the boundary above indicated is vague, it may 
not be real, and the primary chamber may be pear-shaped and 34 by 50yu 
in dimensions. ‘The initial chamber (or chambers) is followed by a coil 
composed of seven chambers. The primary coil is followed by a second 
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coil of nine chambers, and these are succeeded by six chambers which 
form only a part of a coil. There are about two and a half coils, com- 
posed of about 21 chambers, not counting the initial chamber (or chambers). 
The total diameter of the coiled part of the test is about 350u. No defi- 
nite structure of the walls in the chambers of the coil was recognizable. 
There was no evidence of any canal system. The outer chambers of the 
coils pass directly into the equatorial chambers without any sharp dis- 
tinction. 

The equatorial chambers have arcuate outer walls, the amount of the 
curvature varying considerably, but it is rarely or never so much as to 
make the radial diameter as great as the transverse. The size of the cham- 
bers increases from the center toward the periphery. Near the center, 
radial diameter about 504; transverse, about 100u. Near the periphery 
radial diameter, 100 to 120u; transverse diameter, 200 to 250u. The 
transverse is usually about twice the radial diameter. In vertical sec- 
tions the outer walls are nearly vertical; their height near the center, 
about 50 at the periphery as much as 300u. 

The equatorial chambers are arranged in radiating rows and those in 
adjacent rows alternately overlap. Each chamber communicates with 
four other chambers, on each side with a chamber of the next inner circle 
of each adjacent row and with a chamber of the next outer circle of each 
adjacent row. That is, of the four chambers with which any chamber 
communicates two belong to the next inner and two belong to the next 
outer circle. 

The apertures through which adjacent chambers communicate are clear 
in very few chambers of the vertical sections. There may be as many 
as three in a row, with a diameter of about 254. They are clear in some 
of the horizontal sections and are from 15 to 18 in diameter. Except 
these apertures there are no perforations between the equatorial chambers 

The chamber roofs and floors are pierced by cribiform perforations 
which range between 3.5 and 7y in diameter. These perforations are 
excellently shown in a section that passes through or just above the cham- 
ber roofs. They occur above the chamber cavities, but not in the walls 
between the cavities. The structure of the roofs and floors is trabecular, 
similar to that of the walls of the embryonic chambers described in a pre- 
ceding paragraph. 

The walls on the distal sides of the chambers, as seen in both vertical 
and horizontal sections, are similar in structure to those of the species 
already considered in this paper. Communication between chambers of the 
equatorial layer is confined to the stoloniferous apertures described above. 

Localities and Geologic Horizon At numerous places in the Island of 
St. Bartholomew, French West Indies, in the upper Eocene St. Bartholo- 
mew limestone. 
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Conclusion.—The orbitoidal foraminifera discussed in this paper are 
of interest both from the standpoint of stratigraphic geology and that 
of the possible phylogenetic relations of the group of organisms to which 
they belong. The new genus to which I am applying the name Orbito- 
cyclina occurs in the Upper Cretaceous of eastern Mexico. Polylepidina 
is known only from deposits of Eocene age. Both Orbitocyclina and Polylepi- 
dina are very close to Orbitoides, from which they differ in the features 
of the embryonic chambers of the megalospheric forms. In Orbitoides 
the entire embryonic apparatus is enveloped in a thick shell, within which 
are several chambers. The embryonic apparatus of Polylepidina is what 
Hofker (1927) has designated as a ‘“‘mulberry mass’’ and is very different 
from that of the Camerindae (‘‘Nummulitidae’’). Furthermore there 
is no canal system in the forms here considered. The ancestry of the orbi- 
toids must be sought in a phylum probably represented by Planorbulina, 
or some allied form, as Hofker has suggested. 
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